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ABSTRACT 


DIGITAL  CONTROLLER  DESIGN  FOR 
ALIASED  SAMPLED-DATA  SYSTEMS 

David  S.  Hansen 

Department  of  Mechanical  Engineering 
Doctor  of  Philosophy 

Many  of  the  dynamic  systems  in  use  today  consist  of  an  analog  plant  controlled  by  a 
digital  controller.  Design  of  these  digital  controllers  assumes  a  system  with  no  aliasing, 
and  most,  if  not  all,  current  design  methods  employ  means  to  assure  no  aliasing  is 
present.  The  drawbacks  to  these  methods  are  unobservable  and/or  uncontrollable  high 
frequency  modes. 

This  research  effort  provides  a  means  to  design  a  digital  controller  for  systems  with 
aliased  modes.  A  model  is  developed  in  which  both  continuous  and  discrete  portions  of 
sampled-data  system  are  included.  From  this  model,  the  unique  digital  controller 
function  that  will  excite  the  analog  plant  is  found.  Moreover,  as  a  result  of  the  discrete- 
to-continuous  manner  of  the  model,  the  sample  period  allowing  aliasing  becomes  a 


design  parameter.  Sampling  frequencies  providing  improved  response  are  determined 
and  responses  at  these  frequencies  are  explored. 

With  aliasing  allowed  in  the  discrete-to-continuous  system,  stability  for  some  systems  is 
actually  introduced  and  improved  over  the  same  systems  in  which  aliasing  is  not  allowed, 
as  evidenced  by  increased  gain  margins  and  improved  time  domain  responses. 

A  deadbeat  response  controller  is  designed  using  proportional-integral  control  by  placing 
one  pole  of  a  second-order  system  at  the  optimal  position  of  the  origin  in  the  z-plane. 

The  response  of  this  controller  is  such  that  the  desired  position  is  obtained  within  a  single 
sample  period.  Robustness  of  this  controller  is  examined  with  favorable  results. 

In  a  fourth-order  system  in  which  one  mode  is  aliased,  control  of  the  aliased  mode  is 
accomplished  simply  by  selecting  the  sample  rate  correctly.  The  resulting  effects  of  the 
aliased  mode  are  negligible,  and  design  of  the  controller  for  the  non-aliased  mode  follows 
traditional  theory. 

Disturbance  rejection  for  both  the  second-order  and  fourth-order  systems  was  examined 
and  found  to  be  excellent  for  disturbances  with  frequencies  below  that  of  the  sampler. 
These  systems  were  unable  to  follow  input  or  reject  disturbances  with  frequencies  greater 
than  that  of  the  sampler.  Experiments  on  electric  circuits  matched  the  predicted  model. 
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Chapter  One.  Introduction 


1.0  Background 

As  computers  have  become  the  norm  in  use  for  control  of  analog  systems,  much  research 
has  been  conducted  into  methods  of  digital  control  of  these  systems.  Many  of  the 
dynamic  systems  in  which  control  is  utilized  consist  of  an  analog  plant  controlled  by  a 
digital  controller  (computer). 

Several  methods  are  currently  employed  to  design  digital  controllers  [3].  All  these 
methods,  however,  assume  a  system  with  no  aliasing,  ignore  the  continuous  response 
between  samples,  and  ignore  the  computational  delay — each  assuming  one  can  sample 
and  compute  as  fast  as  is  necessary.  A  brief  review  of  these  design  methods  will  follow. 
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1.1  Current  Digital  Controller  Design 

One  of  the  easiest  design  methods  to  understand  is  that  of  emulation.  Emulation  involves 
designing  the  controller  in  the  continuous  (s-plane)  domain  using  traditional  root  locus  or 
frequency  response  design  methods,  and  then  mapping  the  controller  design  from  the 
continuous  domain  into  the  discrete  (z-plane)  domain.  Several  methods  for  the  mapping 
have  been  developed.  Included  in  these  methods  are  numerical  integration  methods,  the 
zero-pole  matching  method,  and  hold  equivalent  methods  [3:187-208, 215-222]]. 

Another  method  of  designing  digital  controllers  involves  direct  design  by  the  root  locus 
methods  in  the  z-plane.  Performance  specifications  familiar  in  the  s-plane  are  translated 
into  the  z-plane,  and  the  mechanics  of  root  locus  plotting  introduced  by  Evans  remain  the 
same.  This  method  requires  a  complete  discretization  of  the  model  prior  to  any  design 
[3:222-234]. 

Design  in  the  digital  domain  is  also  accomplished  using  frequency  response  function 
methods.  However,  the  effects  of  the  sampling  and  zero-order  hold  make  a  computer 
necessary  for  the  design.  The  simple  hand-plotting  techniques  for  continuous  systems  no 
longer  hold  true  in  the  discrete  domain.  Interpretations  of  phase  and  gain  margins  are 
also  more  difficult,  causing  the  need  for  more  analysis  [3:234-264]. 
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Both  emulation  and  direct  digital  domain  design  can  be  used  in  state  space  design,  as 
well.  Additionally,  other  methods  such  as  LQR  optimal  control,  Kalman  filters,  and 
Luenberger  observers  are  powerful  tools  in  the  design  of  digital  controllers  using  state 
space  design  methods  [3:270-351]. 

As  mentioned  above,  all  these  digital  design  methods  make  assumptions  as  to  the 
necessary  sample  rate,  as  well  as  the  time  between  sampling,  computation,  and  control 
output.  Sample  rates  of  20  to  30  times  the  system  natural  frequencies  must  be  employed 
in  most  cases.  Moreover,  computation  time  is  generally  assumed  to  be  as  fast  as 
necessary  in  order  for  it  to  be  neglected  in  the  design  of  the  controller.  Physical  systems, 
however,  may  make  such  assumptions  invalid. 


1.2  Aliasing 

In  sampled-data  systems,  where  the  plant  is  continuous  and  the  controller  is  digital,  a 
digitization  or  sampling  of  the  analog  plant  must  be  performed.  One  well-observed  effect 
of  sampling  analog  data  is  the  ambiguity  that  occurs  when  a  signal  of  high  frequency  is 
sampled  at  a  low  frequency.  Shannon’s  well-know  sampling  theorem  [1]  states  that  any 
signal  whose  frequency  components  are  greater  than  one-half  the  sampling  frequency, 
known  as  the  folding  or  Nyquist  frequency  [2],  will  be  aliased  to  a  lower  frequency. 
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These  aliased  frequencies  will  become  indistinguishable  from  the  system’s  response  in 
the  lower  frequency  range. 

Aliasing  causes  errors  in  the  characterization  of  the  system  or  plant.  These  errors  can  be 
severe  if  the  sampled  signal  contains  enough  high-frequency  components.  Because  of  the 
sampling  operation,  the  signal’s  high  frequency  components  will  not  be  correctly 
identified.  Li  fact,  the  high-frequency  content  will  incorrectly  appear  as  lower-frequency 
components  in  the  sampled  signal.  A  simple  example  will  demonstrate  this: 

Suppose  one  wishes  to  sample  a  10  Hertz  signal.  If  samples  are  taken  every  0.0833 
seconds  (fs  =  12  Hz)  the  sampled  signal  will  appear  as  a  sine  wave  of  only  2  Hertz,  as 
shown  in  figure  2.2.  The  10  Hertz  signal  has  been  aliased  to  a  2  Hertz  signal  by  the 
sampling  operation.  This  is  an  example  of  aliasing  of  a  pure  sinusoidal  signal.  Aliasing 
of  a  pure  sinusoid  can  also  be  thought  of  as  a  “folding”  operation.  The  explanation  of  this 
follows. 


Figliola  and  Beasley  [18:273]  define  the  Nyquist  frequency  as 


Figure  1.1  Aliasing  of  10  Hz  Signal 


Any  content  in  the  signal  above  the  Nyquist  frequency  will  be  aliased  and  will  falsely 
appear  as  frequency  components  of  the  signal  below  the  Nyquist  frequency.  The  Nyquist 
frequency  is  also  known  as  the  folding  frequency  since,  in  the  frequency  domain,  aliased 
frequencies  can  be  thought  of  as  having  been  folded  back  at  that  frequency  and 
superimposed  on  the  signal  at  lower  frequencies.  This  is  illustrated  for  the  simple 
example  given  above  in  figure  1.2. 
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Figure  1.2  Aliasing  Shown  as  Folding  in  Frequency  Domain 


Recall  the  sampling  frequency  was  12  Hz,  so  the  folding  frequency  is  6  Hz.  The 
frequency  band  from  dc  (f  =  0)  to  the  folding  frequency  (f  =  f^2)  is  known  as  the  primary 
band.  Any  frequency  components  above  the  Nyquist  frequency  will  fold  as  many  times 
as  necessary  to  be  within  the  primary  band.  Thornhill  &  Smith  [19:7.7-7.1 1]  and  Figliola 
&  Beasley  [18:273-275]  offer  excellent  examples  of  this  folding  technique  and  its  aid  in 
visualizing  aliasing. 

Likewise,  aliasing  of  signals  not  purely  sinusoidal  can  be  thought  of  as  this  same  folding 
operation.  The  component  of  the  signal  found  above  the  folding  frequency  is  folded  back 
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to  its  corresponding  frequency  within  the  primary  band  and  added  to  the  components  of 
frequency  already  present  there.  Based  on  the  band  from  which  the  frequency  component 
is  folded,  a  complex  conjugation  might  also  be  necessary.  This  folding  operation  is 
linear,  and  the  aliased  FRF,  given  the  original  FRF  is  linear,  remains  linear.  Dumey 
offers  an  excellent  analysis  of  this  operation  in  his  thesis  [6:15-18]. 

One  current  practice  used  to  avoid  aliasing  is  to  precede  the  sampling  operation  with  low 
pass  filters.  This  will  remove  the  spectral  content  above  the  Nyquist  frequency,  rendering 
these  modes  unobservable.  Another  practice  is  to  smooth  the  output  of  the  controller 
such  that  the  high  frequency  components  are  not  excited.  One  drawback  to  these 
methods,  however,  is  that  these  modes  thus  become  uncontrollable. 

One  may  also  avoid  aliasing  errors  by  sampling  at  a  rate  greater  than  twice  the  largest 
frequency  content  present  in  the  system  response.  This,  however,  requires  a  greater 
knowledge  of  the  system  to  be  controlled,  and  may  require  expensive  hardware  to  sample 
at  the  required  frequencies.  In  other  cases  it  may  simply  be  impossible  or  undesirable. 

Kuo  [4:21]  states  that  a  signal  can  be  completely  defined,  even  if  sampled  at  less  than 
twice  the  folding  frequency,  if  the  derivatives  of  the  signal  are  also  sampled  at  the 
sampling  instants.  The  downside  to  this  is  an  increase  in  the  number  of  sensors  needed  to 
characterize  the  signal,  or  the  difficulty  associated  with  determining  the  derivatives. 
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In  physical  systems  there  may  be  a  need  to  control  modes  that  would  be  rendered  either 
unobservable  or  uncontrollable  by  filters.  Cost  limitations  always  impact  the  number  of 
sensors  or  the  speed  of  those  sensors  or  the  hardware  necessary  to  interpret  their  output. 
Aliasing  at  times  is  unavoidable  and  may  even  be  the  desirable  alternative  for  some  of  the 
reasons  mentioned  above.  This  research  effort  proposes  to  re-examine  aliased  systems 
and  determine  some  basic  means  and  methods  for  designing  control  systems  for  such 
aliased  systems. 


1.3  Previous  Research  IN  Aliasing 

Very  little  research  has  been  found  in  the  area  of  aliasing.  This  must  be  due,  in  a  large 
part,  to  the  assumptions  mentioned  above  regarding  sampling  rates  and  low-pass  filters. 
However,  a  few  efforts  have  been  completed  in  this  direction.  They  are  briefly 
summarized  here. 

Ehrlich,  Taussig,  and  Abramovitch  [5]  introduce  a  method  to  determine  the  dynamics  of  a 
continuous-time,  linear,  time-invariant  (LTI)  system  to  frequencies  well  beyond  the 
Nyquist  frequency.  Given  a  stable  closed-loop  configuration,  using  single-frequency 
inputs,  and  measuring  the  input  and  output  spectral  amplitudes  using  a  spectrum  analyzer, 
they  identify  the  open-loop  plant  transfer  function  of  the  LTI  system,  Gp(jw),  from  the 
closed-loop  response.  They  take  the  aliased  transfer  function  and  divide  its  normalized 
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spectrum  by  the  normalized  spectrum  of  the  zero-order  hold  function  to  give  an  estimate 
of  the  discrete-time  closed-loop  transfer  function.  The  natural  frequency  of  the  system  of 
their  example  is  four  times  the  sampling  frequency.  The  effects  of  the  zero-order  hold  on 
aliasing  are  demonstrated  and  used  in  their  research  as  they  will  be  in  this  research  effort. 

Dumey  [6]  demonstrates  that  a  robot  controller’s  on-board  hardware  and  sensors  can  be 
effectively  used  to  identify  the  system’s  dynamics,  thus  allowing  for  on-line  control  law 
design.  In  his  work,  he  develops  a  model  for  aliased  frequency  response  functions 
(AFRF)  and  shows  that,  given  an  original  linear  plant,  these  AFRFs  are  also  linear.  The 
linearity  of  the  AFRF  will  be  used  in  this  research  effort. 

Kang  and  Yang  [8]  examine  the  effects  of  time  delay  in  stabilization  of  noncollacted 
control  of  beam  vibrations.  In  their  research,  they  show  the  effects  on  stability  of 
frequency  folding  when  there  is  aliasing.  They  introduce  a  controller  phase  plot  folded  at 
the  Nyquist  frequency  and  show  that  it  is  possible  to  stabilize  the  infinite  number  of 
modes  introduced  by  sampling.  This  research  revisits  the  infinite  modes  introduced  by 
sampling  and  develops  root  locus  methods  for  design  with  respect  to  these  infinite  modes. 

Yang  and  Hu  [7]  develop  further  Kang  and  Yang’s  controller  phase  plot,  and  analytically 
show  they  can  control  aliased  modes.  If  the  aliased  mode  lies  within  the  folded  phase 
plot,  the  aliased  mode  will  be  stabilized,  and  hence,  controllable.  The  idea  of  controlling 
aliased  modes  is  developed  more  thoroughly  in  this  research.  Furthermore,  a  drawback  to 
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their  research  is  that,  other  than  stability,  other  system  response  characteristics  cannot  be 
obtained  from  the  phase  plots. 

Sims  [9]  presents  the  use  of  aliased  frequency  response  functions  in  controller  design. 
His  work  shows  AFRFs  can  predict  closed-loop  stability,  and  AFRFs  are  useful  tools  in 
the  design  of  controllers  implemented  without  smoothing  or  antialias  filters,  actually 
showing  improved  performance  over  such  systems.  Sims  shows  that  systems 
performance  can  actually  be  improved  at  times  with  aliased  systems,  a  fact  which  also 
becomes  evident  in  this  research  effort. 


Dormido,  Lopez-Rodriguez,  Fdez-Marron,  and  Aranda  [10]  present  a  method  of  plotting 
a  root  locus  using  the  sample  period  as  the  parameter  of  interest.  While  not  directed 
specifically  at  aliasing,  their  results  led  to  the  use  of  the  sample  period  as  a  design 
parameter  in  the  aliased  systems  of  this  research  effort. 

Lyman  [10a]  further  develops  the  method  of  root  locus  plotting  using  the  sample  period 
as  the  parameter  of  interest.  He  allows  aliasing  and  determines  its  affect  on  the  z-plane 
root  locus  of  the  systems.  He  further  examines  the  effects  of  time  delay  on  these  root  loci 
plots.  His  root  locus  method  is  useful  in  designing  controllers  of  sampled-data  systems. 
His  work  is  directly  applicable  to  this  research  effort,  as  the  effects  of  aliasing  and  sample 
period  selection  on  the  root  loci  of  the  mixed  continuous  and  discrete  systems  of  this 
research  is  used  extensively. 
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1.4  Related  Research 


Although  not  directly  related  to  efforts  on  designing  controllers  for  systems  in  which 
aliasing  is  allowed,  much  current  research  is  being  conducted  in  areas  of  peripheral 
interest  to  this  research  effort.  Some  of  these  efforts  are  summarized  below. 

Wang,  Wang,  and  Lee  [33]  explore  the  effects  of  sampling  time  on  the  higher-order 
digitizations  used  to  convert  a  continuous  system  into  the  discrete  domain.  The  two 
higher-order  digitizations  include  the  Madwed  and  Boxer-Thaler  methods  [34]. 
Comparisons  are  made  between  these  digitizations  and  the  Tustin  method.  They 
introduce  the  z-plane  root  locus  based  on  sample  period,  which  they  call  the  “polynomial 
root  locus,”  for  all  three  digitization  methods.  Stability  is  determined  by  the  unit  circle  in 
the  z-plane.  Additionally  they  developed  a  new  method  to  transform  the  system  to  the  w- 
plane,  allowing  use  of  the  Routh-Hurwitz  stability  criterion  of  a  continuous  system  in  the 
determination  of  the  maximum  sample  time. 

Diduch  and  Doraiswami  [35]  perform  a  study  on  the  influence  of  the  sample  period  on  a 
system  in  which  the  controller  is  a  simple  servomechanism.  They  look  at  how  the  sample 
period  affects  both  the  sensitivity  of  the  design  and  the  performance  of  the  system.  Their 
simulations  show  that  the  performance  of  the  closed-loop  system  degrades  monotonically 
with  the  sample  period,  leading  to  the  obvious  conclusion  that  a  sample  period  of  zero  is 
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optimal.  They  also  demonstrate  a  decrease  in  robustness  and  an  increase  in  sensitivity 
with  an  increase  in  the  sample  period. 

Hirata  and  Powell  [36]  examine  the  effects  of  sample  rate  on  broadband  disturbance 
rejection  in  digital  control  systems.  They  show  that  the  selection  of  the  correct  sample 
rate  has  a  great  effect  on  the  system’s  ability  to  reject  disturbances.  They  show  that  with 
a  sample  rate  as  low  as  eight  times  the  closed-loop  bandwidth  there  is  a  50%  degradation 
in  the  discrete  system’s  disturbance  rejection  in  comparison  to  a  continuous  controller. 
Hirata  also  discusses  peaks  found  in  the  cost  survey  of  a  linear  quadratic  Gaussian  (LQG) 
approach  used  to  evaluate  performance  [37:43-47].  These  peaks  occur  at  integer  fractions 
of  twice  the  natural  frequency  of  the  systems  he  researched.  He  states  that  at  these 
frequency  ratios  the  systems  he  studies,  in  particular  the  mass-spring-mass  system,  are 
uncontrollable,  since  their  cost  functions  peak  quite  sharply.  In  this  dissertation, 
however,  at  these  frequency  ratios  the  systems  prove  to  be  controllable  given  the  fact  they 
include  damping.  As  damping  approaches  zero,  the  cost  of  controlling  at  these  frequency 
ratios  do  indeed  rise  drastically. 

Kranc  introduced  the  concept  of  analysis  of  multirate  feedback  systems  in  1957  [27].  He 
developed  a  general  technique  that  extended  the  z-transform  methods  to  sampled-data 
systems  in  which  synchronized  switches”  did  not  operate  with  the  same  sampling 
frequency  as  the  overall  system.  This  method  could  be  used  to  examine  system  output 
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between  sampling  instants  by  introducing  a  sampler  operating  at  a  rate  greater  than  that  of 
the  system. 

Tomero,  Gu,  and  Tomizuka  [23]  analyze  multi-rate  controllers  that  update  the  controller 
output  faster  than  the  measurement  sampling  frequency  by  an  integer  factor.  The  systems 
they  control  are  continuous  and  the  controller  is  discrete  with  a  zero-order  hold.  As  a 
result  of  the  multi-rate  controller  the  control  input  becomes  oscillatory.  They  theorize, 
however,  that  the  multi-rate  controller  may  enable  stability  of  a  system  that  is  unstable 
with  a  single-rate  controller. 

Tangirala,  Li,  Patwardhan,  Shah,  and  Chen  [24]  also  examine  multi-rate  controllers  in 
which  the  controller  receives  input  at  a  greater  frequency  than  the  output.  They,  however, 
introduce  the  concept  of  lifting  techniques  [25]  in  which  a  fast-rate  signal  is  mapped  to  a 
slow-rate  signal  with  an  increased  dimensionality,  and  inverse  lifting  techniques  in  which 
the  mapping  is  from  a  lifted  signal  to  a  fast-rate  one.  They  demonstrate  that  this  discrete 
lifting  will  introduce  intersample  ripples  in  the  outputs  of  the  closed-loop  multi-rate 
systems  if  the  steady-state  gains  of  the  lifted  models  are  non-uniform.  If  the  gains  are 
identical,  however,  there  are  no  oscillations  of  the  outputs.  They  additionally  show  that  a 
fast-rate  integrator  will  eliminate  these  oscillations  no  matter  what  the  gains  are. 

Cavicchi  [26]  introduces  the  concept  of  a  phase-root  locus  and  the  resulting  stability 
measures  than  can  be  examined  from  this  plot.  Cavicchi  develops  the  phase-root  locus  as 
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the  dual  of  the  conventional  root  locus,  it  being  a  plot  of  the  closed-loop  pole  movement 
as  phase  is  added  to  the  open-loop  transfer  function.  The  phase-root  locus,  plotted 
together  with  the  traditional  root  locus,  allow  the  s-plane  to  be  used  to  determine 
robustness,  perform  transient  design,  and  aid  in  stability  analysis. 

Diduch  and  Doraiswami  [38]  show  that  when  unmodeled  process  dynamics  are  included 
in  the  analysis,  an  increased  sample  time  may  actually  improve  the  robustness  of  an 
optimally  designed  sampled-data  system.  This  is  because,  with  the  inclusion  of  the 
unmodeled  dynamics,  the  mapping  of  the  uncertainty  from  the  continuous  to  the  discrete 
domain  shrinks  or  contracts  as  the  sample  period  is  increased.  This  would  tend  to 
indicate  that  stability  may  be  increased  with  an  increase  in  sample  period  as  shown  in  this 
dissertation. 


Braslavsky,  Middleton,  and  Freudenberg  [39]  introduce  a  frequency  domain  lifting 
technique  to  analyze  what  they  call  hybrid  systems — systems  in  which  they  wish  to 
examine  intersample  behavior.  They  demonstrate  the  application  of  their  technique  to  the 
stability  robustness  of  an  LTI  and  non-LTI  system. 

Much  has  been  researched  with  respect  to  the  effects  of  time  delay  on  the  stability  of 
systems.  Perhaps  the  “grand-daddy”  of  this  research  effort  is  that  presented  by  Ts)q)kin 
[28].  He  derived  stability  conditions  for  systems  in  which  the  feedback  is  delayed.  His 
development  allowed  simplification  of  such  systems. 
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Other  research  into  the  effects  of  time  delay  include  the  work  done  by  Lyman  [12]  on  root 
loci  with  time  delay  as  the  parameter  of  variation.  He  verifies,  based  on  the  root  loci 
plots,  Tsypkin’s  three  types  of  behavior:  stability  for  all  delay,  stability  to  a  critical  value 
of  delay,  and  a  finite  alternation  between  stability  and  instability  as  delay  increases. 

Other  research  in  time  delay  effects  includes  stability  criteria  introduced  by  Hou  and  Ali 
[29],  compensation  methods  for  computational  delay  by  Rattan  [30]  and  Ovaska  and 
Vainio  [31],  and  the  effects  of  zeros  in  digital  control  systems  in  which  computational 
time  delay  is  present  by  Kara,  Kondo,  and  Katori  [32]. 

1.5  Dissertation  Outline 

This  research  effort  begins  by  developing  a  model  for  a  system  in  which  the  effects  of 
aliasing  on  continuous  output  can  be  examined.  The  model,  then,  will  consist  of  discrete 
and  continuous  portion.  Due  to  this  mixed  nature,  examination  of  the  system  poles  and 
zeros  must  occur  in  a  common  plane.  The  s-plane  is  valid  for  examination  of  the 
continuous  roots  of  the  system  and  the  z-plane  is  valid  for  examination  of  the  discrete 
roots  of  the  system,  but  neither  is  completely  accurate  for  examination  of  discrete  and 
continuous  roots  together. 
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A  new  plane  for  examination  of  these  poles  and  zeros  is  first  introduced:  the  s*-plane. 
This  plane  allows  for  examination  of  the  continuous  poles  and  zeros  with  their  locations 
determined  as  they  are  in  the  traditional  s-plane.  Discrete  poles  and  zeros  and 
transformed  to  the  primary  band  of  the  s*-plane  from  the  z-plane  and  repeated  infinitely 
into  every  other  band  formed  according  to  sample  rate. 

A  model  of  a  mixed  discrete  and  continuous  system  is  developed  next;  this  model  is 
contrary  to  the  strictly  discrete  or  continuous  models  used  for  common  study  of  sampled- 
data  systems.  In  the  model  developed  here,  the  plant  will  be  a  continuous  plant  and  the 
controller  will  be  a  discrete  controller.  Input  will  be  restricted  to  discrete  input  only.  The 
systems  studied  in  this  research  effort,  then,  are  sampled-data  systems  in  which  the  input 
and  feedback  are  discrete,  but  the  output  is  continuous. 

For  this  model  to  be  valid  it  must  equal  a  discrete  system  in  the  limit  as  sampling  period 
increases,  and  it  must  equal  a  continuous  system  in  the  limit  as  sampling  period  goes  to 
zero.  Limit  analysis  shows  this  to  indeed  be  the  case.  In  strictly  discrete  or  continuous 
systems,  the  open-loop  poles  of  the  plant  are  cancelled  by  a  set  of  the  system’s  closed- 
loop  zeros,  allowing  for  placement  of  the  closed-loop  poles  to  enable  control.  This 
research  shows  this  to  be  the  case  for  the  mixed  discrete  and  continuous  system  under 
consideration,  allowing  control  of  the  continuous  plant  with  the  discrete  controller. 
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With  the  model  developed  and  the  s*-plane  introduced,  the  effects  of  aliasing  on  classical 
measures  of  stability  such  as  the  root  loci  and  frequency  response  function  plots  are 
investigated,  with  a  determination  of  acceptable  aliasing  regions  introduced.  Design 
methods  using  root  locus  tools  are  extended  from  the  traditional  methods  of  non-aliased 
system  design.  The  impulse  and  step  responses  of  these  systems  are  examined,  as  well  as 
how  well  these  systems  follow  a  sinusoidal  input.  The  ability  of  these  aliased  control 
systems  to  follow  a  sinusoidal  input  is  poor  unless  the  frequency  of  the  sine  wave  is  less 
than  the  sampling  frequency. 

This  research  effort  includes  only  second-  or  fourth-order,  linear,  time-invariant  systems. 
First,  a  second-order  plant  with  no  zeros  is  examined  and  a  relationship  between  sampling 
frequency  and  root  loci  locations  is  developed.  Sampling  frequencies  are  found  which 
actually  introduce  stability  with  increased  aliasing.  A  deadbeat  response  of  the  closed- 
loop  system  is  developed  using  simple  proportional-integral  control,  with  the  response 
settling  to  the  final  values  within  one  sample  period.  This  deadbeat  controller  requires 
very  accurate  knowledge  of  the  controlled  system’s  natural  frequency  and  damping  ratio. 
Robustness  is  examined  and  found  to  be  adequate.  Disturbance  rejection  is  also 
examined.  It  was  found  the  system  can  reject  step  disturbances  within  one  sample  period 
of  the  feedback  of  the  disturbance  to  the  controller.  This  system  does  not  have  the  ability 
to  reject  broadband  disturbances. 
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A  simple  fourth-order  plant  with  no  zeros  is  investigated  next.  The  high-frequency  mode 
of  this  plant  is  aliased,  and  the  effects  of  this  high-frequency  are  minimized  by  allowing 
aliasing.  Again,  the  increased  aliasing  actually  improves  response.  With  the  effects  of 
the  aliased  high-frequency  mode  minimized,  the  design  of  control  for  the  low-frequency 
mode  follows  traditional  design  methods.  Again,  very  accurate  information  of  the  aliased 
mode’s  natural  frequency  and  damping  ratio  are  required.  Robustness  and  disturbance 
rejection  results  were  similar  to  those  of  the  second-order  system. 

The  contribution  of  the  research  effort  is  easy  to  determine.  As  a  result  of  the  research, 
aliasing  becomes  acceptable  in  certain  design  regions;  indeed,  aliasing  caused  by  a  slow 
sample  rate  even  becomes  a  design  parameter,  introducing  stability  and  allowing  the 
control  of  the  high-frequency  aliased  modes.  In  areas  where  aliasing  is  unavoidable,  the 
development  of  these  design  tools  will  become  very  valuable.  The  limitation  to  this 
method  is  its  inability  to  follow  high-frequency  inputs  or  reject  high  frequency 
disturbances.  If,  however,  the  input  and  disturbances  are  at  frequencies  below  the 
sampler,  this  method  provides  an  excellent  means  to  control  aliased  modes. 
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Chapter  Two.  Theory 


2.0  Introduction 

In  this  chapter  we  develop  the  basic  theory  and  model  for  the  study  of  systems  in  which 
aliasing  is  allowed.  This  development  depends  on  the  reader’s  understanding  of  Laplace 
transforms,  z-transforms,  complex  logarithms,  root  locus  plots,  the  sampling  operation 
and  the  definition  of  an  ideal  sampler,  and  the  effects  of  a  zero-order  hold  on  a  plant.  A 
brief  review  of  these  principles  can  be  found  in  Appendix  A. 

This  chapter  begins  with  the  introduction  of  the  s*-plane,  a  plane  in  which  both  the 
continuous  and  discrete  poles  and  zeros  can  be  examined.  With  the  introduction  of  the 
s*-plane  and  a  basic  understanding  of  sampling,  aliasing,  and  the  zero-order  hold,  the 
model  used  for  examination  of  aliased  systems  can  now  be  developed.  Here  the  model 
consists  of  both  discrete  and  continuous  parts,  and  is  thus  named  a  discrete-to-continuous 
system.  From  this  model,  a  discrete-to-continuous  transfer  function  is  derived.  The 
model  is  examined  in  the  limits  as  the  sample  period  goes  to  zero  and  infinity  and  found 
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to  reduce  to  the  continuous  and  discrete  system,  respectively.  For  control  of  the  closed- 
loop  discrete-to-continuous  system  to  occur,  the  open-loop  poles  must  be  cancelled  by  a 
set  of  the  closed-loop  zeros.  This  is  found  to  be  the  case  for  the  mixed  system.  And 
finally,  prior  to  plotting  the  s*-plane  root  locus,  an  examination  is  made  of  the  origins  of 
the  closed-loop  poles  and  zeros  and  the  effects  of  the  sample  period  on  these  roots. 


2.1  The  S*Plane 

Traditionally,  root  locus  analysis  has  been  performed  in  the  s-  or  Laplace  domain  for 
continuous  systems  and  in  the  z-domain  for  discrete  systems.  In  this  research  effort, 
however,  examination  will  be  of  a  system  containing  both  discrete  and  continuous 
portions  in  the  transfer  function.  Transforming  the  continuous  poles  and  zeros  to  the 
discrete  domain,  although  mathematically  possible,  means  little  in  terms  of  design  with 
slow  sample  rates.  Conversely,  transforming  the  discrete  poles  and  zeros  to  the  s-plane 
means  little  since  these  discrete  poles  and  zeros  are  repeated  infinitely  due  to  the  sampler. 
To  have  a  common  plane  in  which  to  examine  both  the  continuous  and  discrete  poles  and 
zeros,  we  will  introduce  the  s*-plane.  It’s  development  follows. 

The  mapping  from  the  s-domain  to  the  z-domain  is  common,  and  tables  included  in  most 
digital  control  theory  texts  [1,  31, 32,  33]  have  s-domain  as  well  as  z-domain  forms  of 
basic  time  functions,  allowing  the  transform  from  the  s-domain  to  the  z-domain.  A  word 
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of  caution  is  in  order  however.  When  transforming  from  the  s-domain  to  the  z-domain 
one  must  ensure  the  system  to  be  mapped  is  in  pole-residue  form — the  tables  are 
constructed  based  on  this  fact.  A  z-transform  (Appendix  A)  of  the  poles  and  zeros  using 

z  =  [2.1] 

would  map  the  zeros  to  different  locations  than  finding  the  pole-residue  form,  performing 
the  z-transform,  and  then  recombining  to  calculate  the  zero  locations,  bi  other  words, 
poles  map  using  Equation  2.1,  but  zeros  do  not!  The  reason  for  this  will  be  explored 
later. 

The  mapping  from  the  s-plane  to  the  z-plane  is  a  many-to-one  mapping.  The  reason  for 
this  follows:  The  s-plane  is  first  divided  into  infinitely  many  bands  of  “height”  ,  as 

shown  in  Figure  2.1.  The  primary  band  extends  from  ©=  to  ©=-i-— . 


Using  the  relation  of  Equation  2.1  the  points  |^,  |],  Q,  and  U  are  mapped  to  the  z-plane. 

Corresponding  points  are  shown  in  Figure  2.2. 
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Since 


^(s+jn(S),)T,  _^sT,^jna,T,  _^sT, 


[2.2] 


.  .  .  ,  ,  j  j.  ^  3(0,  -  CO.  3co. 

each  of  the  complementary  bands  extending  from  — ^  to - from  H — ^  to  H - ^ 

2  2  2  2 


and  so  forth,  map  to  the  same  location  in  the  z-plane.  Thus,  this  mapping  is  a  many-to- 
one  mapping.  This  fact  also  explains  the  non-uniqueness  of  the  inverse  z-transform.  hi 
performing  this  inverse  z-transform,  one  cannot  define  to  which  band  the  mapping  should 
occur — ^it  is  a  one-to-many  mapping.  The  mapping  from  the  z-plane  to  the  s*-plane, 
however,  will  account  for  this  fact. 


With  the  above  mentioned  caution  of  transformations  using  the  pole-residue  form  it  is 
relatively  straightforward  to  translate  from  the  s-plane  into  the  z-plane,  and  no  difficulties 
are  encountered  with  the  sample  period,  .  In  transforming  from  the  z-plane  to  the 

s*-plane,  however,  because  of  the  one-to-many  mapping,  there  exist  multiple  bands  to 
which  the  z-plane  poles  can  be  mapped.  Additionally  when  considering  pole-zero  or 
pole-residue  form,  one  must  know  what  form  is  invoked  by  performing  the  inverse 
z-transform. 


An  inverse  z-transform  of  the  pole-residue  form  of  F (z)  gives  the  equation  for  f(kT^) , 
which  may  or  may  not  equal  /  *  (t)  because  of  it’s  non-uniqueness  (see  Equation  A.  14). 
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Taking  the  inverse  z-transform  of  F(z)  to  the  s*-plane  provides  the  pole-residue  form  of 

F*{s),  if  the  poles  and  residues  of  are  mapped  to  the  primary  band  and 

repeated  in  every  band  infinitely.  The  infinite  sum  of  these  poles  and  residues  would  then 
provide  the  pole-zero  form  of  F  *  (s) . 

The  easier  way  of  obtaining  this  pole-zero  form,  however,  is  to  simply  map  the  zeros  and 
poles  from  the  z-plane  to  the  s*-plane  using  the  inverse  relation  of  Equation  2.1: 


[2.3] 


Once  the  poles  and  residues  are  combined  into  the  poles  and  zeros  in  the  z-plane,  the 
infinite  sum  is  already  accomplished  because  of  the  many-to-one  mapping  that  occurs  in 
the  transform  from  the  s-plane  to  the  z-plane.  These  z-plane  poles  and  zeros  are  then 
mapped  to  the  primary  band  of  the  s*-plane  using  Equation  2.3,  and  repeated  infinitely  to 
produce  the  pole-zero  form  of  F  *  (^) . 


This  infinite  repetition  to  produce  the  pole-zero  form  of  F*(s)  is  the  result  calculating  the 
logarithm  of  a  complex  number  (Appendix  A).  Because  of  the  nature  of  calculating  the 
argument  of  z ,  the  mapping  from  the  z-plane  to  the  s*-plane  is  a  one-to-many  mapping. 
Each  of  the  poles  and  zeros  mapped  from  the  z-plane  to  the  principal  or  primary  band  in 
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the  s*-plane  is  repeated  infinitely  to  every  other  band  at  intervals  of  jnco^ ,  where 
n  =  0,±l,±2,...  . 

The  difference  between  the  mapping  of  the  pole-residue  form  and  the  pole-zero  form  is  a 
direct  result  of  aliasing.  In  the  limit  as  7^  goes  to  zero,  the  width  of  the  primary  band  in 

the  imaginary  axis  direction  becomes  infinite.  Thus,  the  effect  of  the  poles  and  residues 
does  not  overlap  across  bands,  and  the  infinite  sum,  which  produces  the  pole-zero  form, 
equates  to  the  infinite  product  of  that  pole-zero  form.  As  7^  goes  to  zero,  the  pole-zero 

and  the  pole-residue  forms  become  the  same.  As  increases,  however,  the  “tails”  of  the 

poles  and  residues  cross  the  band  boundaries,  adding  to  the  effects  of  the  poles  and 
residues  of  other  bands,  and  the  infinite  sum  has  to  include  these  influences.  Thus,  as 
aliasing  increases  the  location  of  the  zeros  will  move,  and  the  poles  and  zeros  found  by 
mapping  the  poles  and  residues  using  Equation  2.3  will  not  produce  the  same  zero 
locations  as  mapping  the  poles  and  zeros  unless  the  infinite  summation  is  performed. 

One  may  notice  the  similarity  in  concept  of  mapping  from  the  z-plane  to  the  s*-plane  as 
the  bilinear  transformation  made  to  map  the  unit  circle  from  the  z-plane  to  the  imaginary 
axis  of  the  w-plane  [34:225-230].  This  w-plane  mapping  allows  the  use  of  the  continuous 
forms  of  the  Routh-Hurwithz  criterion  and  Bode  techniques  in  the  examination  of  a 
discrete  system.  The  difference  here,  however,  is  that  the  w-plane  mapping  is  a  one-to- 
one  mapping  and,  as  such,  does  not  allow  for  the  examination  of  aliasing  that  the  s*-plane 
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mapping,  a  one-to-many  mapping,  does.  Perhaps  a  mapping  closer  to  that  of  the  s*-plane 
would  be  that  one  found  from  using  the  modified  z-transform  [33] 


2.2  Model  Development 

The  model  for  this  research  effort  is  shown  below  in  Figure  2.3.  Parameters  resulting 
from  discrete  components  are  designated  by  the  asterisk  (*);  therefore,  the  sampled 
version  of  the  output  signal  in  the  frequency  domain,  C(.s) ,  is  designated  as  C*  (.y) . 


If  (ZOH)*(Gjj(s))  =  H(s),  then 


Figure  2.3  Sampled-Data  System  Model 
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The  input  to  the  sampled-data  system  will  be  restricted  to  change  only  at  the  sample 
instants,  hence  the  input  is  shown  as  R*  ( 5') .  The  zero-order  hold,  ZOH(s) ,  and  the 
plant,  Gp{s) ,  will  be  combined  into  a  function  called  H(s) .  The  control  law, 
implemented  by  the  computer,  is  discrete  and  is  designated  as  *(s) . 


The  following  two  relations  can  be  derived  from  Figure  2.3: 

C(s)  =  G,*is)His)E*is) 

E*(s)  =  R*(s)-C*(s) 

Substituting  Equation  2.5  into  Equation  2.4  gives 

C(.s)  =  G,*  (s)H(s)[R  *(s)-C*  ( j)] 

=  G,  *  is)H(s)R  *  (s)  -  G,  *  (s)His)C  *  (s) 


[2.4] 

[2.5] 


[2.6] 


Applying  the  sampling  relationship  (Equation  A.  18,  Appendix  A)  to  Equation  2.6  to  find 
an  expression  for  the  sampled  version  of  C(s)  gives 


C*(s)  =  [G,  *  (s)H(s)R  *  (5)  -  G,  *  (s)H(s)C  *  (5)]  * 
=  G,  *  (s)H  *  (s)R  *  (s)  -  G,  *  (s)H  *  (s)C  *  (s) 


[2.7] 
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Collecting  like  terms  and  solving  for  C  *  (^)  gives 


C\s) 


H\s)G:(s) 

1  +  H\s)G:(s) 


R\s)  . 


[2.8] 


Substituting  Equation  2.8  into  Equation  2.6: 


C(.y)= 


H  *  (s)G^  *  (5) 
l  +  H*(s)G,*(s) 


=^G,*is)H(s) 


1— 


H*(s)G,*(s) 
l  +  H*(s)G^*(s) 


\R*is) 


=  G,*is)H(s) 


1+H* (s)G,  (s)G^  * (s) 

l  +  H*(s)G^*(s) 


'R*(s) 


C(s)  = 


G,*(s)H(s) 
l  +  G^*(s)H*is) 


R*(s)  . 


[2.9] 


Equation  2.9  suggests  a  discrete-to-continuous  transfer  function  in  the  form  of  continuous 
output  over  discrete  input: 


Cjs) 

R\s) 


GlisWis) 

l  +  G:(s)H*is) 


=  F*{s)H(s). 


[2.10] 
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Systems  having  both  discrete  and  continuous  signals  are  traditionally  known  as  sampled- 
data  systems.  Such  systems  have  been  studied  extensively.  However,  examination  has 
been  in  either  the  digital  (discrete)  or  continuous  domain.  Here,  the  discrete  and 
continuous  domains  are  examined  together.  The  discrete-to-continuous  transfer  function 
of  Equation  2.10  includes  both  the  continuous  (typically  aperiodic)  and  the  discrete 
(periodic)  frequency  components,  and  as  such  is  not  a  transfer  function  in  the  classical 
sense  of  the  term.  The  frequency  response  function  (FRF)  corresponding  to  this  discrete- 
to-continuous  transfer  function  will  be  examined  below.  The  FRF  of  the  sampled-data 
transfer  function  will  be  termed  a  discrete-to-continuous  frequency  response  function 
(DCFRF). 

A  DCFRF  for  the  system  discussed  earlier  (cOh  =100  rad/sec,  ^  =  0.01)  is  plotted  below  in 
Figure  2.4  for  a  sampling  frequency  of  0)^;  =  90  rad/sec.  The  DCFRF  includes  only 
proportional  control,  with  *  (5)  =  0. 1 .  Examination  of  Figure  2.4  clearly  shows  the 
effects  of  aliasing,  sampling,  and  the  zero-order  hold  (ZOH). 

Although  typical  intuition  cannot  be  applied  to  the  DCFRF,  several  results  can  be  seen  in 
Figure  2.4.  The  aliasing  effect  can  be  seen  as  the  resonant  frequency  is  folded  about  the 
folding  frequency  (cOf  =  45  rad/sec).  All  frequency  components  above  45  rad/sec  are 
folded  into  the  primary  band  (-45  <  co<  45).  Sampling  effects  are  seen  in  the  repetition 
of  the  primary  band  characteristics  at  multiples  of  the  sampling  frequency.  And  finally. 
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Frequency  Ratio  Ws/Wn  =  0.9;  Kp  =  0.1 


Figure  2.4  Discrete-to-Continuous  Frequency  Response  Function 

at  00^  =  90  rad/sec 


the  effects  of  the  ZOH  can  be  seen  as  the  magnitude  decreases  while  frequency  increases, 
and  with  the  magnitude  plot  going  to  zero  at  integer  multiples  of  the  sampling  frequency. 


The  discrete-to-continuous  impulse  response  for  this  system  is  shown  in  Figure  2.5. 
Sampling  points  in  time  are  shown  by  the  asterisks,  and  the  continuous  output  of  the 
sampled-data  system  is  the  solid  line.  Here,  the  aliased  mode  of  the  sampled  system  is 
easily  observable,  as  is  the  natural  response  of  the  system.  It  appears  that  between  sample 
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Impulse  Response 


0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 


Time  (sec) 

Figure  2.5  Impulse  Response  of  Aliased  Sampled-Data  System  at 

oofe  =  90  rad/sec 


periods  the  continuous  system  responds  at  its  open-loop  undamped  natural  frequency. 
This  will  indeed  be  shown  to  be  the  case  later  on. 

In  this  dissertation,  we  will  primarily  be  concerned  with  the  unique  F*(s),  and  as  such 
will  be  using  the  inverse  z-transform  of  the  pole-zero  form  of  F(z) .  Traditionally,  poles 
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and  zeros  are  mapped  to  a  primary  band  existing  between  s  =  — ^  and  s  =  ^.  For  our 

2  2 

analysis,  however,  we  are  also  interested  in  the  poles  and  zeros  that  are  repeated  at  ± 

intervals  infinitely  throughout  the  s-plane  and  their  eflFect  on  the  continuous  response. 
Therefore,  any  plot  of  the  root  locus  in  the  s*-plane  includes  these  repeated  poles  and 
zeros. 

Recall  that  the  transform  of  F(z)  from  the  z-plane  to  the  s-plane  using  the  inverse 
z-transform  method  gives  one  of  an  infinite  number  of  functions  F(j)  that  will  describe 
the  function  f(t)  only  at  the  sample  instants.  The  unique  transform  of  F(z)  to  F*(s), 
however,  gives  the  unique  function,  which  is  also  the  F(s)  function  that  will  excite  the 
analog  plant,  G^is) .  The  plant  is  only  excited  at  the  sample  points  with  impulses  input 
into  the  zero-order  hold.  The  F  *  (j) ,  then  is  the  only  appropriate  controller  function. 


As  explained  before,  the  unique  form  of  F  *  (5)  which  represents  a  discrete  time  function 
is  periodic  in  frequency.  The  plant  His) ,  pre-multiplied  by  the  zero-order  hold,  is 
aperiodic  in  frequency  and  represents  a  continuous  time  function.  When  His)  is 
multiplied  by  F  *  (5) ,  then,  the  resulting  output  of  our  system  will  become  aperiodic  in 
frequency  and  continuous  in  time. 
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2.3  Discrete-To-Continuous  System  In  Sample  Period  Limits 


The  equation  for  the  discrete-to-continuous  transfer  function  is  defined  in  Equation  2.10. 
This  equation  should  reduce  to  the  continuous  transfer  function  in  the  limit  as  the  sample 
period,  ,  goes  to  zero.  The  following  shows  this  is  indeed  the  case. 


Recall  H(s)  is  defined  as  the  plant,  Gp(s) ,  multiplied  by  the  zero-order  hold: 


His)  = 


1  —  e  ^ 


GM)  . 


[2.11] 


The  limit,  then,  as  the  sample  period  goes  to  zero  is 


lim 

r,-»o 


H(s) 

Ts 


=  lim 
r,-»o 


sTs 


Gp(s)  = 


[2.12] 


Applying  L’Hopital’s  rule  gives 


lim 

Ts^O 


H(s) 

Ts 


=  lim 


se 


-sT, 


S 


G,(s)  =  -G,W  =  G,(j) 
s 


[2.13] 


33 


Now  examine  the  H  *  (s)  term  of  the  denominator  of  — ,  Por  our  system  the 

T,R*(s) 


ZOH  is  included  with  the  plant 


H(s)  =  ZOH(s)G^(s)  . 


[2.14] 


Since 


ZOH(s)=— — 
s 


[2.15] 


the  resulting  product  of  the  zero-order  hold  and  plant  is 


H(s)  = 


Np(s) 
Dp(s)  ’ 


[2.16] 


and,  using  the  equation  for  a  sampled  signal,  Equation  A.  16, 


T  ^ 

n-- 


Nj,{s-¥  jnaj 
s  +  jn(d^  Dpis  +  jnoij' 


[2.17] 
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Expanding  the  exponential  term  in  the  numerator  gives 


g-is+jnai,)T, 


The  relationship  between  the  sample  time  and  the  sample  frequency  is 


so 


-jn—T, 

g-Jn<oJs  _  g  T,  _  g-J'imt  _  j 


for  n  =  integer  . 


Substituting  this  into  Equation  2.17,  the  equation  for  H*(s),  gives 


y  Np(s  +  jnoi,) 

Ts  (s  +  jnd),  )Dp(s  +  jn(d, ) 


[2.18] 


[2.19] 


[2.20] 


[2.21] 
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Looking  at  the  limit  of  the  summation: 


^r> 

...  -L 

(s  1 

i  0  J  .  Np(s) 

is 

limT  N{s  +  jn(SiA  _  1  0  j 

Np  s  + 
J _ V  _ 

d(s  ^'^’'1 

1  0  J 

j27C^ 

0  ^ 

y  1 

(  y27i\ 

L 1  0 ; 

% 

1  0  J  J 

Np(s  +  jn(i)^) 


OO 

lim  V  - 

(s  +  jn(0^  )D  (j  +  jndi^ ) 


NAs)  1 


[2.22] 


The  zeros  of  the  infinite  summation  term  occur  because  the  order  of  the  denominator  of 
the  plant  will  always  be  greater  than  or  equal  to  the  order  of  the  numerator  (physical 
realizability).  Thus  with  the  zero-order  hold  term,  the  order  of  the  denominator  terms 

will  always  be  at  least  one  greater  than  that  of  the  numerator,  causing  the  limit  to  be 

OO 

which  is  zero. 


The  limit  of  the  non-summation  term  is 


lim 


1-10 
0  “O  ’ 


[2.23] 
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so  L’Hopital’s  rule  must  be  applied: 


\im— - =lim^^ - =  5  .  [2.24] 

T,^o  1 


Therefore, 


[2.25] 


which  is  the  expected  result. 


Finally,  examine  the  controller,  *  (s) .  For  the  sake  of  comparison  to  the  continuous 
case,  assume  *  (s)  is  simply  the  ideally  sampled  version  of  G^is) ,  the  inverse 

mapping  of  the  controller  G^(z) . 


In  this  case 


G,*(s)  = 


1 


[2.26] 
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Following  the  same  reasoning  as  before,  in  the  limit  as  the  sample  period  goes  to  zero  this 
approaches  the  continuous  controller,  G^is) : 


limG^  *{s)  =  G^(s) 
r,-»o 


[2.27] 


Looking,  then,  at  the  limit  as  the  sample  period  goes  to  zero,  one  can  see  that  the  discrete- 
to-continuous  transfer  function  equals  the  continuous  transfer  function: 


T,^0T^R*(s) 


G,*(s) 


H(s) 


=  lim 


Gc(s)G^(s) 
l  +  G,{s)Gp(s) 


[2.28] 


Now,  examine  the  limit  of  the  discrete-to-continuous  transfer  function  as  the  sample 
period  approaches  infinity.  Begin  by  examining  the  numerator  term,  H  {s) 


lim 

r,->oo 


tm 

Ts 


lim 

r,->~ 


l-e 


sZ 


1-0 


Gpis)  =  0 


[2.29] 


This  will  be  the  limit  of  the  entire  discrete-to-continuous  transfer  function,  if  the 


denominator  terms  do  not  also  equal  zero.  The  limit  of  H*  (j)  is 
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[2.30] 


lim  H*(s)=  lim 

Tj  --400  Tj  ->oo 


1-e 


1 

«=-oo 


Npjs  +  jnG),) 

(s  +  jn(i>^)Dp{s  +  jnG),) 


1-0  y  iVp(5  +  ;nm,) 

~  Ts „~o (-5  +  >0), )T^p{s  +  jndi, ) 


Because  the  limit  of  the  term  is  zero,  the  limit  of  the  denominator  term  of  the  continuous- 
to-discrete  transfer  function  equals  1,  and 


lim 


C(s) 


[2.31] 


This  result  is  expected.  As  the  sample  period  goes  to  infinity,  the  loop  is  never  closed, 
never  allowing  the  controller  to  input  to  the  system,  and  the  system  response  is  zero. 


2.4  Cancellation  of  The  Open-Loop  Poles  by  The  Zeros  from  The 
Closed-Loop  System 


With  a  strictly  continuous  or  discrete  system,  the  process  which  allows  its  control  once 
the  loop  is  closed  is  that  fact  that  the  open  loop  poles  will  be  canceled  by  the  closed  loop 
zeros,  allowing  for  the  placement  of  the  closed  loop  poles  in  desired  locations  for  desired 
system  response.  This  fact  will  be  shown  below.  This  is  also  the  case  with  the  discrete- 
to-continuous  system  of  this  research  effort.  The  mathematical  proof  of  this  fact  follows. 
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Suppose,  in  the  continuous  domain,  the  plant  is  defined  by 


GJs)  = 


Dpis) 


and  a  controller  defined  by 


NM 

DM  ' 


If  the  loop  is  closed  with  negative  feedback  the  resulting  transfer  function  is 


c(s)_  gmgm 
R(s)  l  +  GMGp(s)’ 


where  C(s)  is  the  output  and  R(s)  is  the  input.  Using  Equations  2.32  and  2.33  in 
Equation  2.34  gives 


C(s) 

R(s) 


NM  Np(s) 
Deis)  DM 

DM  Dp(s) 


[2.32] 


[2.33] 


[2.34] 


[2.35] 
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which  after  some  algebraic  manipulation,  and  before  any  cancellation,  becomes 


C(s)  N,(s)N^(s)D,(s)D^(s) 

R(s)  D,(s)Dp(s)[D,is)Dpis)  +  N^(s)Np(s)]  ‘ 


[2.36] 


Zeros  are  those  factors  for  which  the  numerator  equals  zero.  Likewise,  poles  are  those 
factors  for  which  the  denominator  equals  zero.  Notice  that  by  closing  the  loop  we  have 
caused  a  pole-zero  cancellation — ^the  zeros  found  from  the  term  Dp  (s)  in  the  numerator 

cancel  the  poles  of  the  same  form  in  the  denominator.  When  one  thinks  this  through,  this 
indeed  must  occur  if  we  wish  to  control  the  system:  the  open  loop  poles  must  be  rendered 
inactive  by  this  cancellation  to  allow  control,  or  movement  of  the  poles,  through  design  of 
the  controller.  The  resulting  system  is 


C(5)  _  N,(s)Np(s) 

R(s)  D,(s)Dp(s)  +  N,(s)Np(s) 


[2.37] 


In  the  discrete  domain  the  analysis  is  identical,  with  each  ‘V’  in  Equation  2.32  through 
Equation  2.37  replaced  with  a  “z”.  Thus,  the  results  are  the  same:  to  control  a  discrete 
system,  the  discrete  open  loop  poles  must  be  cancelled,  with  new  pole  placement  (and 
thus  system  control)  occurring  through  design  of  the  controller. 
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The  question  that  must  next  be  answered  is  whether  this  same  pole-zero  cancellation 
occurs  in  the  mixed  continuous-discrete  system  of  this  research.  The  following  proof 
shows  that,  indeed,  this  cancellation  occurs.  Given  a  function,  F(s),  the  definition  for  its 
sampled  version  is  given  as  [7:19] 


P*(s)  =  jr'ZF(s  +  jnGi,)  , 

■^s  n=~oo 


[2.38] 


given  the  order  of  the  system’s  numerator  is  less  than  the  order  of  its  denominator,  which 
is  true  for  all  real  systems.  Recall  H*(s),  the  sampled  version  of  H(s) ,  was  found  to  be 


H*is)  =  I; 


Ts  „=^(s  +  yn®,  )D  (s  +  jnca, ) 


[2.39] 


Finding  a  common  denominator  of  the  term  within  the  summation  sign  gives 


oo 

m^n 

oo 

Np(s  +  J  n j^^s)Dp (s  +  ymmj]  > 

m=-oo 

n  b  +  jncOs  )Dp(s  +  j] 


[2.40] 


42 


Assuming  a  similar  form  for  the  sampled  version  of  the  controller  gives 


oo 

m^n 

oo 

I 

Ncis  + 

^  n--oo 

[2.41] 


where  the  controller  includes  no  zero-order  hold. 


The  equation  for  the  discrete-to-continuous  transfer  function  is 


C(.y) 

T,R*is) 


G*(s) 


His) 


l  +  G*is)H*is) 


[2.42] 


Substituting  the  equations  for  *  (.s) ,  H(s),  and  H  *  (s)  into  Equation  2.42  gives 


T. 

\  ^  ) 

^  m^n 

+  ymco,) 

V- 

1 

1, 

_ _ 

■4-1 

V  '  > 

^  m*n  '' 

^  m^n  \ 

^Np(s  +  >(0,)]^  (j  +  jm(A,)D^  (y  +  jnua,) 

n=-<»  m=-«» 

n=-oo 

<  j 

]^  ( j  +  jnco,  (y  +  >co, ) 

n=z-ea 

V  )- 

[2.43] 
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which,  after  algebraic  manipulation  gives 


C(s)  _{l-e  )(iVp  (j))(nM/nl)(nMm2) 
T^R*(s)  sDp{s)(denl  +  den2) 


[2.44] 


where 


m^n 

oo  oo 

numl  jn(0, )  ]][ D, (j  +  jmd)^ )  , 

n=-oo  m^-oo 


[2.45] 


oo 

numl  =  n(^  +  jn({)JDpis  +  jnoij  , 

n-’-oo 


[2.46] 


deni  =  T‘ 


Y 


+  yncoj  n(j  +  ynO),)Z)^(j  +  ynco^) 

n^-oo  ^n=-oo 


,  and  [2.47] 


deni  =  (1  -  €  ■'  0 


V 


X  ^A^+jnfUs) 


X  iVp  ( J  +  ;nco, )  n  (-y  +  ;m(D, ) 


A 


[2.48] 


Notice,  however,  that  the  n  =  0  term  in  the  numl  term  will  cancel  the  sD^  (s)  from  the 
denominator: 
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C(s)  _(l-e  ](Np  (s))(numl)(num3) 

T^R  *  (s)  (deni  +  deni) 


where 


n^O 

oo 

num3=  ^(■y  +  jn(ii^)Dp(s  +  jno)^)  .  [2.50] 

n=-oo 


Thus,  one  set  of  discrete  zeros  will  always  cancel  the  open-loop  continuous  poles  from 
the  plant,  allowing  for  design  of  the  control  for  the  discrete-to-continuous  system. 


2.5  Root  Locus  Development 

The  root  locus  of  a  system  may  be  plotted  in  one  of  two  domains — ^the  continuous  (s-  or 
Laplace  domain)  or  the  discrete  (z-domain).  Because  of  the  peculiar  nature  of  the 
discrete-to-continuous  system  developed  here,  however,  the  root  loci  will  by  nature  be 
“mixed,”  containing  both  continuous  and  discrete  portions.  It  is  advantageous  to  first 
calculate  poles  and  zeros  in  the  appropriate  domain  then  combine  these  poles  and  zeros 
from  both  domains  into  the  Laplace  domain. 
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To  begin  examination  of  the  root  loci  plots,  the  discrete-to-continuous  transfer  function 
previously  introduced  will  be  rewritten  to  reflect  the  domain  in  which  calculations  will 
initially  be  performed: 


Cjs)  ^  C(s) 

TXis)  mz) 


[2.51] 


or 


j-GlUWU)  yGAz)H(.s) 

l+0,(z)H(,z)‘  l+G.(z)H{z) 

The  starred  quantities  are  the  periodic  versions  of  the  sampled  time  function  in  the 
s-domain.  For  example,  the  sampled  form  of  a  function  f(t)  is  /  *(0 ,  which  is  a  series 
of  impulses  at  the  sample  points,  has  a  value  of  zero  between  sample  points,  and  is 
continuous.  Recall  that  H(z)  is  the  z-transform  of  the  plant  Gp(s)  with  a  zero-order 

hold  applied.  Therefore,  H*is)  is  the  periodic  version  of  the  function  H(s)  which 
gives  the  same  values  as  H(z)  at  the  sample  instants. 

Suppose  each  of  the  pieces  of  this  transfer  function  is  split  into  numerator  and 
denominator  terms: 
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G,(Z) = 


G^num(z) 

G^den{z) 


[2.53] 


ZOHjs)  _ 
Ts  ■ 


[2.54] 


Gp(^)  = 


Gpnum(s) 

Gpden(s) 


,  and 


[2.55] 


Hden{z) 


[2.56] 


Substituting  these  back  into  Equations  2.51  and  2.52: 


Gjium{z)\-e~"'^’‘  Gpnum(s) 
Cis)  _  G^deniz)  sT^  Gpden(s) 
TsR(z)  ^  ,  G^num{z)  Hnum(z) 
G^den(z)  Hden{z) 


[2.57] 


Following  algebraic  manipulation  Equation  [2.57]  becomes 


C{s)  _ 

Gjium{z)Hden{z) 

(1  -  )Gpnum(s) 

TMz) 

Gcnum(z)Hnum(z)  +  G^den{z)Hden(z) 

isT,)Gpden(s) 

[2.58] 
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This  equation,  then,  becomes  the  basis  for  determining  pole  and  zero  locations. 

2.5.0  Zeros 

The  zeros  of  Equation  2.58  will  be  located  where  the  numerator  terms  equate  to  zero. 
These  zeros  will  consist  of  four  sets  of  zeros  from  the  following  components  of  the 
discrete-to-continuous  transfer  function:  1)  the  zero-order-hold  numerator,  2)  the 
continuous  plant  numerator,  3)  the  discrete  controller  numerator,  and  4)  the  discrete  plant 
denominator. 

The  zeros  found  from  the  zero-order-hold  numerator  are  defined  by  the  equation 

l-e""^^=0  [2.59] 

in  the  s-plane.  The  solution  to  Equation  2.59  gives 

s  =  0±jna},  n  =  0,1,2,...  [2.60] 

This  corresponds  to  one  zero  at  z  =  1  in  the  z-plane,  which  is  fixed  regardless  of  the 
sample  period,  7^ .  In  the  s*-plane  the  zero-order-hold  produces  zeros  starting  at  j  =  0 

and  repeated  every  to^  in  both  the  positive  and  negative  directions  along  the  imaginary 
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axis.  A  change  in  the  sample  period,  ,  will  only  change  the  spacing  along  the 


imaginary  axis  of  these  zeros. 


The  continuous  plant  numerator  will  introduce  zeros  at  the  solutions  to  the  equation 


Gpnum(s)  =  0 


[2.61] 


Locations  of  these  zeros  in  the  s*-plane  are  fixed,  given  the  equation  of  the  plant  is  fixed. 
These  zeros,  since  they  are  from  the  continuous  plant,  do  not  exist  in  the  z-plane. 


The  discrete  controller  numerator  will  produce  zeros  in  the  z-plane  according  to  the 
solution  of  the  equation 


G^numiz)  =  0 


[2.62] 


These  zeros  may  be  placed  anywhere  in  the  z-plane  according  to  classical  design  theory 
with  the  only  restriction  being  physical  realizability,  which  gives  the  condition  that 


Order  {Gcnum(z))  ^  Order  {G^deniz)) 


[2.63] 
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The  zeros,  once  the  controller  is  designed,  are  fixed  in  the  z-plane  with  no  dependence  on 
the  sample  period,  .  In  the  s*-plane  the  zeros  can  be  found  by  transforming  the  z-plane 
zeros  by 


s  = 


[2.64] 


giving  their  location  in  the  principal  band.  This  transformation  will  be  explained  in  a 
later  paragraph.  These  zeros  are  then  repeated  at  jn(o^  intervals  for  n  =  0,  ±  1,  ±  2, . . . 
Sample  period  only  changes  the  spacing  of  these  zeros  in  the  imaginary  axis  direction. 

The  zeros  found  from  the  discrete  plant  denominator  are  found  from  the  solution  of  the 
equation 


Hden(z)  =  0  [2.65] 

In  the  z-plane  these  correspond  to  the  location  of  the  open-loop  poles  of  the  plant 
transformed  according  to 


z  = 


[2.66] 
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so  their  locations  are  sample  period  dependent.  Just  like  the  controller  zeros,  they  are 
transformed  to  the  s*-plane  and  repeated  in  the  imaginary  axis  direction  at  jnOi^ 

intervals.  Because  they  are  sample  period  dependent  in  the  z-plane,  they  are  sample 
period  dependent  in  the  s-plane,  in  addition  to  the  sample  period  dependence  of  the 
spacing  of  repeated  zeros  in  the  imaginary  axis  direction.  Notice  these  zeros  are  located 
at  precisely  the  same  location  they  were  at  before  the  z-transformation,  only  now  they  are 
repeated  in  every  band. 

A  summary  of  the  location  of  the  zeros  is  presented  in  table  2. 1  below. 


Table  2.1  Summary  of  Zero  Locations 


Zero  name 

Found  from. . . 

Z-plane  T* 
Dependence? 

S*-Plane  Ts 
Dependence? 

1.  ZOH  Numerator 

=0 

No 

Yes 

2.  Gp(s)  Numerator 

N/A 

No 

3.  Gc(z)  Numerator 

No 

Yes 

4.  H(z)  Denominator 

Hdeniz)  =  0 

Yes 

Yes 

2.5.1  Poles 

The  poles  of  Equation  2.58  will  be  located  where  the  denominator  terms  equate  to  zero. 
These  will  consist  of  three  sets  of  poles  from  the  following  components  of  the  discrete-to- 
continuous  transfer  function:  1)  the  zero-order-hold  denominator,  2)  the  continuous  plant 
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denominator,  and  3)  the  combined  term  involving  both  numerator  and  denominator  terms 
from  both  the  discrete  controller  and  discrete  plant. 

The  single  pole  from  the  zero-order-hold  denominator  is  found  from  solving  the  equation 

sT,=0  [2.67] 


giving 


^  =  0  [2.68] 

This  pole  in  the  s*-plane  is  located  at  the  origin  and  cancels  the  zero  from  the  zero-order- 
hold  located  there.  It  is  sample  period,  ,  independent  in  the  s*-plane,  and  has  no  real 
meaning  by  transforming  to  the  z-plane. 

The  locations  of  the  continuous  plant  denominator  poles  come  from  the  solution  of  the 
equation 


Gpden(s)  =  0 


[2.69] 


The  locations  of  these  poles  are  the  locations  of  the  open  loop  plant  poles,  and  the 
locations  do  not  move,  given  the  equation  of  the  plant  is  fixed. 

Poles  from  the  combined  term  are  found  by  solving  the  equation 

G^num(z)Hnum(z)  +  G^den(z)Hden{z)  =  0  [2.70] 

Locations  of  these  poles  in  the  z-plane  are  dependent  both  upon  the  controller  design  as 
well  as  the  sample  period,  ,  because  of  the  discretization  of  the  plant,  Gp(s) ,  which 

gives  the  locations  of  the  Hnum(z)  and  Hden(z)  terms.  The  locations  in  the  s*-plane 
are  determined  by  the  transform  of  Equation  2.70  to  the  s*-plane  using  the  inverse  of 
Equation  2.67,  and  repeated  in  the  imaginary  axis  direction  at  jnco^  intervals.  Because 

they  are  sample  period  dependent  in  the  z-plane,  they  are  sample  period  dependent  in  the 
s*-plane,  in  addition  to  the  sample  period  dependence  of  the  spacing  of  repeated  poles  in 
the  imaginary  axis  direction. 

A  summary  of  the  location  of  the  poles  is  presented  in  table  2.2  below. 
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Table  2.2  Summary  of  Pole  Locations 


Pole  name 

Found  from. . . 

Z-plane  Ts 
Dependence? 

S*-Plane  Ts 
Dependence? 

1.  ZOH 
Denominator 

11 

o 

N/A 

No 

2.  Gp(s) 
Denominator 

Gpden{s)  =  0 

N/A 

No 

3.  Combined 
Term 

G^den(z)Hden(z) 

+  G^num(z)Hnum(z)  =  0 

Yes 

Yes 

2.6  Chapter  Two  Summary 

In  this  chapter,  the  s*-plane  was  introduced.  The  s*-plane  is  the  plane  that  allows 
analysis  of  the  discrete  and  continuous  pieces  of  the  sampled-data  system  together.  The 
model  of  the  discrete-to-continuous  system  and  its  related  equations  were  developed.  The 
model  consists  of  discrete  and  continuous  portions.  In  the  limits  as  the  sample  period 
goes  to  zero  and  infinity  the  model  reduces  to  its  continuous  and  discrete  counterparts, 
respectively.  A  set  of  the  discrete  closed-loop  controller  zeros  was  determined  to  cancel 
the  continuous  open-loop  poles,  allowing  for  control  of  the  discrete-to-continuous  system. 
The  origins  of  the  mixed  system  zeros  and  poles  and  the  effects  of  the  sample  period  on 
these  roots  were  examined  and  determined.  With  the  theoretical  development  complete, 
examination  of  the  effects  of  aliasing  on  the  design  of  the  discrete-to-continuous  system 
can  now  begin. 
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Chapter  Three.  Second-Order  System  Design  With 

Aliasing  -  Proportional  Control 


3.0  Introduction 

With  the  tools  and  model  developed  for  the  discrete-to-continuous  system,  design  issues 
will  now  be  addressed.  First,  the  generalized  root  locus  for  the  system  is  developed  and 
the  locations  of  the  poles  in  the  gain  limited  is  examined.  The  second-order  plant  is 
introduced  and  its  z-transform  performed.  The  effect  of  the  sample  time  on  the  zero-gain 
root  locus  is  determined.  Using  a  simple  proportional  gain  for  the  discrete  controller 
(P-control),  the  effects  of  different  sampling  rates  that  allow  aliasing  of  the  system  are 
examined.  Comparison  between  root  loci,  frequency  response  functions,  impulse  and 
step  responses,  and  responses  to  sinusoidal  inputs  is  made  to  a  system  with  the  same 
controller  gain  sampled  at  a  fast  sample  rate  to  avoid  aliasing.  A  sampling  rate  equal  to 
the  damped  natural  frequency  of  the  second-order  system  is  found  to  give  excellent 
response  and  an  optimal  gain  for  this  system  with  P-control  is  found.  Robustness  of  this 
system  is  investigated. 
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3.1  Root  Locus  Preliminaries 


With  the  additional  allowance  of  aliasing  in  the  system,  the  sample  period  or  time  as  it 

introduces  aliasing  now  becomes  a  design  parameter.  One  must  also  examine  the  effects 

/ 

of  the  design  in  the  s*-plane,  the  plane  in  which  poles  and  zeros  incorporate  both  the 
continuous  plant  and  the  discrete  controller. 


For  convenience  in  design,  Equation  2.58,  is  repeated: 


C(s)  _ 

Gjium{  z)Hden{  z) 

(l—e  ^^’)Gpnum(s) 

TMz) 

G^num(z)Hnum(z)  ■^^den(z)Hden(z) 

(sT^)Gpden(s) 

Now  letting  the  controller  be  written  as  in  Equation  2.53,  but  with  a  controller  gain,  K: 


G,den{z) 


[3.2] 


Grouping  all  the  terms  in  Equation  3.1  which  are  functions  of  z  and  substituting 


Equation  3.2,  define  F{z) : 


F(z)  = 


_ KGjium{z)Hden{z) _ 

KG^num(z)Hnum(z)  +  G^den{z)Hden{z) 


[3.3] 
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Finally,  recalling  the  zero-order-hold  times  the  plant  is  defined  as  H(s) ,  Equation  3.1 
can  now  be  written  as 


C(s) 

TMz) 


=  Fiz) 


Hjs) 

Ts 


[3.4] 


Using  Equations  3.2  and  3.3,  the  controller  design  will  be  accomplished  in  the  digital  or 
z-domain.  However,  the  effects  on  the  overall  discrete-to-continuous  transfer  function 
must  be  examined  as  well.  Therefore,  the  effects  of  z-plane  controller  design  will  be 
combined  with  the  continuous  plant  effects  in  the  s*-plane. 


3.2  Root  Loci  AT  Gain  Limits 

Since  the  root  loci  in  both  the  z-  and  s*-planes  will  be  plotted,  determination  of  trends  as 
the  gain  of  Equation  3.3  approaches  zero  and  infinity  needs  to  be  made. 

For  the  z-domain,  in  the  limit  as  K  approaches  zero,  the  poles  of  F(z)  will  approach  the 
roots  of  G^deniz)  and  Hden{z) .  Likewise,  in  the  limit  as  K  approaches  infinity,  the 
poles  of  F{z)  will  approach  the  roots  of  G^num(z)  and  Hnum(z) .  As  in  either  the 

continuous  or  discrete  cases,  then,  the  poles  of  the  root  locus  begin  at  the  open-loop  pole 
locations  and  migrate  toward  the  open-loop  zeros  of  the  system.  Notice,  however,  that 


this  migration  of  the  discrete  poles  is  not  toward  the  open-loop  zeros  of  the  continuous 
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system  Hnum(s) ,  but  toward  the  open-loop  zeros  of  the  discrete  or  sampled  version  of 
the  plant  Hnum{z) . 

These  results  transform  to  the  s*-plane  as  well,  using  the  transformation  of 
Equation  2.64.  In  other  words,  the  root  locus  of  the  closed-loop  system  begins  at  the 
open-loop  pole  locations,  found  from  G^*den(s)  and  H*den(s) ,  and  migrate  toward 

the  open-loop  zeros  found  from  *num(s)  and  H*num{s) .  In  this  case,  however,  the 

location  of  the  poles  from  Hden(s)  and  H*  denis)  are  the  same,  even  though  the  zeros 
of  Hnumis)  and  H  *  numis)  are  different. 

The  root  locus  in  the  s*-plane  will  therefore  include  the  following  poles  and  zeros: 

1 .  The  poles  and  zeros  from  the  continuous  plant  G^  ( j?) .  These  poles  and  zeros 

are  stationary  and  will  not  move  based  on  either  controller  gain  or  sample 
period  selections. 

2.  The  zeros  from  the  zero-order-hold.  These  zeros  will  not  move  based  on 
controller  gain,  but  will  repeat  at  intervals  dependent  upon  as  described  in 
chapter  three. 
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3.  The  single  pole  of  the  zero-order-hold.  It  remains  fixed  at  =  0 ,  but  is 
always  cancelled  by  one  of  the  zeros  of  the  zero-order  hold  which  is  at  the 
same  location. 

4.  The  zeros  from  the  numerator  of  the  mixed  term  of  the  controller  F(z) 
transformed  to  the  s*-plane  (F*(s) ).  These  zeros  will  not  move  based  on 
controller  gain.  In  the  s*-plane,  they  are  located  at  the  location  of  the  open- 
loop  poles  of  H(s)  and  repeated  at  intervals  dependent  upon  7^  as  described 
in  chapter  three. 

5.  The  poles  from  the  denominator  of  the  mixed  term  of  the  controller  F(z) 
transformed  to  the  s*-plane  (F*(s) ).  These  poles  begin  at  the  location  of  the 
open-loop  poles  and  move  toward  the  location  of  the  open-loop  zeros  of 
H*(s),  which,  again,  is  different  than  H{s) . 


3.3  The  Second-Order  Plant 

Before  determining  the  effect  of  sample  period  on  the  discrete-to-continuous  system  of 
the  model  developed  in  chapter  three,  the  characteristics  of  a  simple  second-order  plant 
will  be  examined. 
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Consider  a  plant  that  is  lightly  damped  with  two  poles  of  the  form 


G,(s)  = 


+2^(D„i:  +  tO„^ 


[3.4] 


Easily  recognizable  are  the  locations  of  the  poles  in  the  s-plane: 


[3.5] 


Note  that  for  damping  less  than  ^  =  1.0  these  poles  are  a  complex  conjugate  pair. 
Labeling  these  poles  -  and  -  ,  Equation  3.4  can  be  rewritten  as 


G,(s)  = 


(s  +  PiXs+p^) 


[3.6] 


Using  partial  fraction  expansion  on  Equation  3.6,  divided  by  s  in  preparation  for  its 
z-transform  with  the  zero-order  hold  applied,  gives 


G,(s) 

s 


s  +  Pi 


+ 


«2 

S+Plj 


[3.7] 
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where 


1 

^0  ~  ’ 
P1P2 


[3.8] 


«i  = 


1 


Px^Pi-Pi) 


,  and 


[3.9] 


1 


<h=- 


PiiPi  -  Pi) 


[3.10] 


Taking  the  z-transform  of  the  plant,  including  the  zero-order  hold  gives 


z{g^  (5)}=  (1  -  z-‘  ^ = (1  -  z“‘  )z|co; 


f 


^0  , 


a. 


—  + - i— + 

5  S+p^  S  +  P2 


z{g,w}=(i-z-‘)(o; 


an 


a, 


(Xr. 


^  ^  2  _  0  PiTj  ^  ^  ^  Pl^s  ^  ^ 


[3.11] 


which,  following  algebraic  manipulation,  becomes 


z{g^(s)}=H{z)  = 


CjZ  +  C^ 

(z+b^Xz  +  b^) 


[3.12] 
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where 


b,=-e-^'‘^‘  ,  [3.13] 

b2=-e~'’^^-  ,  [3.14] 

Ci=®/(«o-«A-«2^2).  and  [3.15] 

C2  =  Ci)„^(ao£»^  -  a,i»2  -  ajbi)  .  [3. 16] 


This  system,  then,  becomes  the  basis  for  the  designs  to  follow. 


3.4  Design  of  A  Proportional  Controller  Using  Sample  Period 
As  A  Design  Parameter 

Since  the  system  consists  of  a  mixture  of  continuous  and  discrete  parts,  it  is  advantageous 
to  determine  the  effect  of  sample  period  Ts  on  the  location  of  the  poles  and  zeros  in  the  z- 
plane,  and  the  effect  of  these  locations  on  the  response  of  the  system.  Thus,  the  sample 
period  ,  or  sample  frequency  (O^ ,  becomes  a  design  parameter. 
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For  the  design  to  follow,  the  controller  will  be  strictly  proportional  control, 


G,(z)  =  G,Hs)-K^ 


[3.17] 


Solving  for  F(z) ,  the  combined  term  of  Equation  3.3, 


F(z)  = 


Kpjz+biXz+b^) 
Kp(CiZ  +  C2)  +  (z+bi)(z  +  b2) 


[3.18] 


where  and  C2  are  as  defined  in  Equations  3.13  through  3.16.  After  algebraic 

manipulation,  this  equation  becomes 


z^  +  djZ  +  ^2 


[3.19] 


where 


=  KpC^  +b^+b2  ,  and 


[3.20] 


d2  =  KpC2  +  b^b2  . 


[3.21] 
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When  the  gain  is  set  to  zero  the  poles  of  Equation  3. 19  are  located  at  -b^  and  -^2  > 

which  also  correspond  to  the  locations  of  the  open-loop  poles,  were  they  transformed  to 
the  z-plane.  Plotting  the  location  of  these  poles,  referred  to  as  zero-gain  poles,  as  sample 
period  7^  is  varied  produces  Figure  3.1  below.  For  purposes  of  non-dimensionalization 

let  T,  =  TjCO„  and  plot  for  changing  .  The  value  of  damping  was  set  to  ^  =  0.1 . 


Figure  3.1  Root  Locus  of  Zero-Gain  Poles  in  Z-Domain 
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In  the  z-plane,  as  increases  the  locations  of  the  zero-gain  poles  move  in  a  logarithmic 

spiral  inward  from  z  =  1  at  =  0  toward  z  =  0  at  =  <» .  Li  the  s*-plane,  a 

relationship  exists  between  the  crossings  of  the  real  axis  and  the  location  of  the  infinitely 
repeated  plant  poles  within  the  bands.  This  relationship  will  be  explored  further  later. 

Li  examining  Figure  3.1,  one  may  question  if,  given  the  need  for  allowing  aliasing,  there 
exist  locations  within  the  z-plane  where  the  discrete-to-continuous  system  response  is 
more  desirable  than  other  locations.  In  other  words,  given  aliasing,  can  the  designer 
select  a  sample  period  that  will  produce  better  results  than  by  selecting  some  other 
sample  period?  To  examine  this,  we  will  look  at  several  classic  tools  used  to  evaluate 
designs,  including  root  loci  in  both  the  z-  and  s*-planes,  magnitude  and  phase  plots  of 
frequency  response  functions,  and  continuous  and  discrete  responses  to  an  impulse,  step, 
and  sinusoidally  varying  inputs.  The  proportional  controller  of  Equation  3.17  will  be 
employed  with  gain  .K'p  =0.15  to  ensure  stability.  Examination  of  this  simple  system 
will  aid  in  seeing  some  trends  otherwise  not  seen  as  clearly. 

To  non-dimensionalize  the  plots  to  follow,  define  the  frequency  ratio  (FR)  as 

FR  =  -^  ,  [3.22] 
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where  is  the  damped  natural  frequency: 


[3.23] 


For  comparison  purposes  of  the  aliased  system.  Figures  3.2  through  3.5  are  plots  of  the 
system  for  FR  =  50,  where  the  discrete  portion  of  Equation  3. 19  is  basically 
continuous — i.e.,  no  aliasing  is  present.  Notice  from  Figure  3.2  the  limited  amount  of 


Z-Plane  Root  Locus;  FR  =  50  Kp  =  0.15 


S*-Plane  Root  Locus;  FR  =  50  Kp  =  0.15 


Figure  3.2  Root  Loci  for  FR  =  50,  Kp  =  0.15 
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Frequency  Ratio  Ws/Wd  =  50;  Kp  =  0.15 


1 


Response  to  sjn(0.05*t)  input;  FR  =  50;  Kp  =  0.15 


Figure  3.5  Response  to  Sinusoidal  Inputs  for  FR  =  50,  Kp  =  0.15 


gain  before  instability  occurs.  Figure  3.3  shows  the  classical  frequency  response  function 
plot  of  a  second-order  system,  with  a  magnitude  peak  at  the  plant  natural  frequency  along 
with  the  phase  passing  through  -  90'  at  that  point,  as  well.  The  impulse  and  step 
responses  of  the  system  are  shown  in  Figure  3.4,  showing  the  characteristics  of  the  lightly 
damped  system.  Finally,  Figure  3.5  shows  how  well  the  closed-loop  system  responds  to 
the  sinusoidal  inputs  shown  as  dashed  lines. 
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Figures  3.6  through  3.9  show  the  root  loci  for  varying  FR,  and  thus  varying  .  These 


correspond  to  frequency  ratios  of  — ,  — ,  — ,  and  — ,  which  are  even  steps  of 

^  ^  1.0  1|  If  2.0 

frequency  ratios  around  the  second  loop  of  Figure  3.1. 

The  symbols  for  these  root  loci  are  as  follows.  The  triangle  and  plus  symbols  represent 
the  zeros  and  pole,  respectively,  of  the  zero-order  hold.  The  square  and  asterisks 
represent  the  zeros  and  poles,  respectively,  of  the  continuous  plant,  (s) .  Although  in 

reality  these  do  not  exist  in  the  z-plane,  they  are  shown  there  for  clarity.  The  diamond 
represents  the  zero  of  the  H(z) .  Again,  this  term  does  not  exist  in  the  s*-plane,  but  is 
shown  there  for  clarity.  The  circles  and  “x’s”  respectively  represent  the  zeros  and  poles 
at  the  given  gain  of  the  mixed  term  F (z)  or  F  *  (j) ,  which  includes  the  controller  and 
the  discretized  plant  with  the  zero-order  hold.  The  solid  line  represents  the  locus  of  the 
poles  of  the  mixed  term  for  all  values  of  .  Finally,  in  the  s*-plane  the  dashed  lines 
represent  the  “band”  boundaries,  which  equate  to  multiples  of  the  folding  frequency. 

Thus,  to  -I-—  represents  the  primary  band,  to  and  +—  to 

2  2  2  2  2  2 

represent  the  secondary  band,  and  so  on. 

With  the  sampling  rate  such  that  the  frequency  ratio  is  2.0  the  damped  natural  frequency 
is  located  at  the  folding  frequency,  and  the  entire  z-plane  root  locus  lies  on  the  negative 
real  axis.  The  plots  of  the  system  root  loci  for  this  frequency  ratio  are  shown  in 
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Figure  3.6.  A  gain  of  =  0.15  was  selected  to  ensure  stability,  and  one  can  see  that, 

although  they  are  close  to  the  stability  boundaries,  the  poles  of  Figure  3.6  do  indeed  lie 
within  the  stability  boundaries  of  the  unit  circle  in  the  z-plane  and  the  negative  left-half 
plane  in  the  s*-plane. 

The  entire  root  locus  lies  entirely  upon  the  real  axis  in  the  z-plane.  Although  not  easily 
discemable,  the  zero  from  H{z) ,  and  therefore  H*{s),  lie  directly  on  top  of  the  zeros  of 

F(z)  and  the  open-loop  poles  of  Gp{s) .  Thus,  one  of  the  poles  of  F{z)  does  not  move 

at  all,  a  phenomenon  that  will  be  discussed  later.  The  remaining  pole  moves  from  the 
closed-loop  zero  toward  infinity,  as  predicted  by  gain  limit  analysis. 

Additionally,  one  can  graphically  see  here,  as  well  as  in  Figures  3.7  through  3.9,  the 
proof  that  one  set  of  the  closed-loop  zeros  cancels  the  open-loop  poles,  as  proved 
mathematically  in  chapter  three  (Equation  2.49),  in  both  the  z-plane  and  s*-plane  root 
loci.  In  the  s*-plane,  the  open-loop  poles,  closed-loop  zeros,  and  the  zeros  of  H{z)  all 
lie  on  the  band  boundaries:  the  frequencies  corresponding  to  multiples  of  the  folding 
frequency.  A  frequency  ratio  of  FR  =  2.0 ,  then,  causes  the  entire  root  locus  to  lie  on  the 
real  axis  of  the  z-plane  and  on  band  boundaries  of  the  s*-plane. 
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Figure  3.6  Root  Locus  -  FR  =  2.0,  Kp  =  0. 15 


71 


One  last  thing  to  recognize  about  Figure  3.6  is  the  limited  amount  of  gain  allowed  before 
the  stability  boundary  is  reached.  With  this  small  gain,  however,  the  closed-loop  pole  has 
moved  quite  a  distance  from  the  closed-loop  zeros.  Both  these  facts  speak  to  the 
robustness  of  this  design. 

Figures  3.7  and  3.8  show  the  root  loci  for  frequency  ratios  of  1.5  and  1.2,  values  of 
frequency  ratios  that  “move”  the  location  of  the  zero-gain  poles  around  the  second  loop 
of  Figure  3.1.  Notice  in  the  s*-plane  the  open-loop  poles  fall  between  the  boundaries  of 
the  secondary  bands,  indicating  that  aliasing  is  occurring.  If  there  were  no  aliasing,  these 
poles  would  remain  within  the  primary  band.  As  further  and  further  aliasing  occurs, 
these  poles  will  move  to  bands  further  and  further  from  the  primary  band.  Moreover,  as 
the  frequency  ratio  decreases  (sample  period  increases)  the  size  of  all  the  bands 
decreases,  in  effect  “moving”  the  open-loop  poles,  and  hence  the  zero-gain  poles,  to 
different  locations  within  the  bands.  One  may  note  these  poles  are  moving  toward  the 
center  of  the  band,  which  corresponds  to  multiples  of  the  sampling  frequency. 

Further  examination  of  the  s*-plane  root  loci  show  that,  indeed,  the  closed-loop  poles 
begin  at  the  location  of  the  closed-loop  zeros — ^the  zeros  of  F(^z)  are  correspondingly  the 
open-loop  poles  of  (s) — and  migrate  toward  the  open-loop  zeros  of  - H(z) . 


Another  effect  of  increasing  the  sampling  period  identifiable  from  the  Figures  is  that  of 
the  decreasing  effect  of  the  gain  ,  since  the  poles,  given  the  same  gain  of  =  0.15 , 
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Figure  3.8  Root  Locus  -  FR  =  1 .2,  Kp  =  0. 1 5 
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move  less  and  less  distance  from  the  closed-loop  zeros.  This  would  tend  to  indicate  an 
increased  robustness  with  increased  aliasing!  This  fact  will  be  explored  further  later. 

Finally,  Figure  3.9  shows  the  root  loci  for  a  frequency  ratio  of  1.0.  Notice  the  similarities 
between  these  root  loci  and  the  root  loci  of  Figure  3.6.  The  entire  root  locus  lies  on  the 
real  axis  in  the  z-plane.  Here  again,  one  of  the  closed-loop  poles  lies  directly  on  top  of 
one  of  the  closed  loop  zeros,  and  remains  at  the  location  of  the  open-loop  zero  of  H{z) . 

The  mathematics  behind  this  stationary  pole  follow:  recall  the  definition  of  H{z)  of 
Equation  3.12,  repeated  here  for  clarity: 


H{z)  = 


iz  +  b^){z■^b^) 


[3.24] 


where  ,  ^2*  Cj  are  as  defined  in  Equations  3.13  through  3.16.  At  the  frequency 

2  2 

ratios  of  interest,  FR= —  and  FR  =  — ,  the  entire  root  locus  lies  on  the  real  axis. 

1.0  2.0 


In  this  case. 


b^=b^=b  . 


[3.25] 
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The  zero  of  H{z)  is  found  from  setting  the  numerator  term  of  Equation  3.24  equal  to 


zero: 


zero  =  ■ 


-c 


2  _ 


=  0: 


{a^b + -  ajb^) 
(o/(a,ij  +  a2fe-ao) 


[3.26] 


Equation  3.26  is  easily  simplified  using  a  relationship  found  from  Equations  3.8  through 
3.10: 


aj+a2=-ao  •  [3.27] 

Substituting  Equation  3.27  into  Equation  3.26  and  canceling  gives: 

_  [3.28] 

a^b  +  -Oq  -  ajb  - 

which  from  Equation  3.19  is  the  value  of  the  closed-loop  zeros.  Hence,  the  open-loop 
zero  of  H{z)  is  located  at  the  same  point  as  the  closed-loop  zeros  and  the  open-loop 
poles,  and  one  of  the  closed-loop  poles  remains  at  this  point. 

The  motion  of  the  other  closed-loop  poles,  however,  is  somewhat  different  in  this  case 
since  the  closed-loop  zeros  are  located  together  on  the  positive  real  axis.  Because  of  this, 
the  second  closed-loop  pole  will  travel  through  the  point  z  =  0 ,  which  corresponds  to  the 
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pole  moving  through  infinity  in  the  s*-plane.  At  that  point,  a  shift  occurs  because  of  the 
nature  of  the  natural  logarithm  of  a  complex  number  (chapter  two).  As  z  approaches 
z  =  0  from  the  negative  side  of  the  real  axis,  its  natural  logarithm  has  an  argument  equal 
to  71 .  Once  z  crosses  to  the  positive  side  of  the  real  axis,  the  argument  of  its  natural 
logarithm  equals  0.  This  fact  is  seen  in  the  s*-plane  as  the  pole  travels  toward  -  «>  along 
the  band  boundaries,  corresponding  to  an  argument  of  n,  then  “returns”  from  -  «>  along 
the  mid-band  line,  corresponding  to  an  argument  of  0.  Hence,  in  the  s*-plane,  the  second 
pole  travels  toward  -oo  along  the  mid-band  line  then  returns  along  the  band  boundary. 
This  mid-band  to  band  boundary  “jumping”  will  occur  anytime  the  zero  of  H(z)  lies  on 
the  positive  real  axis  of  the  z-plane. 


9 

Finally,  recall  the  s*-plane  root  loci  for  the  system  in  which  FR  =  —  was  entirely  on  the 

1.0 

line  separating  bands,  the  lines  corresponding  to  multiples  of  the  folding  frequency. 
Similarly  here  the  open-loop  poles,  the  closed-loop  zeros,  and  the  zero  of  H*  (s)  with  its 
accompanying  closed-loop  pole  all  lie  on  the  mid-band  line:  the  line  corresponding  to 
multiples  of  the  sampling  frequency.  One  can  see  that  unique  things  occur  when  the  root 
loci  lie  on  these  lines. 

These  unique  things  occur  at  specific  frequency  ratios.  The  frequency  ratios  of  interest 
correspond  to  the  real  axis  crossings  of  the  zero-gain  root  locus  of  Figure  3.1.  The 
crossings  of  the  real  axis  occur  at  multiples  of  the  frequency  ratio  defined  by  the 
following  equation: 
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FR 


"crossing 


n 


n  =  1,2,3,... 


[3.29] 


In  the  z-plane,  if  n  of  Equation  3.29  is  odd  the  root  locus  lies  entirely  on  the  negative 
real  axis.  The  zero  of  H{z)  lies  on  top  of  the  open-loop  poles  and  the  closed-loop  zeros, 
which  are  all  located  at  the  same  spot.  One  of  the  closed-loop  poles  also  lies  at  this  point, 
since  it  starts  at  the  closed-loop  zero  and  travels  to  the  open-loop  zero  of  H{z) . 
Additionally,  the  other  pole  travels  along  the  negative  real  axis  toward  infinity,  the  other 
zero  of  H{z) .  Since  both  closed-loop  zeros  are  located  on  the  negative  real  axis,  neither 
pass  through  the  point  z  =  0 . 

If  n  is  even,  the  entire  root  locus  lies  partially  on  the  positive  and  partially  on  the  negative 
real  axis.  Again,  the  zero  of  (z)  lies  on  top  of  the  open-loop  poles  and  the  closed-loop 
zeros,  and,  with  one  of  the  closed  loop  poles,  all  are  located  at  the  same  spot.  The  other 
pole  travels  along  the  real  axis  toward  infinity,  passing  through  the  point  z  =  0 . 

These  same  trends  are  evident  in  the  s*-plane.  However,  if  n  of  Equation  3.29  is  odd, 
the  open-loop  poles,  closed-loop  zeros,  and  closed-loop  poles  are  all  located  on  band 
boundaries,  corresponding  to  integer  multiples  of  the  folding  frequency.  If  n  is  even, 
however,  the  open-loop  poles  and  closed-loop  zeros  are  located  on  the  mid-band  line, 
corresponding  to  integer  multiples  of  the  sampling  frequency.  The  one  traveling  pole 
will  also  move  toward  infinity  along  this  mid-band  line,  but  return  along  the  band 
boundary  line. 
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Movement  of  the  open-loop  poles  within  the  bands  of  the  s*-plane  is  also  dependent  upon 
the  frequency  ratio.  This  “movement”  is  a  measure  of  how  much  the  poles  are  aliased. 
Traditionally,  aliasing  has  been  thought  of  as  the  folding  and  addition  of  components  of 
the  frequency  response  when  these  components  are  not  equal  to  zero  outside  of  the 
folding  frequency.  Here,  however,  we  wish  to  discuss  the  aliasing  of  poles,  or  the 
movement  outside  of  the  primary  band  of  the  poles  causing  the  response.  This  type  of 
aliasing  will  hereafter  be  referred  to  as  “pole-aliasing.”  When  thinking  of  pole-aliasing, 
think  of  the  poles  transformed  from  the  z-plane  to  the  s*-plane  and  repeated  infinity 
through  all  bands. 

Little  aliasing  as  defined  traditionally  is  present  for  frequency  ratios  much  greater  than 
FR  =  2.0 ,  the  aliasing  boundary  according  to  Nyquist.  At  this  boundary,  however, 
aliasing  is  already  occurring,  unless  the  signal  is  a  pure  sinusoid.  Before  this  boundary, 
pole-aliasing  does  not  occur.  These  frequency  ratios  correspond  to  the  first  half-loop  of 
Figure  3.1  from  z  =  1  to  the  point  it  meets  the  negative  real  axis.  The  open-loop  poles  are 
located  within  the  primary  band  for  these  frequency  ratios  as  well.  At  FR  =  2.0 ,  the 
open-loop  poles  lie  on  the  boundary  between  the  primary  and  secondary  bands.  This  is 
the  boundary  for  pole-aliasing. 

As  the  frequency  ratio  is  decreased,  the  open-loop  poles  move  farther  into  the  secondary 
band.  These  changes  correspond  to  traveling  along  the  second  loop  of  Figure  3.1  until 
the  positive  real  axis  is  reached  at  a  frequency  ratio  of  FR  =  1.0.  At  this  point,  the  open- 
loop  poles  are  located  in  the  exact  middle  of  the  frequency  band.  Now  as  the  frequency 
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ratio  is  decreased,  corresponding  to  the  third  loop  of  Figure  3.1,  the  poles  travel  again 
toward  the  band  boundary  and  will  reach  it  when  the  frequency  ratio  is  FR  =  0.67 .  The 
special  characteristics  of  the  response  at  the  frequency  ratios  of  Equation  3.29  will 
continue  to  be  seen  through  the  remaining  plots  and  examined  in  greater  detail  below. 

Figures  3.10  through  3.13  show  the  magnitude  and  phase  plots  of  the  discrete-to- 
continuous  frequency  response  functions  of  the  system  for  the  same  frequency  ratios. 

The  effects  of  aliasing  are  clearly  discemable  in  these  plots.  Figure  3.10,  corresponding 
to  a  frequency  ratio  of  2.0,  shows  the  dominant  peak  of  the  damped  natural  frequency 
located  at  the  folding  frequency.  The  effect  of  the  zero-order  hold  can  be  seen  at 

multiples  of  the  sampling  frequency,  —  =  2,4 . Notice  that  as  the  frequency  ratio 

decreases,  the  aliasing  becomes  prevalent  in  the  additional  peaks  seen  in  the  magnitude 
plots.  Finally,  in  Figure  3. 13  see  the  disappearance  of  any  peak  in  the  magnitude  plot  of 
the  system  with  frequency  ratio  equal  to  1 .0.  It  appears  that  by  selecting  a  correct 
sampling  frequency,  a  response  cancellation  occurs.  This  is  indeed  the  case,  as  will  be 
shown  later. 

Finally,  examination  of  the  phase  plots  of  Figures  3.10  through  3.13  shows  that  increased 
aliasing  tends  to  linearize  the  phase  plots.  The  phase  plot  of  Figure  3. 13  seems  much 
more  linear  than  the  phase  plot  of  Figure  3.10,  although  this  is  also  a  function  of  the  gain 
margin  differences. 
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Phase  (deg)  Magnitude  (dB) 


Frequency  Ratio  Ws/Wd  =  1.2;  Kp=:0.15 


Figure  3.12  Frequency  Response  Function  Plot  -  FR  =  1.2,  Kp  =  0.15 
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Figures  3.14  through  3.17  show  the  impulse  and  the  step  responses  for  the  frequency  ratio 
changes  of  the  systems  represented  by  root  loci  in  Figures  3.10  through  3.13.  Here  the 
solid  line  represents  the  response  of  the  continuous  plant  (or  continuous  output).  The 
asterisks  represent  the  discrete  response,  or  the  response  the  controller  sees.  Responses 
here  show  the  same  trends  visible  in  the  root  loci  and  frequency  response  function  plots. 

Figure  3.14  shows  the  responses  of  the  system  with  a  frequency  ratio  of  2.0.  Based  on  a 
location  of  the  poles  in  the  z-plane  root  locus  the  response  of  the  discrete  system  should 
be  a  lightly  damped  sinusoidal  response,  which  is  indeed  the  case.  For  both  the  impulse 
and  step  responses,  the  response  is  much  worse  than  the  “continuous”  responses  of 
Figure  3.3.  The  magnitude  of  the  initial  pulse  of  the  impulse  response  is  greater  by  a 
factor  of  almost  20.  Although  the  magnitude  of  the  step  response  is  only  slightly  larger, 
the  number  of  oscillations  to  obtain  its  final  value  is  very  much  greater. 

As  the  frequency  ratio  decreases  in  Figures  3.15  through  3.17,  both  impulse  and  step 
responses  improve.  Note  the  magnitudes  of  the  responses’  initial  “pulses”  remain  the 
same,  but  the  normalized  settling  time  improves  greatly.  Additionally,  discrete  responses 
(dashed  lines)  are  as  would  be  predicted  given  the  pole  locations  in  the  z-plane 
[33:351-353].  The  continuous  responses,  however,  are  those  of  the  continuous 
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Response  Response 


Impulse  Response  at  Ws/Wd  =  1 ;  Kp  =  0.15 


Step  Response  at  Ws/Wd  =  1 ;  Kp  =  0.1 5 


Figure  3.17  Impulse  and  Step  Responses  -  FR  =  1 .0,  Kp  =  0. 15 
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plant.  In  both  cases,  the  continuous  response  is  that  of  step  inputs  at  the  sample  periods 
with  the  added  initial  conditions  caused  because  of  previous  plant  response. 

The  continuous  response,  then,  is  the  open-loop  response  of  the  plant  to  a  step  input  with 
the  given  initial  conditions.  This  open-loop  response  is  the  damped  sinusoidal  response 
given  from  the  open-loop  poles  located  where  they  are  in  the  s*-plane.  How  can  this  be, 
given  these  poles  are  essentially  cancelled  by  one  set  of  closed-loop  zeros? 

This  phenomenon  appears  to  be  similar  to  Gibb’s  phenomenon  [35:740-741].  In  taking 
an  infinite  Fourier  series  of  a  function,  a  “ringing”  appears  at  the  discontinuity  produced 
by  truncating  the  series.  Taken  one  step  farther,  if  one  term  were  eliminated  from  this 
series,  the  result  would  be  a  sinusoid  of  the  frequency  of  the  eliminated  term.  This  is 
illustrated  for  the  following.  The  Fourier  series  for  a  square  wave  is  given  as 


f(x)  =  —  ^  sin(wA:)  [3.30] 

This  equation  is  plotted  in  Figure  3.18  with  200  (essentially  oo)  terms  included.  It  is 
plotted  again  in  Figure  3.19  with  the  tenth  term  missing.  The  result  is  the  sine  wave 
superimposed  on  the  square  wave.  The  absence  of  the  term  in  the  infinite  sum  makes  the 
presence  of  the  term  in  the  response. 
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Similarly,  with  the  closed-loop  zeros  repeated  infinitely  in  the  s*-plane  their  effect  is  to 
“eliminate”  the  frequency  at  that  point.  When  these  zeros  are  cancelled  by  the  open-loop 
poles,  the  effect  is  to  eliminate  all  the  harmonics  of  that  frequency,  except  that  one  which 
will  show  up  in  the  continuous  response  of  the  system. 

Examining  Figures  3.14  through  3.17  show  the  effect  changing  the  sample  period  has  on 
the  continuous  response.  It  is  also  easy  to  see  that  certain  frequency  ratios  give  better 

2 

responses  than  others.  Indeed,  frequency  ratios  corresponding  to  FR  =  —  for  n  =  even 

n 

give  better  response  than  those  in  which  n  =  odd  .  These  correspond  to  closed-loop  zeros 
located  on  the  positive  real  axis  in  the  z-plane,  with  the  corresponding  location  of  these 
zeros  on  the  mid-band  line  in  the  s*-plane. 

Figure  3.20  shows  a  comparison  of  the  impulse  and  step  responses  of  the  discrete-to- 
continuous  system  at  the  frequency  ratios  discussed  above.  Notice  that  by  slowing  the 
sample  rate  one  may  actually  improve  the  system  response,  a  fact  opposite  to  common 
practice  of  improving  response  by  sampling  faster. 

The  fact  that  the  continuous  response  between  samples  is  that  of  the  open-loop  plant  is 
more  easily  seen  as  the  sample  rate  is  further  decreased,  thus  allowing  more  aliasing. 
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Figure  3.21  shows  the  responses  of  the  discrete-to-continuous  system  when  the  frequency 


ratio  is  set  to  FR  = — ,  corresponding  to  open-loop  poles  located  in  the  quaternary  band. 

6.0 

Easily  discemable  in  the  impulse  and  step  responses  are  three  cycles  of  the  open-loop 
response  between  sample  periods.  Thus  as  stated  before,  the  continuous  response 
consists  of  a  series  of  step  responses,  given  the  initial  conditions  which  correspond  to  the 
open-loop’s  response  at  the  sample  period.  Since  the  controller  has  no  capability  of 
controlling  between  sample  periods,  this  free  response  of  the  system  will  always  be 
present  with  no  possibility  of  control  between  sample  periods. 

Figures  3.22  through  3.25  show  the  response  of  the  discrete-to-continuous  system  to  two 
different  frequency  sine  waves.  The  discrete  response  remains  essentially  unchanged 
from  the  “continuous”  (or  non-aliased)  case  of  Figure  3.5,  while  the  continuous  response 
seems  to  worsen  for  the  smaller  frequency  ratios.  One  can  see  the  increased  magnitude 
of  the  continuous  response  as  the  frequency  ratio  decreases  from  FR  =  2.0  to  Fi?  =  1.0 . 

The  sinusoidal  responses  show  the  drawback  of  allowing  aliasing.  The  controller  can 
only  control  the  system  at  the  sample  periods.  It  can  do  nothing  for  the  continuous 
response  between  samples.  The  continuous  output  between  sample  periods,  as  discussed 
before,  is  the  damped  natural  response  of  the  open-loop  plant.  Nothing  can  be  done  to 
control  this  response  with  aii  aliased  discrete-to-continuous  system,  and  the  responses  of 
the  system  to  sinusoidal  input  seen  in  Figures  3.22  through  3.25  again  reaffirm  this. 
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Figure  3.22  Sinusoidal  Responses  -  FR  =  2.0,  Kp  =  0.15 
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Response  to  sin(0.09*t)  Input;  FR  =  1.5;  Kp  =  0.15 


5  10  15  20  25  30  35  40 

Normalized  Time  (t/Ts) 


Figure  3.23  Sinusoidal  Responses  -  FR  =  1.5,  Kp  =  0.15 


Response  Response 


Response  to  stn(0.05*t)  input;  FR  =  1 ;  Kp  =  0.15 


Figure  3.25  Sinusoidal  Responses  -  FR  =  1 .0,  Kp  =  0. 15 
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3.5  Optimal  Proportional  Controller  Design  at  FR  =  1.0 


With  the  preliminary  examination  of  sample  period  effects  on  the  discrete-to-continuous 
system  with  proportional  control  complete,  it  is  obvious  that,  indeed,  the  designer  who 
must  allow  aliasing  can  judiciously  select  a  sample  rate  that  will  improve  performance. 

2 

Particularly,  the  sample  rates  such  that  FR  =  —  ,  n  =  even ,  give  superior  results  to  those 

n 

in  which  the  sample  rates  are  such  that  n  =  odd  .  In  both  cases  the  entire  root  loci  lie 
entirely  on  the  real  axis  in  the  z-plane.  However,  only  when  n  =  even  does  the  root  locus 
lie  partially  on  the  positive  real  axis.  Additionally,  in  the  s*-plane  for  these  cases  the  root 
loci  lie  on  either  the  band  boundaries  or  the  mid-band  line,  with  n  =  even  equating  to 
lying  on  the  mid-band  line  and  giving  the  better  response.  Also,  in  the  frequency  domain 
the  frequency  response  function  plots  show  that  when  n  =  even  the  sampling  to  damped 
natural  frequency  ratio  is  such  that  the  zeros  of  the  zero-order  hold  are  located  at  the 
same  frequencies  as  the  damped  natural  frequency,  another  reason  for  improved 
response. 

The  drawback  when  aliasing  is  allowed  is  the  controller  has  no  control  over  system 
response  between  sample  periods,  as  discussed  earlier.  The  more  aliasing  allowed  in  the 
system,  the  more  the  system  will  respond  between  sample  periods. 

Given  the  restraints  above,  it  appears  the  best  possible  sample  time,  giving  that  aliasing 
must  be  allowed,  is  that  which  will  make  the  frequency  ratio  equal  to  1 .0.  This  allows  for 
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the  benefits  from  allowing  n  =  even  and  avoids  the  pitfalls  of  open-loop  system  response 
for  too  many  system  periods.  With  this  in  mind,  a  further  exploration  of  the  discrete-to- 

2 

continuous  system  will  be  made  at  FR  = —  =  1.0. 

2.0 


As  discussed  earlier,  the  discrete  response  follows  what  is  expected  from  classical  theory 
(see  [  1 : 126]).  With  this  in  mind,  then,  a  pole  located  on  the  positive  real  axis  of  the 
z-plane  (but  within  the  unit  circle)  should  produce  an  exponentially  decaying  response, 
similar  to  an  overdamped  system.  This  was  indeed  the  case,  as  examination  of 
Figure  3.17  demonstrates.  A  pole  on  the  negative  real  axis  should  cause  exponentially 
decaying  oscillations,  similar  to  an  underdamped  system.  Again  this  was  the  case,  shown 
in  Figure  3. 14  for  a  FR  =  2.0.  Will  this,  however,  still  be  the  case  when  F/f  =  1 .0  if  the 
one  moving  pole  were  placed  there?  Theory  also  predicts  that  a  pole  placed  at  the  origin 
of  the  z-plane  will  produce  a  “deadbeat”  response  [36:598-601].  Again,  this  needs  to  be 
investigated  for  the  case  when  aliasing  is  allowed  and  FF  =  1.0 . 

Begin  by  placing  the  one  moving  pole  at  z  =  -1 .  This  should  correspond  to  a  non¬ 
decaying  oscillatory  response  of  the  discrete  system.  To  find  the  value  of  the  gain  that 
will  do  this,  use  the  denominator  of  Equation  3.19,  setting  the  poles  equal  to  z  =  -1  since 
that  is  where  one  pole  is  to  be  placed,  and  z  =  -b  since  one  sits  on  top  of  the  closed-loop 
zeros: 


(z  +  l)(z+b)  =  z^+d^z  +  d2  . 


[3.31] 
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Expanding  the  left-hand  side  of  Equation  3.31  and  substituting  the  values  of  and  , 
Equations  3.20  and  3.21,  into  the  right-hand  side  of  Equation  3.31  gives 


+(]b+\)z+b  =  z^  +  {K^c,  +  2b)z^-{K^c^  +  b'^)  .  [3.32] 


This  gives  two  equations  from  which  may  be  found: 


b  +  1  —  K  C-^  -J"  2b  ,  and 


[3.33] 


b  =  Kc^+b'^ 

P  2 


[3.34] 


Solving  Equation  3.33  for  K  : 


l-b 


[3.35] 


Solving  Equation  3.34  for  K  : 


b(l-b) 


[3.36] 
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These  equations  appear  to  give  different  results  until  one  solves  for  the  ratio  of  ^  from 

2 

Equations  3.15  and  3.16,  setting  b^=b2=b  since  FR  =  —  : 

n 


^2  _  (^0^^  -a^b-  ajb)  , 

2  ^ 
{a^-a^b-aj}) 


[3.37] 


Substituting  this  result  into  either  Equation  3.35  or  3.36  shows  the  equivalency  of  these 
two  equations.  Setting  the  value  of  as  given  in  Equation  3.35  or  3.36  produces  the 

results  of  Figures  3.26  and  321.  These  plots,  a  set  of  root  loci  and  the  impulse  and  step 
responses,  are  for  the  same  values  of  ©„  and  ^  as  previously. 

Examination  of  Figure  3.26  shows  that  the  single  non-stationaiy  root  is  located  on  the 
stability  boundary  in  both  the  z-plane  and  the  s*-plane.  In  both  planes,  this  fact  would 
predict  a  non-damped  oscillatory  behavior.  Looking  at  Figure  3.27  shows  this  to  be  the 
case,  both  the  discrete  and  continuous  portions  of  the  response  show  this  oscillatory 
response.  The  continuous  response  shows  the  aperiodic  nature,  however,  resulting  from 
the  combination  of  the  discrete  and  continuous  systems. 

Comparison  of  Figure  3.27  with  Figure  3.14  shows  that  the  response  can  be  driven  to  the 
extreme  edge  of  stability  where  non-damped  oscillations  occur.  Notice  the  oscillation  of 
the  control  effort,  shown  as  the  dashed-dot  line  of  Figure  3.27.  Furthermore,  comparing 
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Z-Plane  Root  Locus;  FR  =  1  Kp  =  3.2717 


Figure  3.26  Root  Loci  for  FR  =  1.0  With  Non-Stationary  Pole  at  z  =  -1 
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Figure  3.27  Impulse  and  Step  Responses  for  FR  =  1.0 
with  Non-Stationary  Pole  at  z  =  -l 
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magnitudes  of  response,  it  would  appear  the  case  in  which  FR  =  \.0  is  worse,  until  one 
considers  the  magnitude  increase  in  the  gain  between  the  two  cases. 


The  final  case  of  interest  when  FR  =  1.0  is  that  of  the  system  when  the  non-stationary 
pole  is  located  at  z  =  0.  Equation  3.31  is  rewritten  with  this  condition: 


iz  +  0)iz+b)  =  z^  +  diz+d2  . 


[3.38] 


Solving  this  as  before  gives  the  value  of  the  gain  sought  for: 


[3.39] 


Figures  3.28  and  3.29  show  the  results.  The  z-plane  root  locus  of  Figure  3.28  shows  the 
non-stationary  root  is  located  at  z  =  0 ,  as  desired.  This  root  cannot  be  seen  on  the 
s*-plane  root  locus  since  it,  and  all  the  band-repeated  roots,  are  located  at  negative 
infinity. 

Perhaps  most  interesting  are  the  results  shown  in  the  impulse  and  step  responses  of 
Figure  3.29.  Not  only  does  the  discrete  response  show  the  deadbeat  response  predicted 
by  root  location  in  the  z-plane,  but  the  continuous  response  is  a  deadbeat  response  as 
well.  The  impulse  shows  a  step  “up”  during  the  first  sample  period  with  a  step  “down” 
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Figure  3.28  Root  Loci  for  FR  =1.0  With  Non-Stationary  Pole  at  z  =  0 
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Response  ,  ,  Response 


Impulse  Response  at  Ws/Wd  =  1;  Kp  =  1 .1358 
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Figure  3.29  Impulse  and  Step  Responses  for  FR  =  1 .0  with  Non-Stationary  Pole  at  z  =  0 
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during  the  second  period.  This  down-step  is  exactly  equal,  but  opposite,  to  the  up-step. 
In  other  words,  the  continuous  response  of  the  second  sample  period  is  equal  and 
opposite  to  the  continuous  response  of  the  first  time  period,  causing  the  deadbeat 
response  seen.  The  control  effort  for  this  response  is  plotted  as  the  dash-dot  lines  of 
Figure  3.29.  For  the  impulse  response,  the  control  effort  of  the  second  sample  period  is 
exactly  equal  and  opposite  to  the  control  effort  of  the  first  sample  period.  For  the  step 
response,  the  control  effort  does  not  oscillate  in  any  way,  and  is  not  large  in  magnitude. 
Similarly,  here,  the  step  response  only  takes  one  sample  period  to  obtain  the  final  value 
of  the  step  with  no  oscillations  whatsoever  after  the  initial  time  period. 

An  interesting  note  made  by  Astrom  and  Wittenmark  [38:254-255]  is  one  of  the 
drawbacks  of  traditional  deadbeat  controllers.  They  mention  that  the  magnitude  of  the 
control  signal  increases  with  a  decreasing  sample  period,  thus  giving  the  deadbeat 
controller  a  bad  reputation.  An  examination  of  the  deadbeat  controller  here  shows  very 
little  control  effort  is  required.  By  allowing  aliasing,  a  pure  deadbeat  response  is 
obtained  much  easier.  With  aliasing,  the  response  of  the  system  is  not  forced  to  a  much 
greater  level  than  the  system’s  own  natural  response.  Here,  advantage  is  taken  of  that 
natural  response  and  used  for  simple  and  elegant  control.  This  is  contrary  to  control 
required  for  deadbeat  response  without  allowed  aliasing. 

Research  continues  in  the  area  of  deadbeat  controllers.  Some  recent  works  include  [39], 
[40],  [41],  [42],  and  [43]. 


110 


To  answer  concerns  as  to  the  robustness  of  this  system,  examine  parameters  that  may 
change,  affecting  the  response  of  the  discrete-to-continuous  system.  Begin  first  by  again 
examining  Equation  3.39.  Recall  b  is  defined  as 


[3.40] 


for  the  case  of  FR  =  1 .0 .  One  can  easily  see  from  Equations  3.4  and  3.6  that 


PlP2=^n 


[3.41] 


when  FF  =  1.0 .  Given  this  fact,  the  exponent  from  Equation  3.40  becomes 


rr  271 
PxTs=Pi  —  =Pi 
CO. 


27C 


27Cp, 


■  =  PiTs 


2np2 


V(FiF2)1-C^ 


[3.42] 


Since  the  terms  for  the  poles  are  interchangeable  in  Equation  3.42,  the  term  b  becomes 
independent  of  the  plant’s  natural  frequency  at  this  frequency  ratio.  Equation  3.15  gives 
the  definition  for  c, ,  and  using  Equations  3.8  through  3.10  in  Equation  3.15  leaves 


Cl 


1 

PxPi 


b  b 

- 1 - 

PxiP\-Pi)  PiiPi-Pi) 


[3.43] 
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which,  following  algebraic  manipulation  becomes 


Cl 


Pi-  Pi-bp2+bp^ 

(Pl-P2)0  +  b) 

PlPliPx-P2) 

.  ®/(Pl-P2)  . 

=  \+b  . 


[3.44] 


Substituting  Equations  3.44  and  3.40  into  Equation  3.39: 


[3.45] 


Equation  3.45  along  with  Equation  3.42  show  for  the  deadbeat  response  is  dependent 
only  on  the  damping  ratio  ^ ,  and  not  the  undamped  natural  frequency.  The  model  for 
this  deadbeat  response,  therefore,  is  dependent  on  the  undamped  natural  frequency  (or 

damped  natural  frequency  since  co^  =  )  and  the  damping  ratio  for  calculation 

of  the  needed  sampling  frequency,  and  the  damping  ratio  for  calculation  of  the 
proportional  gain  .  Robustness,  then,  will  be  examined  for  changes  in  these  two 
parameters. 

Figure  3.30  shows  plots  of  the  impulse  and  step  responses  to  changes  in  the  natural 
frequency.  If  the  natural  frequency  of  the  open-loop  plant  were  to  change,  or  were 
incorrect  from  the  modeling,  the  frequency  ratio  would  change  as  a  result,  since 
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Impulse  Response  at  Kp  =  1 .1358 


Step  Response  at  Kp  =  1 .1358 


Figure  3.30  Robustness  Comparison  for  Change  in  Natural  Frequency 
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FR  = 


[3.46] 


Figure  3.30  shows  comparisons  for  a  ±  5%  change  in  the  natural  frequency. 


Finally,  Figure  3.31  shows  comparisons  for  a  ±5%  change  in  the  damping  ratio,  which 
again  affects  the  FR  according  to  Equation  3.46,  as  well  as  the  correct  proportional  gain 
for  the  deadbeat  response  according  to  Equation  3.45. 

Figures  3.30  and  3.31  demonstrate  the  robustness  of  this  system.  The  most  extreme 
changes  occur  with  changes  in  the  natural  frequency  of  the  open-loop  plant.  However  for 
this  case,  even  with  a  change  of  ±  5% ,  a  large  variation  in  a  parameter,  the  system  step 
response  settles  to  a  final  value  within  six  sample  periods.  The  change  in  responses  of 
the  system  with  changes  in  the  natural  frequency  is  greater  than  for  changes  in  (or 
damping). 

These  changes  of  the  system  response  from  the  optimal  value  demonstrate  one  of  the 
classical  problems  with  deadbeat  systems — ^the  cancellation  of  a  desired  frequency 
requires  accurate  knowledge  of  the  frequency.  Without  this  accuracy  it  becomes  difficult 
to  cancel  the  mode.  Given  the  accuracy  of  the  system  model,  however,  this  type  of 
deadbeat  response  is  difficult  to  match  with  conventional  methods. 
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Response 


Impulse  Response  at  Kp  =  1 .1358 


Step  Response  at  Kp  =  1 .1358 


Figure  3.31  Robustness  Comparison  for  Change  in  Damping  Ratio 
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As  is  always  the  case  with  a  control  system  consisting  of  proportional  control  only,  the 
steady-state  error  is  perhaps  unacceptable.  Traditionally,  this  error  is  taken  care  of  by 
introducing  an  integrator.  The  next  chapter  will  examine  such  a  system  to  try  and 
achieve  comparable  results  as  the  proportional  controller,  but  without  the  steady-state 
error  seen  here. 


3.6  Chapter  Three  Summary 


In  chapter  Three  the  root  locus  for  the  discrete-to-continuous  system  was  developed  and 
the  movement  of  the  poles  in  the  gain  limits  for  the  system  were  examined.  They  were 
found  to  start  at  the  locations  of  the  closed-loop  discrete  zeros,  which  correspond  to  the 
open-loop  continuous  poles,  and  move  toward  the  locations  of  the  open-loop  zeros  of 
H*(s).  A  second-order  system  with  proportional  control  only  was  developed  and 

examined  at  different  frequency  ratios  allowing  aliasing.  The  frequency  ratio  of  FR  =  — 

n 

was  found  to  cause  the  entire  z-plane  root  locus  to  lie  on  the  real  axis,  as  well  as  fixing 
one  of  the  z-plane  plant  poles.  When  n  is  odd,  part  of  the  root  locus  lies  on  the  positive 
real  axis  and  one  can  place  the  non-stationary  pole  at  the  origin.  This  produces  a 
deadbeat  step  response  with  a  simple  control  effort  contrary  to  tradition  deadbeat 
controllers.  It  was  seen  that  the  system  with  aliasing  will  not  follow  a  sinusoidal  input 
very  well,  giving  a  limitation  of  this  model — the  frequency  of  the  input  must  be  lower 
than  the  frequency  of  the  sampler. 
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Chapter  Four.  Second-Order  System  Design  with 

Aliasing  -  Proportional-Integral 
Control 


4.0  Introduction 

As  introduced  in  chapter  three,  the  sample  period  turns  out  to  be  a  powerful  tool  in  the 
design  of  a  discrete-to-continuous  control  system  with  a  proportional  controller  in  which 
aliasing  is  present.  The  deadbeat  response  seen  in  Figure  3.29  is  an  attractive  design  goal 
for  certain  applications,  given  that  the  steady-state  error  could  be  reduced.  Here,  the 
search  for  such  a  system  will  begin,  starting  with  an  examination  of  the  integral  controller 
and  moving  on  to  the  proportional-integral  controller.  Comparison  between  root  loci, 
frequency  response  functions,  impulse  and  step  responses,  and  responses  to  sinusoidal 
inputs  will  again  be  examined.  Again,  the  frequency  ratio  allowing  for  the  best  response 
with  aliasing  present  is  found  to  be  F/?  =  1.0 .  A  deadbeat  controller  with  zero  steady- 
state  error  is  designed,  and  robustness  and  disturbance  rejection  for  this  controller  is 
examined. 
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4.1  Design  of  an  Integral  Controller  Using  Sample  Period  As  a 
Design  parameter 


Similar  to  the  investigation  of  the  proportional  controller  in  chapter  three,  here  examine 
the  behavior  of  the  integral  controller  as  sample  period  changes.  The  integral  controller 
is  defined  by 


GAz)  = 


K,z 

z-1 


[4.1] 


Solving  for  F{z)  from  Equation  3.3  gives 


F(z)  = 


_ Kizjz+biXz+b^) _ 

Kiz(CiZ  +  C2)  +  (z-l){z+bi)(z+b2) 


[4.2] 


where  the  constants  are  as  defined  previously  in  Equations  3.13  through  3.16.  This 
equation  becomes 


z^+d,z^  +  d^z  +  d,  ’ 


[4.3] 


where 


d,  —  +  ^>2  ~  1  j 


[4.4] 
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When  Ki  is  equal  to  zero  the  poles  of  Equation  4.3  are  located  at  z  =  1 ,  -b^,  and  -  , 

which  correspond  to  the  pole  of  the  controller  and  the  poles  of  the  open-loop  plant,  were 

they  transformed  to  the  z-plane.  At  Kj  =  «>  the  poles  move  to  z  =  0 ,  -  ,  and  — ^ ,  the 
location  of  the  controller  zero  and  the  zeros  of  H{z) .  This  is  as  expected. 

For  comparison  purposes,  the  discrete-to-continuous  system  consisting  of  the  second- 
order  plant  of  chapter  three  (Equation  3.6)  and  the  controller  of  Equation  4.1  with  a 
frequency  ratio  of  FR  =50  is  used  for  Figures  4.1  through  4.4.  These  correspond  to 
plots  in  which  no  aliasing  is  present — ^i.e.,  a  system  for  which  traditional  design  methods 
would  be  used.  A  very  low  gain  of  =  0.01  is  used  for  stability.  Notice  the  very  small 
gain  margin  of  this  system,  a  direct  result  of  using  an  integral  controller  by  itself. 

Now  allow  aliasing  as  before  in  chapter  three.  The  frequency  ratio,  for  non- 
dimensionalization,  will  again  be  as  defined  in  Equation  3.22.  Figures  will  be  presented 

2  2  2  2 

for  frequency  ratios  of  FR  = — ,  — ,  — ,  and  — . 

1.0  1|  If  2.0 


Figure  4.1  Root  Loci  for  FR  =  50,  Ki  =  0.01 
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Frequency  Ratio  Ws/Wd  =  50;  Kl  =  0.01 


Figure  4.2  Frequency  Response  Function  Plot  for  FR  =  50,  Ki  =  0.01 
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Response  ^  ^  Response 


Figure  4.3  Step  and  Impulse  Responses  for  FR  =  50,  Ki  =  0.01 


X 1 0"^  Impulse  Response  at  WsA/Vd  =  50;  KI  =  0.01 


Step  Response  at  Ws/Wd  =  50;  Ki  =  0.01 
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Response  Response 


Figure  4.4  Response  to  Sinusoidal  Inputs  for  FR  =  50,  Ki  =  0.01 
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Figures  4.5  through  4.8  show  the  root  loci  plots  for  the  discrete-to-continuous  second- 
order  system  with  an  integral  controller.  Recall  the  definition  of  symbols  on  the  plots: 
triangle  and  plus  -  zero-order  hold  zeros  and  pole,  square  and  asterisks  -  plant  zeros  and 
poles,  diamond  -  pole  of  H{z) ,  and  circles,  “x’s”  and  solid  line  -  zeros  and  poles  of 
F(z) .  Notice  the  increase  in  allowable  gain  before  the  stability  boundaries  are  crossed. 
Again,  by  slowing  the  sample  rate,  stability  is  introduced,  contrary  to  current  design 
practices. 

Figure  4.5  shows  the  entire  root  locus  for  a  frequency  ratio  of  FR  =  2.0  lies  entirely  on 
the  real  axis  of  the  z-plane.  Similar  to  that  of  the  proportional  controller,  it  is  different, 
however,  in  that  a  portion  of  it  lies  on  the  negative  real  axis,  and  another  portion  lies  on 
the  positive  real  axis.  The  portion  that  lies  on  the  positive  real  axis  consists  of  a  pole  that 
begins  at  z  =  1 ,  the  open-loop  controller  pole  location,  and  moves  toward  z=0,  the 
open-loop  controller  zero  location.  This  portion  of  the  root  locus  lies  entirely  within  the 
stability  range.  The  portion  that  lies  on  the  negative  real  axis  consists  of  one  pole  fixed  at 
the  location  of  the  zeros  of  F(z)  found  from  the  denominator  of  H(z)  and  one  pole  that 
travels  toward  negative  infinity.  Once  this  pole  passes  z  =  —  1 ,  the  system  goes  unstable. 

The  portion  of  the  z-plane  root  locus  that  lies  on  the  positive  axis  corresponds  to  the 

portion  of  the  s*-plane  root  locus  that  lies  on  the  mid-band  line.  These  poles  start  at 

5*  =  0  and  move  toward  negative  infinity,  the  point  corresponding  to  z  =  0 .  The  other 

non-stationary  poles  starts  at  the  values  of  s*  corresponding  to  the  zeros  of  F  *  (j)  and 

move  along  the  band  boundaries  to  infinity. 
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Figure  4.5  Root  Loci  for  FR  =  2.0,  Ki  =  0.3 
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Figures  4.6  and  4.7  show  the  root  loci  for  frequency  ratios  of  FR  =  —  and  — , 
respectively.  Here,  there  is  no  stationary  pole,  and  the  same  general  trends  are  observed. 


Figure  4.8  shows  the  root  loci  for  the  special  frequency  ratio  of  F/?  =  1.0 .  Note  the 
unique  nature  of  these  root  loci,  even  though  there  are  similarities  of  characteristics  with 
the  like  root  loci  of  the  system  with  proportional  control  of  Figure  3.9.  At  this  frequency 
ratio  one  pole  remains  fixed  at  the  location  of  the  zeros  of  F(z)  found  from  the 
denominator  term  of  H (z)  (see  Equation  3.3).  The  second  pole  starts  at  these  zeros  and 
travels  toward  the  pole  that  starts  at  z  =  1 ,  the  controller  pole  location,  and  travels  toward 
the  second  pole.  Once  these  poles  meet,  they  travel  in  a  circle  of  constant  radius  around 
the  zero  found  at  the  origin,  meet  at  the  real  axis,  and  then  split,  with  one  pole  traveling 
toward  negative  infinity  and  the  other  traveling  toward  z  =  0 ,  the  location  of  the 
controller  zero. 


This  constant  radius  circle  in  the  z-plane  corresponds  to  a  constant  value  of  the  real  part 
of  the  locus  in  the  s*-plane.  This  can  be  seen  as  the  change  of  only  the  imaginary  portion 
of  the  root  on  the  locus  in  the  s*-plane  of  Figure  4.8.  These  repeated  roots  in  the  s*-plane 
travel  to  the  band  boundaries,  where  they  “meet”  the  root  coming  from  the  adjacent  band, 
then  separate  with  one  traveling  to  negative  infinity  and  the  other  traveling  to  the  location 
of  the  zero  of  ff(z)  transformed  to  the  s*-plane  and  repeated  infinitely. 
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Z-Plane  Root  Locus;  FR  =  1 .5  Ki  =  0.3 


S*-Plane  Root  Locus;  FR  =  1 .5  Ki  =  0.3 


Figure  4.6  Root  Loci  for  FR  =  1 .5,  Kj  =  0.3 
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Z-Plane  Root  Locus;  FR  =  1 .2  Ki  =  0.3 


S*--Plane  Root  Locus;  FR  =  1 .2  Ki  =  0.3 


Figure  4.7  Root  Loci  for  FR  =  1.2,  Ki  =  0.3 
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Figure  4.8  Root  Loci  for  FR  =  1.0,  Ki  =  0.3 
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Figures  4.9  through  4.12  show  the  frequency  response  function  magnitude  and  phase 
plots  of  the  second-order  discrete-to-continuous  system  with  integral  control.  Similar 
trends  to  those  discussed  in  chapter  three  for  the  proportional  controller  can  be  observed 
here. 

Figures  4.13  through  4.16  show  the  impulse  and  step  responses  of  the  system  at  the 
frequency  ratios  defined  above.  Again,  the  discrete  responses  follow  the  response 
predicted  by  the  location  of  the  poles  in  the  z-plane.  The  response  at  a  frequency  ratio  of 
F/?  =  1 .0  is  much  better  behaved  than  that  at  a  frequency  ratio  of  FR  =  2.0 .  The  step 
response  at  FF  =  1 .0 ,  however,  does  not  appear  to  be  the  best  of  the  four  frequency 
ratios— the  step  response  at  FF  =  1.5  seems  to  be  best  overall.  A  comparison  of  these 
responses  is  plotted  in  Figure  4.17. 

Figures  4.18  through  4.21  show  the  responses  of  the  second-order  discrete-to-continuous 
system  to  sinusoidal  inputs.  Unlike  those  for  the  proportional  controller,  differences  can 
be  seen  when  integral  control  is  employed.  Here,  there  is  a  definite  degradation  of 
performance  as  the  sample  rate  decreases. 

Looking  now  at  the  case  where  FF  =  1.0 ,  as  with  the  proportional  controller,  an  optimal 
gain  will  be  calculated.  It  would  seem  this  gain  would  be  such  that  the  non-stationary 
poles  are  at  the  same  location  in  the  z-plane,  meeting  on  the  positive  z-axis.  Calling  this 
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Magnitude  (dB) 


Frequency  Ratio  Ws/Wd  =  1 ,5;  Kl  =  0.3 


Figure  4.10  Frequency  Response  Function  Plot  -  FR  =  1 .5,  Kj  =  0.3 
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Frequency  Ratio  Ws/Wd  =  1 .2;  Ki  =  0.3 


Figure  4.11  Frequency  Response  Function  Plot  -  FR  =  1 .2,  Ki  =  0.3 
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Magnitude  (dB) 


Response  ,  ,  Response 


Impulse  Response  at  Ws/Wd  =  2;  Ki  =  0.3 


Step  Response  at  Ws/Wd  =  2;  Ki  =  0.3 


Figure  4.13  Impulse  and  Step  Responses  -  FR  =  2.0,  Kj  =  0.3 
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Response  Response 


Impulse  Response  at  Ws/Wd  =  1 .5;  Ki  =  0.3 


Step  Response  at  Ws/Wd  =  1 .5;  Ki  =  0.3 


Figure  4.14  Impulse  and  Step  Responses  -  FR  =  1.5,  Ki  =  0.3 
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Response  ^  Response 


Impulse  Response  at  WsA/Vd  =  1.2;  Ki  =  0.3 


Step  Response  at  Ws/Wd  =  1 .2;  Ki  =  0.3 


Figure  4.15  Impulse  and  Step  Responses  -  FR  =  1.2,  Ki  =  0.3 
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Response  Response 


Impulse  Response  at  Ws/Wd  =  1 ;  Ki  =  0.3 
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Responi 


Comparative  Impulse  Response  at  Ki  =  0.3 


Comparative  Step  Response  at  WsA/Vd  =  Ki  =  0.3 


Normalized  Time  (t/Ts) 


Figure  4.17  Comparison  of  Impulse  and  Step  Responses  for  Varying  Sample  Rate 
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Response  Response 


Response  Response 


Response  to  sin(0.05*t)  input;  FR  =  1 .5;  Ki  =  0.3 


Response  to  sln(0.09*t)  input;  FR  =  1 .5;  KI  =  0.3 


Figure  4.19  Sinusoidal  Responses  -  FR  =  1.5,  Kj  =  0.3 
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Response  to  sin(0.05*t)  input;  FR  =  1.2;  Ki  =  0.3 


Figure  4.20  Sinusoidal  Responses  -  FR  =  1 .2,  Ki  = 


Response  Response 


Response  to  sjn(0.05*t)  input;  FR  =  1 ;  Ki  =  0.3 


Response  to  sin(0.09*t)  input;  FR  =  1 ;  KI  =  0.3 


Figure  4.21  Sinusoidal  Responses  -  FR  =  1 .0,  Kj  =  0.3 
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location  on  the  z-axis  ,  and  recalling  that  b^=b2=b,  the  denominator  of  Equation  4.3 
may  be  set  equal  to  the  desired  location  of  the  poles: 


iz-bjf(z+b)  =  z^ +d^z^ +d2Z  +  d^  . 


[4.7] 


Substituting  in  the  values  of  Equations  4.4  through  4.6  and  performing  the  multiplication 
gives 


+  {b-2b^)z^  +  (b,^ - 2bb^)z  +bb^^  =  z^  +  +2b- l)z^  +  (K.c^  +b^- 2b)z - b^  . 


[4.8] 


Equating  same-order  terms  of  Equation  4.8  gives  three  equations  with  the  two  unknowns, 
^3  and  Kr. 


b  —  2^3  —  "I"  —  1 


[4.9] 


b^  —  2bbj  =  K1C2  +b^  —2b  ,  and 


[4.10] 


[4.11] 
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Solving  Equation  4.1 1  for  gives 


b,=4^  . 


[4.12] 


Substituting  Equation  4.12  into  Equation  4.9  and  solving  gives 


\-b-2^ 

Cl 


[4.13] 


Substituting  Equations  4.12  and  4.13  into  Equation  4.10,  and  recalling  Equation  3.37 
(c2  =bc^)  gives  independent  proof  of  correct  solution  for  the  two  unknowns. 

Figures  4.22  and  4.23  are  the  plots  of  the  root  loci  and  impulse  and  step  responses  when 
is  as  given  in  Equation  4.13.  Examination  of  Figure  4.22  shows  placement  of  the 

poles  is  as  desired.  Figure  4.23  shows  a  much  improved  impulse  response,  the  best  of  the 
previous  plots,  and  a  step  response  with  no  discrete  overshoot.  Again,  this  step  response 
is  not  the  “best”  when  compared  to  Figures  4.14  and  4.15.  No  results  such  as  deadbeat 
response  for  the  proportional  controller,  then,  are  seen  here,  though  this  was  not 
expected. 
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Response  „  Response 


Figure  4.23  Impulse  and  Step  Responses  -  FR  =  1.0,  Ki  =  0.15657 
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4.2  Design  of  a  Proportional-Integral  Controller  at  FR  =  1.0 


Having  a  firmer  understanding  of  the  integral  controller  in  which  aliasing  is  allowed  now 
allows  the  design  of  a  proportional-integral  (PI)  controller  in  which  aliasing  is  allowed. 
The  plant  used  in  this  design  will  again  be  the  second-order  plant  of  Equation  3.4. 

The  controller,  for  PI  control  is  defined  as 


z-l 


{K^  +  K,)z-K^ 
z-\ 


[4.14] 


Using  this  in  Equation  3.3  and  solving  for  F{z)  gives 


F{z)  = 


[{K^  +  K,)z-K^]<,z+b,)iz+b^) 

[{K^  +  K,)z - K^]{c,z  +  cj)  +  (z - l)(z  +b,){z^b^)  ’ 


[4.15] 


and  again  the  constants  b^,  b2,  c^,  and  Cj  are  defined  in  Equations  3.13  through  3.16. 
Multiplying  the  denominator  gives 


^  [(Kp  +  K.)z  -  Kp ](z  +  by){z  +^2) 
z^ +d^z^ -^d^z  +  d^ 


[4.16] 
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where 


d,={K^  +  K,)c^+b^+b^-\  ,  [4.17] 

d^=iK^  +  K,)c^-K^c^+bp^-b,-b^  ,and  [4.18] 

d^=-KpC^-bJb^  .  [4.19] 

Since  the  deadbeat  response  is  the  aim  here,  select  the  frequency  ratio  of  FR  =  1 .0  and 
the  gain  of  the  proportional  portion  of  the  controller  to  be  equal  to  that  which  gives  the 
deadbeat  response: 


K,=--  .  [4.20] 

Cl 

To  begin,  suppose  the  gain  for  the  integral  portion  of  the  controller  is  as  found  in 
Equation  4.13.  Figures  4.24  and  4.25  show  the  root  loci  and  impulse  and  step  responses 
for  these  values  of  the  gain.  Examination  of  Figure  4.24  shows  the  two  “plant”  poles 
located  on  top  of  the  zeros  as  expected  and  the  integrator  pole  close  to  the  point  z  =  1 . 
This  seems  to  be  causing  the  slow  integrator  response  seen  in  the  step  response  of  Figure 
4.25.  Here  again,  the  control  is  shown  as  the  dash-dot  line.  The  impulse  response 
appears  to  be  deadbeat,  but  closer  examination  shows  small  oscillations  for  many  sample 
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-0.5  -0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5 


Figure  4.24  Root  Loci  -  FR  =  1.0,  Kp  =  1.1358,  Ki  =  0.15657 
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Response  ,  ,  Response 


Impulse  Response  at  Ws/Wd  =  1;  Kp  =  1.1358;  Ki  =  0.15657 


Step  Response  at  Ws/Wd  =  1 ;  Kp  =  1 . 1 358;  KI  =  0. 1 5657 


Figure  4.25  Impulse  and  Step  Responses  -  FR  =  1.0,  Kp  =  1.1358,  Ki  =  0.15657 
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periods.  The  step  response  shows  the  characteristics  of  the  deadbeat  response,  but  with 
the  addition  of  the  integrator,  this  response  is  no  longer  deadbeat. 

To  try  and  increase  the  integrator  response,  place  the  pole  for  the  integrator  closer  to 
z  =  0 ,  but  with  it  remaining  on  the  positive  axis.  To  place  it  directly  at  z  =  0  the 
denominator  of  F(z)  must  be 

z^{z+b)  =  z^ +  d^z^ +d^z  +  dj^  .  [4.21] 

Using  Equations  4.17  through  4.19  in  Equation  4.21  and  solving  for  the  values  of  and 
gives 


[4.22] 


which  is  the  same  as  it  was  for  the  deadbeat  response  (Equation  3.39)  and 


[4.23] 


Figures  4.26  and  4.27  show  the  results  of  using  these  gains.  Two  roots  are  indeed  placed 
at  z  =  0  corresponding  to  =  -oo ,  with  the  third  root  placed  at  the  location  of  the  zeros 
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Z-Plane  Root  Locus;  FR  =  1 ;  Kp  =  1 .1 358;  Ki  =  2.1 358 


S*-Plane  Root  Locus;  FR  =  1 ;  Kp  =  1.1358;  Ki  =  2.1358 


Figure  4.26  Root  Loci  -  FR  =  1.0,  Kp  =  1.1358,  Ki  =  2.1538 
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Response  ,  Response 


Impulse  Response  at  Ws/Wd  =  1;  Kp  =  1 .1358;  Ki  =  2.1358 


5. 


found  from  the  denominator  term  of  H{z)  ■  As  a  result,  the  response  is  closer  to  the 
sought-after  deadbeat  response.  The  impulse  response  returns  to  zero  following  the  third 
sample  period,  and  the  step  response  is  at  a  value  of  one  after  two  sample  periods. 
Overshoot  increases  by  quite  a  lot  in  this  case.  Control  effort  is  also  plotted. 

Perhaps  placing  the  pole  on  the  location  of  the  other  stationary  pole  will  give  better 
response.  This  zero  is  located  at 

z  =  -b  .  [4.24] 

Equation  4.21  now  becomes 

z{.z+bf  =  z^ +d{z^ +  d2Z  +  d.^  ,  [4.25] 


and  the  solution  process  yields 


[4.26] 


^  ^(Cl+Cz)^ 

CjC2 


[4.27] 
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Equation  4.27  simplifies  with  the  following.  Recall  Equation  3.15,  repeated  here  for 
clarity: 


Cl  =  OJ„^(ao  - .  [4.28] 

But,  at  F/?  =  1 .0 ,  b^-b2=b  and  Equation  4.28  becomes 

Cl  =  -  b{a^  +a^)]  .  [4.29] 

Since  a^+a^^-a^  (see  Equations  3.8  through  3.10)  Equation  4.29  becomes 

Cl  =  ®/(ao  +  K)  =  + b)  .  [4.30] 

From  Equation  3.8,  ^ ,  so 

PiPi 

Ci=l+^  •  [4.31] 

Likewise 

c^  =  cp  =  b{l+b)  .  [4.32] 

Substituting  Equation  4.31  and  4.31  into  Equation  4.27  gives 
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[4.33] 


^  _  (C;  +C2)b  _  \X+b+b{l  +  b)]b  ^ 

C1C2  ”  a  +  b)b(l  +  b) 

Using  the  values  of  from  Equation  4.26  and  Kj  =1.0,  Figure  4.28  shows  the 

placement  of  the  poles  and  controller  zero.  Notice  the  location  of  the  controller  zero 
exactly  opposite  the  zeros  found  from  the  denominator  of  ff(z)  in  the  z-plane.  In  the 
s*-plane,  these  zeros  all  lie  at  the  same  real  value  in  the  Jco  direction.  Two  poles  are 
located  at  these  zeros  with  the  third  pole  at  z  =  0  or  s*  =  -oo. 

Figure  4.29  shows  the  desired  deadbeat  response  in  both  the  impulse  and  step  responses. 
Control  effort  is  also  plotted  with  the  dash-dot  line.  There  is  a  large  amount  of  overshoot 
in  both  cases.  However,  the  response  is  truly  deadbeat,  and  with  the  addition  of  the 
integrator,  the  final  desired  value  of  both  the  impulse  and  step  is  attained  with  no 
oscillations.  Additionally,  control  effort  to  attain  this  response  is  simple  and  not  overly 
large.  This  method  of  design,  then,  produces  some  exciting  possibilities  for  control. 

In  examining  the  impulse  response,  one  can  see  that  the  sample  rate  is  such  that  the 
response  resulting  from  the  unit  impulse  at  time  t  =  0,  actually  a  step  function  to  the 
continuous  system  once  the  zero-order  hold  is  applied,  is  exactly  equal  and  opposite  to 
the  impulse  resulting  from  the  impulse  returning  to  zero  at  time  f  =  T^.  The  correct 

selection  of  the  sample  period  allows  this  “phasing”  of  the  response  and  its  affect  on  the 
control  effort.  This  method  of  control  is  reminiscent  of  the  idea  behind  input  shaping. 
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-0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4 


Figure  4.28  Root  Loci  -  FR  =  1.0,  Kp  =  1.1358,  Ki  =  1 .0 
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Response  Response 


Impulse  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .1358;  Ki  =  1 


Normalized  Time  (t/Ts) 


Step  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .1 358;  Ki  =  1 


0  2  4  6  8  10 


Normalized  Time  (t/Ts) 


Figure  4.29  Impulse  and  Step  Responses  -  FR  =  1.0,  Kp  =  1.1358,  Kj  =  1.0 
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Currently  much  work  is  being  accomplished  on  a  method  of  control  called  input  or 
command  shaping,  including  the  recent  examples  cited  here  [44],  [45],  [46],  [47],  and 
[48].  Input  shaping  is  a  feedforward  technique  used  to  eliminate  vibration  from  the 
system  under  control.  It  consists  of  convolving  a  sequence  of  impulses  with  the  desired 
command  to  the  system,  with  the  resulting  shaped  input  used  to  drive  the  system.  This 
input  effectively  cancels  the  vibrating  poles  regardless  of  the  desired  system  input.  One 
can  see  the  similarities  between  input  shaping  and  the  results  of  the  simple  Pl-controller 
of  this  dissertation. 

Examining  the  response  of  the  deadbeat  system  found  above  to  a  sinusoidal  input  shows 
again  the  limitation  of  such  a  system.  Figure  4.30  shows  the  response  to  a  sinusoid  with 
frequencies  lower  than  the  sampling  frequency.  The  controller  is  able  to  follow  the  input 
adequately  at  the  discrete  points,  with  the  overshoot  seen  in  the  deadbeat  controller 
coming  between  each  sample  period.  Figure  4.3 1  shows  the  response  of  the  system 
when  the  inputs  frequencies  are  greater  than  the  sampling  frequency.  Here,  the  controller 
cannot  follow  the  input  signal.  Further  showing  of  the  limitation  of  this  controller— input 
frequencies  cannot  be  greater  than  that  of  the  sampler. 

Robustness  of  this  system  is  again  examined  for  five  percent  changes  in  the  system 
parameters.  Figure  4.32  plots  the  impulse  and  step  responses  as  the  natural  frequency  is 
varied,  and  Figure  4.33  plots  the  responses  for  changes  in  the  damping  ratio.  Results  of 
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Response  Response 


Response  to  sin(0.05*t)  Input;  FR  =  1 ;  Kp  =  1.1358;  Ki  =  1 


Response  to  sin{0.09*t)  input;  FR  =  1 ;  Kp  =  1 .1358;  Ki  =  1 


Figure  4.30  Response  to  Sinusoidal  Inputs  at  Low  Frequencies 
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Response  Response 


Response  to  sin(0.5*t)  input;  FR  =  1 ;  Kp  =  1 .1358;  Ki  =  1 


Response  to  sin(0,9*t)  input;  FR  =  1;  Kp  =  1 .1358;  Ki  =:  1 


Figure  4.31  Response  to  Sinusoidal  Inputs  at  High  Frequencies 
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Comparison  of  Impulse  Responses  at  FR  =  1.0 


Comparison  of  Impulse  Responses  at  FR  =  1 .0 


Figure  4.32  Robustness  Comparisons  for  Change  in  Natural  Frequency,  Kj  =  1.0 
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Response 


Comparison  of  Impulse  Responses  at  FR  =  1.0 


Comparison  of  Impulse  Responses  at  FR  =  1 .0 


Figure  4.33  Robustness  Comparisons  for  Change  in  Damping,  Ki  =  1.0 


6. 


these  plots  correspond  closely  to  those  of  the  system  controlled  by  proportional  gain  only. 
Since  has  a  constant  value  for  the  optimal  response,  no  plots  for  changes  in  jRT,-  are 

shown.  Here,  again,  the  changes  caused  by  a  change  in  natural  frequency  of  the  system 
are  greatest.  Recall  the  dual  dependence  on  natural  frequency — correctly  identifying  the 
sampling  frequency  for  a  frequency  ratio  of  FR=\ .0  and  correctly  calculating  the 
proportional  gain  . 


4.3  Disturbance  Rejection  of  Proportional-Integral 
Controller 

One  concern  about  a  controller  design  is  its  ability  to  reject  disturbances.  An 
examination  of  this  controller’s  disturbance  rejection  capability  will  now  follow.  The 
disturbances  to  be  introduced  are  a  step  disturbance  and  a  normally  distributed  random 
disturbance. 

Figure  4.34  shows  the  impulse  and  step  responses  of  proportional-integral  (PI)  with  a  unit 
step  disturbance  at  a  normalized  time  of  —  =  5.0 ,  or  following  five  sampling  periods. 

'^s 

Here  the  disturbance  is  introduced  at  the  sampling  instant.  Hence,  the  controller,  shown 
as  the  dashed  line,  does  not  yet  recognize  the  disturbance.  One  can  see,  however,  that  as 
soon  as  the  disturbance  is  recognized,  the  controller  rejects  the  disturbance,  and  does  so 
within  one  sample  period. 
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Response  Response 


Impulse  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .1 358;  Ki  =  1 


Step  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .1358;  Ki  =  1 


Figure  4.34  Unit  Step  Disturbance  Rejection  —  Disturbance  at  Sample  Point 
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Figure  4.35  shows  the  result  when  the  disturbance  is  between  sample  periods — here  65% 
of  the  discrete  sample  period  has  passed.  Again  the  controller  rejects  the  disturbance 
once  it  “realizes”  it  is  there,  and  does  so  within  one  sample  period. 

Figure  4.36  shows  the  same  disturbance  rejection  for  a  greater  magnitude  of  the  step 
disturbance.  Notice  that  overshoot  of  the  continuous  response  increases,  but  this  is  to  be 
expected. 

Finally,  Figure  4.37  shows  the  disturbance  rejection  when  the  step  is  in  the  negative 
direction.  Examination  of  Figures  4.34  through  4.37  shows  that  this  system  has  excellent 
step  disturbance  rejection  capabilities.  This  can  be  expected  since  the  step  is  basically  a 
zero-frequency  disturbance,  and  the  controller  handles  input  at  frequencies  lower  than  the 
sampling  frequency  with  ease. 

The  controller’s  ability  to  reject  a  broadband  disturbance  will  now  be  examined.  The 
disturbance  to  be  introduced  next  will  be  a  normalized  random  input,  with  vector  values 
of  the  disturbance  vector  generated  by  Matlab’s  “RANDN”  command  [49].  The  resulting 
pseudo-random  numbers  are  normally  distributed  with  zero  mean  and  a  variance  equal  to 
one.  The  random  number  generator  is  reset  for  each  run,  meaning  the  random  input  for 
each  of  the  following  figures  is  the  same. 
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Response  Response 


Impulse  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .1358;  Ki  =  1 


Step  Response  at  Ws/Wd  =  1;  Kp  =  1 .1358;  KI  =  1 


Figure  4.35  Unit  Step  Disturbance  Rejection 
Disturbance  65%  of  Sample  Period  Past  Sample  Point 
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Response  Response 


Impulse  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .1358;  Ki  =  1 


Step  Response  at  WsA/Vd  =  1 ;  Kp  =  1 .1358;  KI  =  1 


Normalized  Time  (t/Ts) 


Figure  4.36  Step  Disturbance  Rejection  -  Step  Size  =  5.0 
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Figure  4.38  shows  the  system  with  the  random  noise  input  starting  just  after  the  fifth 
sample  period  and  continuing  until  just  before  the  seventh  sample  period.  Again,  the 
control  effort  is  shown  as  the  dashed  line.  During  the  random  disturbance,  the  controller 
can  only  control  for  the  point  it  “saw”  at  the  sample  point.  Thus,  little  if  any  control  is 
accomplished.  Following  the  disturbance,  however,  within  one  sample  period  the  system 
returns  to  the  correct  steady-state  values. 

Figure  4.39  shows  similar  results  for  the  disturbance  beginning  between  sample  periods. 
Here,  one  can  see  that  the  location  of  the  start  of  the  random  disturbance  affects  the 
response. 

The  ability  of  this  system  to  control  a  disturbance  with  a  broadband  frequency  range  is 
poor,  at  best.  This  again  follows  what  was  found  from  the  system’s  ability  to  follow  a 
sinusoidal  input  -  good  for  frequencies  below  the  sampling  frequency,  but  poor  for 
frequencies  on  the  order  of  or  greater  than  the  sampling  frequency. 


4.4  Chapter  Four  Summary 


In  chapter  four  we  examined  the  design  of  a  proportional-integral  (PI)  controller  for  a 
system  in  which  aliasing  is  allowed.  The  chapter  began  with  an  examination  of  the 
integral  controller  alone.  Similar  trends  were  seen  as  in  chapter  three.  With  the  PI- 
controller,  the  deadbeat  response  was  obtained  with  zero  steady-state  error  and 
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Response  ,  Response 


Impulse  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .1358;  Ki  =  1 


Step  Response  at  Ws/Wd  =  1 ;  Kp  =  1 ,1358;  Ki  =  1 


Figure  4.39  Random  Disturbance  Rejection 
Disturbance  Beginning  Between  Sample  Points 
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reasonable  control  effort.  The  disturbance  rejection  of  this  controller  is  good  when  the 
frequency  of  the  disturbance  is  less  than  the  order  of  the  frequency  of  the  sampler. 
Finally,  robustness  characteristics  for  this  Pl-controller  were  found  to  be  adequate. 
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Chapter  Five.  Fourth-Order  System  Design  With 

Aliasing 

5.0  Introduction 

With  the  success  of  implementing  aliasing  as  a  design  tool  found  in  chapters  three  and 
four  of  this  dissertation  one  may  ask  if  the  theory  can  be  extended  to  higher-order 
systems,  and  if  this  extension  will  be  of  any  benefit.  Therefore,  a  limited  analysis  of  a 
fourth-order  system  follows,  demonstrating  the  applicability  of  sample  period  as  a  design 
parameter  to  higher-order  systems  is  completed  in  this  chapter.  We  start  by  examining 
the  simple  fourth-order  system  with  two  sets  of  complex-conjugate  poles  and  no  zeros. 
Sample  frequency  ratios  are  such  that  the  high  frequency  mode  is  aliased,  and  a 
determination  of  “good”  frequency  ratios  is  made.  A  design  of  a  controller  for  the  low- 
frequency  mode  is  performed,  since,  by  selecting  the  correct  sampling  rate,  the  high- 
frequency  mode  can  be  virtually  eliminated  from  the  response.  Finally,  disturbance 
rejection  is  again  examined. 


5.1  The  Fourth-Order  Plant 


Consider  a  plant,  consisting  of  two  second-order  systems  found  in  series  with  each  other. 
This  plant  may  be  written  as 


GM  = 


+  a)„/)(5"  +  2C2to„2.s  +  (dj) 


[5.1] 


where  and  are  the  damping  ratios,  and  ro„,  and  ©„2  are  the  natural  frequencies  of 

the  two  modes  of  the  fourth-order  plant.  Locations  of  the  poles  of  this  plant  may  be 
found  using  Equation  3.1,  and  the  fourth-order  system  can  be  divided  by  s ,  and 
expanded  by  partial  fractions  to  enable  its  z-transform: 


GM 


( 


an 


«1 


a. 


Oj 


a. 


J  j  +  p,  S+P2  S  +  P3  s  +  p^ 


[5.2] 


where  the  constants  of  the  numerator  terms  can  be  found  from  a  series  of  algebraic 
manipulations  of  the  poles  using  the  following  constants: 


Ps  =  Pi  P2  +  P1P3  +  Pi  P4  +  P2P3  +  P2P 4  +  P3P4  . 

[5.3] 

P6  =  P2P3+P2P4  +  P3P4  ’ 

[5.4] 

Pl  =  PlP3+PlP4-^P3P4  . 

[5.5] 
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P&=PlP2  +  P\P4  +  P2P4  . 

[5.6] 

P9  =  P\P2  +  PiP3  +  P2P2  ’ 

[5.7] 

PlO  =  P\P2P3  +  P1P2P4  +  P\P3P4  +  P2P3P4  ’ 

[5.8] 

Pn=P2P3P4  > 

[5.9] 

A2  =  PiP3P4  . 

[5.10] 

P\3  =  P\P2P4  . 

[5.11] 

Pu  =  PlP2P3  ’ 

[5.12] 

do= - ^ - 

[5.13] 

P1P2P3P4 

j  _  -doPi 

Ui  —  , 

[5.14] 

P1-P2 

U2  —  , 

[5.15] 

P1-P2 

d  =^4-A 

“3  -  9 

[5.16] 

P1-P2 

II 

1 

1 

[5.17] 

^  1  ^2  * 

[5.18] 

+ 

T-^ 

II 

[5.19] 

+d^p>j  +  d^p^  , 

[5.20] 

^iPi  ^$P(i  > 

[5.21] 

1 

dg=Pg+d^PT+d^P^  , 

[5.22] 
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[5.23] 


II 

C 

VO 

[5.24] 

^12  ™  ^4  ^5^10  » 

[5.25] 

^13  ~  ”  ^5^11  » 

[5.26] 

^14  =  +  ^2^10  » 

[5.27] 

^15  —  <^3  +  ^2^11  • 

[5.28] 

With  these  constants  defined,  the  constants  of  Equation  5.2  can  be  defined  as 


^0  ^0  ’ 


Qj  —  d^2  ^13^4  J 


^2  ~  ^14  ^15^4  ’ 


^  =^10  ■'■^11^4  »2[nd 


[5.29] 


[5.30] 


[5.31] 


[5.32] 


^  —  ^oPio  ^iqPi3  ^uPn  ^uPm 

Pu  '^^nPi3  ■^^uAi  "^^isPn 


[5.33] 


Taking  the  z-transform  of  Equation  5.2,  applying  the  zero-order  hold,  gives 
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Hiz)  =  +  flZ"  +  fxZ  +  /o) 

(z + ^,)(z + + ^3)U + ^4) 


[5.34] 


where 


[5.35] 


/o  =  aohhhK  -  (hhhK  -  -  ^3^4  -  aAhh  > 


[5.36] 


/i  =  «o  (^A^3 + ^Ah + hhh + hhh) 

“1“  Cty  (^2^3^4  ”  ^2^3  "*  ^2^4  ”*  ^3^4  ^  ^2  (^1^3^4  ”  ^1^3  ^1^4  ^3^4  )  * 

+  a^{bp.Jy^-bp2  “^2^4)  +  ^4(^A^3  ~^A  “^^3  “^2^3) 


[5.37] 


A  =OQ(bb  +  bb  +  bb  +  bb  +  bb  +  bb) 

+  a,  +  ^2^4  +  ^3^4  -  ^2  ~  ^3  “  ^4  )  +  ^2  +  ^3^4  “  ^1  “  ^3  “  ^4  )  » 

+  a3(Z>A  +i>A  +^>2^4  -£?,  -i>2  “^4)  +  «4(^A  +^2^3  “^1  “^2  “^3) 

[5.38] 


/3  =flo(fc,  +Z?2  +^3  +fc4)  +  ai  A  +^3  +^4  “1)  +  «2A  +^3  +^4  ~1) 

+  ^(3  (fcj  +  ^>2  +  &4  ~  1)  +  ^3  (Z>j  +  ^>2  "I"  ^3  ~  1) 
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5.2  Design  of  Proportional  Controller  For  Fourth-Order 
System  In  Which  One  Mode  is  Aliased 

One  of  the  concerns  in  designing  systems  in  which  aliasing  is  present  is  the  effect  of 
aliasing  on  only  part  of  the  system’s  modes.  Suppose  a  system  exists  with  two  modes, 
the  fourth-order  system,  and  that  only  one  of  these  modes  is  aliased.  In  other  words, 
sampling  is  done  such  that  one  mode  is  aliased,  but  the  other  is  not.  The  following 
examines  the  effect  of  aliasing  on  such  a  system. 

Begin  by  allowing  the  controller  to  be  proportional  control  only.  Solving  for  F(z)  from 
Equation  3.3: 


F{z)  = 


Fp{z-\-b^){z  +  b2)iz  +  b.^){z  +  b^) 


z*  +  {KJ,  ■¥b,)z^  +  +b,)z^  +  (KJ,  +b,)z  +  KJ,  +b. 


[5.40] 


where 


b^-b^+b2+b^+b^  , 


h  =  +  hh  +  ^1^4  +  Kh  +  Kh  +  b^b^  , 


bj  =  bp2bj,  +  ^>,^2^4  +  b^b2b^  -I-  Z>2^3Z>4  , 


[5.41] 


[5.42] 


[5.43] 
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From  Equation  5.40,  the  root  loci  (z-plane  and  s*-plane)  of  a  fourth-order  system  in 
which  one  mode  is  aliased  can  be  constructed.  Additionally,  the  magnitude  and  phase  of 
the  frequency  response  functions  and  the  impulse  and  step  responses  will  be  plotted  for 
the  fourth-order  plant  of  Equation  5.1  combined  with  the  mixed  controller  function  of 
Equation  5.40.  Parameters  for  the  fourth-order  system  are 

a)„i  =1.0  Hz 
Cl  =0.5 
©„2  =  10.0  Hz 

^2=0.01  [5.46] 

Sampling  period  will  be  such  that  the  second  mode  is  aliased  but  the  first  mode  is  not. 

For  comparison  purposes.  Figures  5.1  through  5.3  show  the  root  loci,  frequency  response 
function  plots,  and  impulse  and  step  responses,  respectively,  when  the  frequency  ratio  is 
FR  =  100  — a  frequency  ratio  for  which  little  or  no  aliasing  is  present  for  either  mode. 

Symbols  on  the  root  loci  plots  are  as  follows.  Diamonds  and  pluses  (+)  represent  the 
zeros  and  poles  of  the  plant,  respectively.  Realize  again  that  these  zeros  and  poles  do  not 
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Phase  (deg)  Magnitude  (dB) 
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Response 


really  exist  in  the  z-plane — they  are  merely  shown  for  a  more  complete  understanding  of 
the  system.  Circles  and  “x’s”  represent  the  zeros  and  poles  of  the  mixed  term,  F(z) , 
respectively.  Finally,  triangles  and  asterisks  (*)  represent  the  zero-order  hold  terms. 

Examination  of  Figure  5.1  shows  the  effects  of  sampling  very  rapidly  with  respect  to 
natural  frequencies  of  the  system.  The  poles  and  zeros  of  the  discrete  portion  of  the  root 
locus  lie  very  close  to  z  =  1 .  Figure  5.2  shows  the  classical  frequency  response  function 
plot  of  a  fourth-order  system.  The  lightly  damped  poles  at  ®„2  =  I®-®  Hz  show  up  as  the 
sharp  peak  with  the  180°  phase  shift,  and  the  slower  poles  at  (0„,  =  1.0  Hz  show  up  with  a 

less  sharp  peak  and  phase  passing  through  -90°.  The  effect  of  the  higher  frequency  poles 
is  clearly  visible  in  the  impulse  and  step  responses  of  Figure  5.3. 

Now,  allowing  aliasing  of  the  second  (or  faster)  poles  select  the  sampling  frequency  with 
respect  to  the  damped  natural  frequency  of  these  poles: 

CO,  =  co^2  =  co„2Vi-C2^  •  [5.47] 

Here,  we  will  only  examine  frequency  ratios  of  FR  =  2.0  and  FR  =  1.0,  taking  advantage 
of  the  lessons  learned  with  the  second-order  system  designs.  These  frequency  ratios 
should  show  the  two  “extremes”  of  aliased  design,  ff  the  design  of  the  fourth-order 
system  with  one  aliased  mode  follows  that  of  the  second-order  system,  the  poles  and 
zeros  of  the  aliased  mode  should  lie  on  the  real  axis  in  the  z-plane. 
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Figures  5.4  and  5.5  show  this  is  indeed  the  case — at  these  frequency  ratios  the  aliased 
poles  lie  on  the  real  axis.  When  FR  =  2.0  the  poles  and  zeros  lie  on  the  negative  real 
axis,  and  when  Fi?  =  1.0  they  lie  on  the  positive  real  axis.  Notice  the  poles  will  move 
toward  the  zeros  of  the  plant,  represented  by  the  diamonds.  Again,  in  both  cases,  one 
pole  remains  fixed  at  the  location  of  one  zero  of  F(z) ,  which  is  the  same  location  as  one 
of  the  zeros  of  Hiz) ,  while  the  other  moves  toward  another  zero  of  H(z) .  The 
remaining  poles  of  F(z) ,  those  not  aliased,  move  in  a  large  circular  path  toward  their 
meeting  point  on  the  negative  real  axis,  where  one  will  then  move  toward  negative 
infinity  and  the  other  will  move  toward  the  remaining  zero  of  H(z) . 

Notice  that  as  the  frequency  ratio  decreases  (sample  period  increases)  the  non-aliased 
poles  move  away  from  the  point  z  =  1 ,  Again,  this  is  as  expected.  The  same 
relationships  between  frequency  ratio  and  locations  in  the  s*-plane  exist  for  the  aliased 
poles  as  were  found  in  chapter  four:  at  FR  =  2.0  the  aliased  poles  lie  on  the  band 
boundaries  and  at  FR  =  1.0  they  lie  on  the  mid-band  lines.  Another  thing  to  note  is  the 
pole  aliasing  of  the  non-aliased  poles — i.e.,  the  non-aliased  poles  are  repeated  in  each 
band  infinitely  in  the  s*-plane. 

At  a  frequency  ratio  of  FF  =  1 .0  the  poles  and  zeros  of  the  aliased  mode  lie  on  the 
positive  real  axis  of  the  z-plane,  and  it  appears  they  are  clustered  together  at  the  same 
point.  This  also  appears  to  be  the  case  in  the  s*-plane.  A  closer  examination  shows  there 
is  some  separation,  as  Figure  5.6  shows.  This  appears,  however,  to  be  very  nearly  a  pole- 

zero  cancellation,  a  fact  whose  significance  will  become  apparent  soon. 
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Figure  5.4  Root  Loci  of  Fourth-Order  System  with  One  Aliased  Mode 

FR  =  2.0,Kp=1.5 
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Z-Plane  Root  Locus;  FR  =  1;  Kp  =  1 .5 


S*“Plane  Root  Locus 


Figure  5.6  Close-Up  of  Aliased  Pole  and  Zero  Locations  for  FR  =  1,0 
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Figures  5.7  and  5.8  show  the  magnitude  and  phase  plots  of  the  discrete-to-continuous 
frequency  response  functions  for  frequency  ratios  of  FR  =  2.0  and  FT?  =  1 .0 , 
respectively.  At  a  frequency  ratio  of  FR  =  2.0 ,  the  second  mode  sits  right  at  the  folding 
frequency.  Thus,  its  effect  is  seen  in  the  magnitude  plot  as  a  sharp  peak  with  a  magnitude 
somewhat  reduced  from  that  of  when  it  was  not  aliased  (Figure  5.2).  For  a  frequency 
ratio  of  FF  =  1 .0 ,  however,  the  sampling  frequency  is  equal  to  the  damped  natural 
frequency  of  the  second  mode,  and  the  effects  of  the  zero-order  hold  cancel  this  mode. 
Notice  the  pole  aliasing  of  the  low  frequency  mode  as  it  is  “folded”  above  the  folding 
frequency  to  9.0  Hz  and  repeated  in  the  other  frequency  bands. 

Figures  5.9  and  5.10  show  the  impulse  and  step  responses  of  the  fourth-order  system  with 
one  aliased  mode.  Here,  the  effect  of  allowing  more  aliasing  is  clearly  seen  when 
compared  to  the  impulse  and  step  responses  of  Figure  5.3.  At  a  frequency  ratio  of 
FR  =  2.0 ,  seen  in  Figure  5.9,  both  the  impulse  and  step  responses  show  improvement 
over  their  counterparts  in  Figure  5.3.  The  magnitudes  of  the  higher-order  oscillations  are 
less  and  the  oscillations  themselves  seem  to  settle  sooner,  although  the  overall 
magnitudes  of  both  responses  are  larger  than  for  the  continuous  case. 

Figure  5.10  shows  these  responses  for  a  frequency  ratio  of  FF  =  1.0 .  As  predicted  by 
both  the  root  loci  plots  and  the  Frequency  Response  Function  plots,  the  effects  of  the  high 
frequency  mode  are  practically  invisible.  The  responses  of  Figure  5.10  appear  to  be  those 
of  a  second-order  system  only,  although  closer  examination  shows  some  of  the  higher- 
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Figure  5.9  Impulse  and  Step  Responses  of  Fourth-Order  System  with 
One  Aliased  Mode;  FR  =  2.0,  Kp  =  1.5 
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Response 


0.5 


Impulse  Response  at  WsA/Vd  =  1 ;  Kp  =  1 .5 


Step  Response  at  Ws/Wd  =  1 ;  Kp  =  1 .5 


Figure  5.10  Impulse  and  Step  Responses  of  Fourth-Order  System  with 
One  Aliased  Mode;  FR  =  1.0,  Kp  =  1.5 
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order  mode  effects  are  still  present.  However,  the  responses  of  Figure  5.10  are  clearly 
superior  to  those  of  Figures  5.3  and  5.9. 

Again,  by  allowing  aliasing  the  system’s  stability  has  actually  increased,  and  in  the  case 
of  the  fourth-order  system  the  high  frequency  mode  has  been  rendered  virtually 
ineffective.  This  leads  to  the  question  of  designing  a  controller  for  the  remaining  second- 
order  response.  Will  design  of  a  controller  for  this  system  in  which  the  higher-order 
mode  is  aliased  to  the  point  that  it  no  longer  affects  the  response  be  possible? 


5.3  Design  of  Proportional-Integral-Derivative  Controller 
For  Fourth-Order  System  In  Which  One  Mode  is  Aliased 


With  the  higher-order  mode  aliased  as  developed  above  (FR  =  1.0)  design  of  a  controller 
for  the  remaining  second-order  response  will  be  examined.  The  design  will  be  a  simple 
proportional-integral-derivative  (PID)  controller,  designed  in  the  z-plane,  to  improve  that 
response. 

The  equation  for  a  PID  controller  in  the  z-plane  is  defined  as 


G,(z)  =  K,+K,^^  +  K, 


z-l 

z 


K^z^  +  K^z  +  K, 
ziz-1) 


[5.48] 
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where 


K^=K^+K,  +  Kj  ,and  [5.49] 

K,=-K^-2K,  .  [5,50] 

Equation  3.3  becomes 

p/  N  _  (^1^^  +  +  Kj)(z  +  b^ )iz + b^)(z + b^)iz + b^) 

654^0  9  I  1  I 

2  +85Z  +84Z  +83r  +  82Z  +81Z+80 

where,  recalling  Equations  5.35  through  5.39  and  Equations  5.41  through  5.45, 

8o  =  KjKJo  ,  [5.52] 

8i=K^{KJ^  +  KJ^)-b^  ,  [5.53] 

82  =  K^{KJf^  +  KJ^-^KJ^)  +  b^-b^  ,  [5.54] 

82=^K^{KJ^  +  KJ^  +  KJ^)  +  bT-b^  ,  [5.55] 

8a,  =  Kg (KJ2  +  ^2/3) +h-h  >  and  [5.56] 
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g,=K^KJ,+b,  . 


[5.57] 


Using  these  equations  a  PID  controller  was  designed  for  the  fourth-order  system  in  which 
the  high-frequency  mode  was  aliased  by  setting  the  frequency  ratio  to  F/?  =  1.0 .  The 
results  of  this  design  are  shown  in  Figures  5.11  and  5.12.  Figure  5. 1 1  shows  the  root  loci 
for  this  system.  The  added  poles  and  zeros  from  the  controller  are  visible,  and  their  effect 
on  the  “fixed”  poles  and  zeros  of  the  aliased  mode  is  inconsequential.  Figure  5.12  shows 
the  impulse  and  step  responses  of  the  system  with  the  aliased  mode.  Here  the  effect  of 
aliasing  on  the  higher-order  mode  is  clearly  visible,  as  its  contribution  to  the  responses  is 
nearly  negligible.  Very  little  of  the  effects  are  visible  in  the  impulse  response,  as  well  as 
in  the  step  response,  although  this  high-order  response  is  visible  at  the  response  settles  to 
the  final  desired  value. 

With  the  design  of  the  PED-controller  for  the  fourth-order  system  accomplished,  an 
examination  of  the  controller’s  ability  to  reject  disturbances  in  again  in  order.  Figure 
5.13  shows  the  system’s  responses  to  a  step  disturbance  at  the  sample  point.  Again,  the 
system  is  able  to  reject  the  disturbance  once  the  controller  “sees”  it  at  the  following 
sample  point.  The  settling  of  the  system  to  the  correct  steady-state  value  can  be  seen 
quite  clearly. 

Figure  5.14  shows  the  same  characteristic  for  the  system  when  the  step  disturbance 
begins  between  sample  periods.  Finally,  Figure  5.15  shows  the  system’s  responses  to  a 
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Z-Plane  Root  Locus,  FR  =  1 


~1  -0.5  0  0.5  1 


Kp=:0.7;  Kl  =  0.6;  Kd=1.4 


S*-Plane  Root  Locus 


Figure  5.11  Root  Loci  of  Fourth-Order  System  with  One  Aliased  Mode 
FR  =  1.0,  Kp  =  0.7,  Ki  =  0.6,  Kd  =  1.4 
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Impulse  Response,  FR  =  1 


Step  Response;  Kp  =  0.7;  Ki  =  0.6;  Kd  =  1 .4 


Figure  5.12  Impulse  and  Step  Responses  of  Fourth-Order  System 
with  One  Aliased  Mode;  FR  =  1.0,  Kp  =  0.7,  Kj  =  0.6,  Kd  =  1.4 
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Impulse  Response,  FR  =  1 


Figure  5.13  Fourth-Order  System  Unit  Step  Disturbance  Rejection 
Disturbance  at  the  Sample  Point 
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Normalized  Time  (t/Ts) 
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random  disturbance.  Again,  the  inability  of  this  design  to  control  anything  changing 
faster  than  the  sampler’s  sample  rate  is  evident.  Following  the  disturbance,  however,  the 
system  returns  to  the  desired  steady-state  value  as  required. 

Design,  then,  for  a  system  in  which  one  mode  is  aliased  can  be  accomplished  using 
traditional  methods  for  the  remaining  un-aliased  modes,  with  the  aliased  mode’s  effects 
being  negated  simply  by  selecting  the  correct  sampling  frequency.  This  selection  is  done 
such  that  the  ratio  of  the  sampling  frequency  to  the  damped  natural  frequency  of  the  high 
frequency  mode  is  equal  to  one.  By  so  doing,  the  designer  may  effectively  eliminate  the 
response  of  that  mode. 


5.4  Chapter  Five  Summary 

In  this  chapter  the  design  of  a  controller  for  a  fourth  order  system  in  which  the  high- 
frequency  mode  is  aliased  was  designed.  Once  again,  the  frequency  ratio  of  FR  =  \  .0 
was  found  to  have  favorable  effects  on  the  control  of  the  aliased  mode.  In  this  case,  the 
mode  was  basically  suppressed,  and  design  for  the  low-frequency  mode  was 
accomplished  using  traditional  design  methods  with  little  response  resulting  from  the 
aliased  high-frequency  mode.  Disturbance  rejection  for  this  controller  is  as  was  found  for 
the  other  designs  in  which  aliasing  was  allowed:  good  for  disturbances  whose  frequencies 
are  below  the  sampling  frequency  and  poor  for  broadband  disturbances. 
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Response 


Impulse  Response,  FR  =  1 


Step  Response;  Kp  =  0.7;  Ki  =  0.6;  Kd  =  1 .4 


Figure  5.15  Fourth-Order  System  Random  Input  Disturbance  Rejection 
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Chapter  Six.  Experimental  Verification 

6.0  Introduction 

Given  the  tools  and  processes  developed  theoretically  in  chapters  three  through  five,  it 
becomes  important  to  verify  these  design  methods  experimentally.  Experimental 
verification  will  be  accomplished  through  applying  the  concepts  of  design  to  simple 
electronic  circuits.  The  comparisons  of  experimental  results  to  theoretical  model  results 
are  excellent. 


6.1  Second-Order  Circuit 


To  demonstrate  the  accuracy  of  the  design  method  developed  in  chapters  four  and  five 
the  circuit  show  in  figure  6.1  was  constructed.  This,  then,  becomes  the  plant  for  the 
discrete-to-continuous  system.  Given  the  values  of  resistance  and  capacitance  seen  in 
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figure  6.1  the  natural  frequency,  system  damping  and  system  gain  values  are  as  found  in 


table  6.1. 


Table  6.1  Second-Order  System  Parameters  From  Calculations 


System  Parameter 

Value 

Natural  Frequency 

0)„  =106.2116  Hz 

Damping  Ratio 

C=  0.08294 

System  Gain 

^  =  0.8106 

With  the  circuit  built,  a  swept  sine  wave  was  introduced  into  the  circuit.  The  sample  rate 
was  set  at  4000  Hz  and  a  Hamming  window  was  applied  to  both  the  input  and  output. 
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The  number  of  data  points  was  16,384  samples,  with  30  averages  of  the  resulting  input 
and  output  auto  spectra  taken.  Coherence  over  the  range  of  interest,  0.1  to  200  Hz,  was 
essentially  equal  to  one,  with  the  lowest  value  at  a  single  point  equal  to  0.982,  The 
values  of  magnitude  and  phase  from  the  resulting  frequency  response  function  were 
saved  to  file. 

Using  Matlab,  a  least  squares  error  fit  between  the  experimental  data  and  model 
parameters  was  accomplished  to  determine  experimental  natural  frequency,  damping 
ratio,  and  system  gain.  The  results  of  this  curve  fit  are  shown  in  figure  6.2,  the  solid  line 
representing  the  experimental  data  and  the  dashed  line  representing  results  of  the  curve 
fit — ^the  frequency  response  function  using  model  values  from  the  curve  fit.  Values 
found  from  this  curve  fit  are  as  found  in  table  6.2. 


Table  6.2  Second-Order  System  Parameters  From  Curve  Fit 


System  Parameter 

Value 

Natural  Frequency 

Q)„  =98.562  Hz 

Damping  Ratio 

C  =  0.1188 

System  Gain 

Ji:  =  0.8085 
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Magnitude  Curve  Fit,  K = 0.8085;  wn  =  98.562;  zeta  =  0.1 188 


Frequency  (Hz) 


Figure  6.2  Least  Squares  Error  Curve  Fit  of  Model  to  Experimental  Data 


Finally,  an  impulse  was  applied  to  the  system  and  10  averages  were  taken.  From  the 
impulse  responses,  the  natural  frequency  and  damping  ratio,  using  logarithmic 
decrement,  were  calculated.  Results  are  listed  in  table  6.3. 
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Table  6.3  Second-Order  System  Parameters  From  Impulse  Responses 

I  “■  I  1!“!;  1  I 


System  Parameter 

Value 

Natural  Frequency 

co„  =  92.7803  Hz 

Damping  Ratio 

§ 

o 

II 

System  Gain 

N/A 

The  values  in  tables  6.2  and  6.3  vary  from  those  of  table  6.1  because  of  differences  in 
theoretical  and  actual  systems.  Physical  systems,  although  close,  rarely  exhibit  exact 
calculated  values.  For  the  experimental  verification  to  follow  values  found  from  the 
impulse  response  will  be  used. 


6.2  Experimental  Verification  of  Second-Order  System  with 
Proportional  and  Proportional-Integral  Control 

With  the  second-order  system  parameters  experimentally  determined,  the  validity  of  the 
model  developed  in  chapters  three  and  four  can  now  be  determined.  Using  the  same 

model  of  chapter  three,  values  for  the  sampling  frequency  (Q)^  =  and  the 

proportional  gain  were  calculated  and  used  in  the  program  for  the  digital  controller. 

The  A/D  card  used  for  data  acquisition  was  simple  with  limited  options  and  a  maximum 

sample  rate  of  2000  Hz.  Sample  rate  was  “infinitely”  variable,  with  a  change  of  sample 

period  possible  to  the  ninth  decimal  place.  Because  of  the  magnitude  of  the  system 

natural  frequency  and  the  limitations  of  the  A/D  sample  rate,  however,  only  20 
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“continuous”  points  were  taken  between  the  sample  points  sent  to  the  controller.  Even 
so,  this  is  a  sufficient  number  of  points  to  map  the  continuous  response  of  the  sampled 
data  system. 

The  results  of  the  first  series  of  runs  are  seen  in  figure  6.3.  Only  proportional  control  was 
employed.  An  average  of  five  data  sets  was  used  to  eliminate  noise  from  individual  runs 
of  the  system.  The  match  between  the  model  results,  the  solid  line,  and  the  continuous 
output  of  the  system,  the  dashed  line,  are  not  very  good,  indicating  a  problem  with  the 
model  parameters.  Examination  of  the  sample  points,  the  asterisks  on  figure  6.3,  shows  a 
difference  between  the  sampling  frequency  and  the  damped  natural  frequency  of  the 
system  with  the  experimental  frequency  being  lower.  Damping  of  the  model  also  appears 
to  be  greater  than  that  of  the  experimental  system.  Decreasing  the  natural  frequency 
found  from  the  impulse  response  by  0.51%  and  increasing  the  damping  ratio  by  15.88% 
gives  the  results  of  figure  6.4.  Here,  the  correspondence  between  experimental  and 
model  data  is  much  improved,  and  one  can  see  that  with  even  further  changes  made,  an 
exact  match  between  the  model  and  experimental  data  could  be  obtained. 

Using  these  same  values  of  C0„=  92.31  Hz.  and  ^  =  0.1211  in  the  proportional-integral 

controller,  the  resulting  figure  is  shown  in  figure  6.5.  One  can  see  again  the  good  match 
between  the  model  and  the  experimental  data. 
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Step  Response  at  Ws/Wd  =  1.0;  Kp  =  1 .1036;  wn  =  98.562;  zeta  =  0.1 188;  K  =  0.80854 


0  0.005  0.01  0.015  0.02  0.025  0.03  0.035  0.04  0.045  0.05 

Time  (sec) 


Figure  6.3  P-Controller  Experimental  Step  Response 
FR=1.0,  ©„=  98.562 Hz,  ^  =  0.1188,  ^:  =  0.8085, =  1.1036 
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step  Response  at  FR  =  1 .0;  Kp  =  1 .0733;  Ki  =  1 .2368;  wn  =  92.31 ;  zeta  =  0.1 21 1 ;  K  =  0.80854 
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Time  (sec) 


Figure  6.5  Pl-Controller  Experimental  Step  Response 
=  1.0,  co„=  92.31  Hz,  C  =  0.1211,  .«:=  0.8085, 1.0733,  Ki=  1.2368 
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Examination  of  figures  6.4  and  6.5  show  that,  indeed,  the  experimental  method  matches 
the  mathematical  model  developed  in  chapters  three  and  four.  It  also  shows  the  difficulty 
in  pole  placement,  as  small  changes  in  the  natural  frequency  or  damping  ratio  causes  the 
system  to  respond  differently  than  predicted.  This  can,  however,  be  predicted  using  the 
model. 


6.3  Fourth-Order  Circuit 


To  demonstrate  the  accuracy  of  the  model  of  the  fourth-order  system  developed  in 
chapter  five,  the  circuit  shown  in  figure  6.6  was  constructed.  The  fourth-order  system 
consists  of  two  second-order  circuits  similar  to  that  of  section  6.1,  To  get  the  fourth  order 
system,  the  signal  of  eom  from  the  top  system  of  figure  6.6  is  input  as  ei„  to  the  bottom 
system  of  figure  6.6.  Given  the  values  of  measured  resistance  and  capacitance  seen  in 
figure  6.6,  the  natural  frequency,  system  damping,  and  system  gain  values  were 
calculated.  These  calculated  values  are  shown  in  table  6.4,  below.  These  values  were 
selected  such  that  one  mode  was  easily  aliased  to  the  frequency  range  of  the  other.  Also 
the  aliased  mode’s  damping  ratio  is  such  that  it  is  very  small  in  comparison  to  that  of  the 
non-aliased  mode  so  that  it’s  effect  might  be  easily  observed  in  the  experimental  data. 
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Table  6.4  Fourth-Order  System  Parameters  From  Calculations 


System  Parameter 

Value 

Natural  Frequency 

tD„,  =106.21  Hz 

0)„2  =9.343  Hz 

Damping  Ratio 

C,  =0.00999 

C,  =0.1318 

System  Gain 

^  =  0.8097 

IT  =  1.0022 

With  the  two  second-order  circuits  built,  a  digital  impulse  was  introduced  into  each  and 
the  natural  frequencies  and  damping  ratios  of  each  of  the  circuits  were  determined  using 
the  logarithmic  decrement  method.  The  results  of  these  experiments  are  shown  in  table 
6.5  below.  Again,  the  physical  parameters  fairly  closely  match  the  calculated  parameters 


Table  6.5  Fourth-Order  System  Parameters  From  Impulse  Response 


System  Parameter 

- - - ^ - x  _ 

Value 

Natural  Frequency 

co„,=  108.71  Hz 

G)„2  =10.10  Hz 

Damping  Ratio 

=0.01498 

C,  =0.1330 

Using  these  values  in  the  experimental  setup,  the  fourth-order  system  was  tested  with  a 
digital  proportional-integral-derivative  (PID)  controller.  The  results,  shown  in  figure  6.7 
show  excellent  agreement  of  experimental  data,  the  dashed  line,  to  model  data,  the  solid 
line.  The  fourth-order  model  accurately  predicts  experimental  results.  Close 
examination  shows  the  response  of  the  aliased  mode,  as  well,  although  because  of  the 
prudent  selection  of  sampling  rate  the  aliased  mode  introduces  very  little  response. 
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step  Response;  Kp  =  0.2;  Kd  =  0.4;  Ki  =  0.1 
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Time  (sec) 


Figure  6.7  PID-Controller  Experimental  Step  Response 
Kp  =  0.2,Kd  =  0.4,Ki  =  0.1 
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To  show  the  effects  of  designing  for  the  lower  order  mode,  and  how  this  will  not  affect 
the  effects  of  the  aliased  mode  on  the  response,  different  values  for  the  gains  were 
selected.  The  result  is  shown  in  figure  6.8.  Again,  experimental  response  matches 
closely  to  that  predicted  by  the  model.  Moreover,  the  aliased  mode  has  very  little  effect 
on  the  design  of  control  for  the  non-aliased  mode. 


6.4  Chapter  Six  Summary 


Experimental  verification  of  the  model  developed  in  this  dissertation  was  accomplished 
in  this  chapter.  Experimental  results  show  excellent  agreement  with  the  model.  They 
also  show  the  sensitivity  of  the  controller’s  response  to  accurate  modeling  of  the  plant’s 
dynamics. 
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step  Response;  Kp  =  0.05;  Kd  =  2.1 ;  Ki  =  0.5 


Figure  6.8  PID-Controller  Experimental  Step  Response 
Kp  =  0.05,  Kd  =  2.1,  Ki  =  0.5 
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Chapter  Seven.  Conclusions  and  Recommendations 


7.0  Conclusions 

This  research  effort  begins  by  developing  a  model  in  which  aliasing  is  allowed.  The 
model  developed  includes  both  discrete  and  continuous  sections  of  the  system.  Unlike 
the  traditional  sampled-data  systems  in  which  analysis  is  completed  in  either  the  discrete 
or  continuous  domains,  this  model  allows  analysis  for  systems  in  which  both  discrete  and 
continuous  models  exist  concurrently.  The  resulting  transfer  function  is  called  a  discrete- 
to-continuous  transfer  function. 

In  the  limit,  as  the  sample  period  of  the  discrete  portion  of  the  model  goes  to  zero,  the 
model  becomes  a  continuous  model,  as  would  be  expected.  To  mathematically  obtain 
this  limit,  however,  the  system  must  be  physically  realizable,  or  the  order  of  the 
numerator  of  the  system  must  be  equal  to  or  less  than  the  order  of  the  denominator  of  the 
system.  This  condition  holds  true  for  all  physical  systems.  The  limit,  as  the  sample 
period  increases  to  infinity,  is  zero,  signifying  a  system  in  which  no  feedback  is  ever 


221 


seen.  Again,  this  is  the  expected  result,  since  as  the  sample  period  goes  to  infinity  in  a 
discrete  system  no  feedback  occurs.  The  discrete-to-continuous  model,  then,  accurately 
models  both  continuous  and  discrete  systems  in  the  limit  as  the  sample  periods  go  to  zero 
and  infinity,  respectively. 

To  aid  in  analysis,  a  “new”  plane  was  introduced — the  s*-plane.  This  is  the  plane  in 
which  the  continuous  poles  and  zeros  exist,  and  to  which  the  poles  and  zeros  from  the 
discrete  portion  of  the  system  are  transformed,  repeating  infinitely  through  all  bands.  It  is 
the  plane  in  which  the  only  appropriate  controller  function  can  be  designed.  This  unique 
controller  function  F  *  (s)  is  the  only  function  representing  the  unique  time  function 
/  *  (0  that  excites  the  plant  only  at  the  sample  points.  The  time  function  /  *  (t)  is 
modeled  as  a  series  of  impulses  at  the  sample  points  and  zero  everywhere  else. 

The  s*-plane  allows  one  to  examine  the  effects  of  aliased  poles  and  zeros — ^poles  and 
zeros  located  in  bands  outside  the  primary  band — and  their  effects  on  the  response  of  the 
discrete-to-eontinuous  system.  This  definition  of  aliased  poles  (and  zeros)  in  the  s*-plane 
is  perhaps  different  than  the  traditional  view  of  aliasing. 

In  the  s*-plane,  a  set  of  zeros  from  the  discrete  controller  will  always  cancel  the  poles 
from  the  continuous  plant.  These  zeros,  however,  are  repeated  in  each  of  the  infinite 
bands,  other  than  the  band  in  which  the  continuous  poles  are  located,  and  as  a  result,  the 
continuous  system’s  response  between  sample  periods  is  that  of  the  open-loop  damped 
natural  frequency  of  the  plant  given  a  step  input  and  the  initial  conditions  found  at  each 
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sample  point.  This  continuous  response  between  samples  cannot  be  completely 
controlled.  It  can,  however,  be  used  in  designing  an  optimal  control  of  the  system. 

Examination  of  the  responses  of  the  system  shows  the  discrete  response  follows  that 
which  would  be  expected  based  on  the  pole  placement  in  the  z-plane.  When  the  poles  of 
the  controller  term  are  located  on  the  negative  real  axis,  the  discrete  response  is  an 
exponentially  decaying  oscillation.  Likewise,  when  located  on  the  positive  real  axis,  the 
response  is  exponentially  decaying  with  no  oscillation. 

Because  of  the  discrete-to-continuous  nature  of  the  developed  model,  the  sample  rate  or 
sample  period  becomes  a  design  parameter.  Allowing  aliasing,  or  slowing  the  sample 
rate  such  that  aliasing  occurs,  has  a  stabilizing  effect  on  a  second-order  plant.  This 
stabilizing  effect  is  seen  in  the  increased  gain  margins  at  certain  frequency  ratios,  as  well 
as  in  improved  impulse  and  step  responses.  The  frequency  ratio  is  defined  as  the  ratio  of 
the  sampling  frequency  to  the  damped  natural  frequency  of  the  aliased  mode. 

Sample  rates  were  found  such  that  the  entire  root  locus  of  a  second-order  system  with 
proportional  control  lies  on  the  real  axis  in  the  z-plane.  These  sample  rates  correspond  to 
frequency  ratios  of 


n 


n  =  integer. 


[7.1] 
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At  these  frequency  ratios,  two  of  the  controller  zeros  collapse  to  a  single  point  in  both  the 
z-  and  s*-planes,  causing  one  of  the  controller  poles  to  remain  at  the  same  point,  no 
matter  the  controller  gain.  At  this  frequency  ratio,  the  second  pole  of  this  pair,  the 
stationary  pole’s  “mate,”  travels  only  along  the  real  axis  in  the  z-plane  as  the  controller 
gain  is  increased. 

When  n  =  odd  is  used  in  the  frequency  ratio,  part  of  the  root  locus  lies  on  the  positive  real 
axis  and  part  on  the  negative  real  axis  of  the  z-plane.  This  corresponds  to  part  lying  on 
the  band  boundaries  and  part  lying  on  the  mid-band  lines  in  the  s*-plane.  The  impulse 
and  step  responses  seen  at  these  frequency  ratios  show  the  most  desired  characteristics  of 
the  aliased  responses.  Selecting  a  frequency  ratio  of  FR  =  \  .0  gives  the  best  response  of 
this  class  of  responses,  giving  only  one  complete  cycle  of  the  damped  natural  response  of 
the  open-loop  plant. 

Additionally,  when  the  frequency  ratio  is  set  to  one,  placement  of  the  non-stationary  pole 
mentioned  above  can  be  such  that  it  lies  at  the  origin  of  the  z-plane,  which  corresponds  to 
negative  infinity  in  the  s*-plane.  With  the  controller  designed  such  that  the  pole  is 
located  here,  a  deadbeat  response  results.  Using  a  simple  proportional-integral  (PI) 
controller,  the  perfect  deadbeat  response  is  obtained  with  very  little  problem:  within  one 
sample  period,  the  step  response  is  at  its  final  desired  value.  For  this  response,  the 
proportional  gain,  ,  is  a  function  only  of  the  damping  and  the  integral  gain.  A", ,  is 

always  equal  to  one.  Experimentation  shows  the  model  to  produce  accurate  results,  with 
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experimental  data  matching  the  response  predicted  by  the  model.  This  deadbeat 
response,  then,  is  obtainable  in  physical  systems. 

The  deadbeat  response  of  the  discrete-to-continuous  system  is  different  than  conventional 
deadbeat  responses  in  that  here,  with  aliasing,  we  are  taking  advantage  of  the  system’s 
natural  response,  phasing  the  control  such  that  the  response  is  exactly  cancelled  following 
one  sample  period.  It  does  not  force  response  beyond  the  system’s  own  natural  response. 

The  second-order  system  with  proportional  integral  control  has  excellent  disturbance 
rejection  capability,  if  the  frequency  of  the  disturbance  is  lower  than  the  frequency  of  the 
controller.  The  opposite  is  also  true — ^this  controller  has  poor  disturbance  rejection  when 
the  disturbance’s  frequency  is  broadband. 

Robustness  of  the  system  is  examined  and  the  system  is  found  to  be  fairly  robust. 
Although  the  deadbeat  response  is  lost,  and  this  was  to  be  expected,  the  system  settles  to 
final  values  within  six  sample  periods  given  a  5%  change  in  natural  frequency  or 
damping  of  the  plant.  Clearly,  the  plant  to  be  controlled  here  must  be  well-identified  and 
modeled. 

One  of  the  drawbacks  of  this  system  is  the  great  overshoot  present  in  the  deadbeat 
response.  Another  is  its  poor  following  of  non-stationary  input  signals.  The  response  of 
the  discrete-to-continuous  system  to  sinusoidal  inputs  is  poor,  at  best,  unless  the  input 
sinusoid  has  a  very  low  frequency  relative  to  the  sampling  frequency.  This  type  of 
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control,  then,  would  not  be  appropriate  for  rapidly  changing  systems  with  respect  to  the 
sampling  frequency. 

The  power  of  this  type  of  control  can  be  seen  when  considering  systems  in  which 
accurate  position  is  required,  with  minimum  time,  and  with  little  care  for  overshoot  of  the 
system.  Such  systems  might  include  the  positioning  of  a  crane  or  a  disk  drive  head  in  a 
minimum  amount  of  time.  These  systems  need  accurate  positioning,  have  allowances  for 
overshoot,  and  require  low  settling  times.  The  aliased  discrete-to-continuous  system 
developed  in  this  dissertation  will  meet  these  requirements. 

The  fourth-order  system  in  which  one  of  the  modes  is  aliased  was  examined.  It  was  seen 
that  by  selecting  the  sample  rate  such  that  the  frequeney  ratio  was  F/?  =  1.0 ,  with  respect 
to  the  natural  frequency  of  the  mode  to  be  aliased,  the  effects  of  the  aliased  mode  are 
minimized.  In  fact,  the  response  of  this  aliased  mode  is  almost  not  seen.  Control  for  the 
non-aliased  mode  can  be  accomplished  using  traditional  methods,  with  the  resulting 
controller  affeeting  the  aliased  mode  very  little.  Again,  disturbance  rejeetion  is  good  for 
disturbances  with  frequencies  below  that  of  the  sampler.  This,  then,  becomes  an  effective 
means  for  controlling  aliased  modes— simply  select  the  sampling  rate  to  be  equal  to  that 
of  the  aliased  mode. 
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7.1  Recommendations  for  Future  Study 


To  further  improve  the  design  methods  found  in  this  research,  additional  studies  should 
be  performed  with  respect  to  physical  systems.  This  effort  examined  only  electric 
circuits  to  verify  the  theoretical  model.  Systems  for  which  this  method  would  be  useful 
should  be  explored  and  implemented. 

The  further  effects  of  zero  placement  in  the  s*-plane  need  to  be  researched  as  well.  Zero- 
placement  was  not  a  parameter  of  the  design  in  the  theoretical  modeling.  The  researcher 
simply  selected  zero  locations  for  ease  of  experimental  verification.  Such  zero-placement 
may  in  some  manner  aid  or  detract  from  the  ideal  designs  examined  here. 

Further  research  should  also  aim  to  tie  this  theory  for  higher-order  systems  into  state- 
space  domain,  so  that  the  powerful  techniques  available  there  could  be  employed  in  the 
design  of  aliased  controllers  as  well. 

Another  area  of  further  research  should  be  the  examination  of  the  effects  of  time  delays 
in  systems  in  which  aliasing  is  allowed.  Are  the  time  delays  detrimental  as  conventional 
theory  would  suggest,  or  will  they  aid  in  stabilization  as  was  the  case  with  aliasing  only 
of  this  research  effort?  Richard  Lyman  [13]  demonstrated  time  delays  can  sometimes 
stabilize  the  system. 
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One  of  the  first  assumptions  in  performing  this  research  was  that  of  an  ideal  sampler — 
one  in  which  the  period  required  to  sample  the  signal  is  very  much  smaller  than  the 
sample  period  of  the  sampler.  Another  avenue  of  study  would  be  examination  of  the 
effects  of  relaxing  this  assumption.  As  the  sampler  takes  a  more  significant  amount  of 
time  compared  to  the  sample  period,  the  ideal  sampler  degenerates  to  a  step-on,  step-off 
mechanism,  causing  some  phase  delay  and  smaller  magnitude  in  the  frequency  response 
function  (FRF)  as  the  frequency  increases.  The  results  become  the  multiplication  of  the 

FRF  by  a  sine  function — similar  to  multiplying  one  more  zero-order  hold  times  the  entire 
FRF. 

Additionally,  the  effects  of  optimizing  the  parameters  found  for  the  controllers  could  be 
examined.  Will  optimization  techniques  work?  And  what  are  their  effects  on  the  aliased 
modes,  either  the  deadbeat  response  found  for  the  second  order  system  or  the  aliased 
mode  found  for  the  fourth-order  system? 
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Appendix  A.  Theoretical  Review 


A.0  Introduction 

This  appendix  will  review  and  discuss  briefly  the  concepts  necessary  to  develop  the 
model  used  in  this  research.  The  review  will  include  the  mathematical  equations  of  the 
Laplace-  and  Z-transforms  and  the  natural  logarithm  of  a  complex  number,  a  brief 
explanation  of  root  loci  and  frequency  response  function  plot  and  their  uses,  the  concepts 
behind  and  the  assumptions  made  for  the  ideal  sampler,  the  result  of  sampling  on  a 
continuous  system,  and  a  review  of  the  zero-order  hold  and  its  affect  on  a  continuous 
system. 


A.1  The  Laplace  Transform 

When  dealing  with  continuous  systems  for  which  control  is  to  be  accomplished,  it  is  often 
easier  to  deal  in  the  s-  or  Laplace  domain,  rather  than  in  the  time  domain.  The  definition 
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of  the  one-sided  Laplace  transform  is  given  by  Zill  and  Cullen  [35:209]  as 


^{f{t)}  =  F{s)=re^f{t)dt  =  \im\le-^^f{t)dt  .  [A.l] 

b-^oo  JO 


The  one-sided  Laplace  transform  given  in  Equation  A.l  assumes  that  f(t)  is  zero  for  all 
time  less  than  zero. 


The  inverse  Laplace  transform  simply  returns  the  original  time  signal  [35:217]: 


mF{_s)}  =  f{t)  . 


[A.2] 


The  inverse  Laplace  transform  of  a  function  may  not  be  unique.  However,  if  /j  and  /j 
are  continuous  on  [0,  <»)  and 


[A.3] 


then  /i  =  /j  on  that  interval. 


Tables  of  Laplace  transforms  for  many  time  functions  can  be  found  in  any  advanced 
mathematics  text  [35:inside  front  cover],  [50:299-300].  The  Laplace  domain  is  also 
known  as  the  frequency  domain. 


A.2  The  Z-Transform 

When  working  with  sampled-data  systems,  a  transformation  between  the  sampled-time 
domain  and  z-domain  aids  in  analysis.  The  definition  of  the  one-sided  z-transform  is 
given  by  Jury  [3 1 :2]  as 


Z[f{nT,)}  =  F{z)  =  Y^f{nT,)z-’^  . 

n-0 


[A.4] 


Again,  tables  of  z-transforms  are  available  to  aid  in  this  calculation  [Linside  front  cover], 
[31:278-296],  [32],  [33:49-50].  The  inverse  z-transform  gives: 


Z-^{F{z))  =  hnT,)  . 


[A.5] 


However,  with  the  inverse  z-transform  one  must  be  careful  as  to  its  interpretation.  Ogata 
[33:69]  states,  “. .  .only  the  time  sequence  at  the  sampling  instants  is  obtained  from  the 
inverse  z-transform.  Thus,  the  inverse  z-transform  of  X  (z)  yields  a  unique  x{k) ,  but 
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does  not  yield  a  unique  x(t) Hence,  the  “hat”  notation  on  f{nT^) — a  reminder  that  it 
represents  any  one  of  an  infinite  number  of  continuous  time  functions,  fit) ,  with 
equivalent  values  at  the  sample  instants. 

Since  both  the  s-  and  z-transformations  are  representations  of  a  function  of  time  one  may 
question  whether  there  is  any  relation  between  the  two.  Indeed,  there  is.  Franklin, 
Powell,  and  Workman  [1:127]  explain  that  the  discrete  system  can  be  found  from  the 
continuous  system  by  mapping  the  s-plane  poles  into  the  z-plane  using  the  relation 

^  =  «  •  [A.6] 

One  must  be  careful  not  to  map  zeros  using  this  same  relationship,  as  the  result  is 
incorrect.  To  map  a  system  with  poles  and  zeros,  one  must  use  the  partial  fraction 
expansion  of  the  system,  expanding  the  ratio  of  polynomials  associated  with  the  pole-zero 
form  into  pole-residue  form.  Since  the  poles  found  in  either  the  pole-residue  form  or  the 
pole-zero  form  are  identical,  they  map  correctly  either  way.  The  zeros,  however,  do  not. 
This  will  be  further  explored  as  the  theory  for  this  research  effort  is  developed. 
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A.3  Natural  Logarithm  of  a  Complex  Number 


The  natural  logarithm  of  a  complex  number,  z ,  is  defined  as  [35:962-963] 


ln(z)  =  ln|z|  + j(0  +  2n7t)  n  =  0,±l,±2,...  ,  [A.7] 


where  0  is  the  argument  of  z  [35:937]: 


0  =  arg(z)  =  tan 


-1 


2 

X 


[A.8] 


given  that 


z  =  x+/y 


[A.9] 


Because  of  the  nature  of  the  tangent  function,  the  principal  value  of  the  argument  of  z  is 
defined  to  lie  in  the  interval  (-7U,  Tt]. 
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A.4  Root  Loci 


Given  a  system  or  plant,  Gp(s) ,  the  plot  of  the  locus  of  its  closed-loop  pole  locations  in 

the  s-plane  as  a  parameter  K  changes  gives  what  is  known  as  the  root  locus  of  the 
system.  It  is  the  graph  of  all  possible  roots,  assuming  single-loop  negative  feedback  of  a 
system,  of 

\  +  KG^{s)Gpis)  =  0  [A.  10] 

as  K  changes  [51:244-246].  Here  the  gain  K  is  shown  as  the  varying  parameter,  and 
Gcis)  is  the  continuous  controller.  The  root  locus  can  be  used  to  design  feedback 

characteristics  of  a  system  in  either  the  continuous  or  discrete  domains.  The  equation  for 
a  discrete  domain  root  locus  is  given  [1:227-234]  as 

\-k-KG^{z)Gp{z)  =  0,  [A.ll] 


where  G^.(z)  is  the  discrete  or  digital  controller  and  G^  (z)  is  the  discrete  representation 


of  the  plant.  Again,  K ,  here  a  proportional  gain,  is  the  parameter  of  interest  that  will  be 
varied  to  develop  the  root  locus  plot. 
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Stability  in  the  continuous  domain  corresponds  to  pole  locations  located  in  the  left-half  of 
the  s-plane  ( s  <  0) ,  with  the  j(0  axis  being  the  stability  boundary.  In  the  discrete  (or 
digital)  domain  the  unit  circle,  I  z  I  <  1 ,  is  the  stability  boundary.  Discrete  roots  located 
within  the  unit  circle  indicate  a  stable  system. 


A.5  Frequency  Response  Functions 

The  frequency  response  function  (FRF)  is  another  means  of  observing  characteristics  of 
dynamic  systems.  Given  a  system,  Gp(s) ,  the  FRF  consists  of  a  magnitude  and  phase, 

often  plotted  in  Bode  form,  where  the  magnitude  is  found  by 


magnitude  =  =  |Gp  (7®)| 

=  ^{Re[Gp(;(D)]}^  +{Im[Gp(7a))]}^ 


[A.12] 


and  the  phase  is  found  by 


phase  =  ZGp(j(0)  =  tan 


'lm[GpU(0)] 

[mGpUm' 


[A.13] 
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By  plotting  the  open-loop  FRF,  Gp(s) ,  with  s  taking  on  values  of  y(0  we  are  basically 

determining  the  stability  of  the  closed-loop  system,  since  the  jto  axis  is  the  stability 

boundary.  Two  measures  of  stability  come  from  the  FRF— the  gain  and  the  phase 
margins. 

The  gain  margin  (GM)  is  that  factor  by  which  the  gain  can  be  increased  before  instability 
occurs  [51:375].  It  is  found  from  the  FRF  magnitude  plot  by  measuring  the  distance  from 
the  |^G(;a))|  curve  and  the  |^G(;Q))|  =  1  line  at  the  frequency  where  ZGiM  = -1S0° 
on  the  FRF  phase  plot.  For  stability,  GM  must  be  greater  than  one.  The  phase  margin 
(PM)  is  the  amount  by  which  the  phase  of  G(jGi)  at  |^G(7C0)|  =  1  exceeds  -180°.  For 
stability,  PM  must  be  positive. 


A.6  Sampling 

Kuo  [7.16-19]  defines  an  ideal  sampler  as  one  in  which  the  sampling  duration  is  very 
much  smaller  than  both  the  sampling  period  and  the  largest  time  constant  of  the  input 
signal.  The  output  of  this  ideal  sampler  is  written  as 


k=0 


[A.  14] 
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where  f(t)  is  the  input  signal,  T,  is  the  sample  period,  and  5  is  the  Dirac  delta  function. 
The  notation  of  an  asterisk  (*)  will  denote  a  sampled  (or  discrete)  version  of  a  continuous 
signal. 

Equation  A.  14,  the  equation  of  a  sampled  signal,  is  just  that — a  mathematical  model  of 
the  real  system.  Ogata  [33:153]  explains  that  by  assuming  a  zero-time  duration  of  the 
real  sampler  pulse,  one  can  assume  that  the  sequence  of  pulses  output  from  the  sampler 
becomes  a  sequence  of  impulses  with  strengths  equal  to  the  continuous-time  signal  at  the 
sampling  instants.  This  assumption,  then,  allows  for  the  simple  analysis  of  discrete-time 
systems.  The  assumption  is  valid  if  the  sampling  duration  is  very  small  in  comparison  to 
the  time  between  samples  and  if  some  sort  of  hold  is  connected  to  the  output  of  the 
sampler  to  drive  a  continuous  plant  with  a  sampled  signal. 

Several  properties  of  the  sampled  signal,  become  apparent.  It  is  a  unique 

function.  However,  because  it  is  unique,  that  does  not  mean  the  function  f(t)  sampled 
to  produce  /  *  (t)  is  unique.  Any  number  of  functions  may  produce  /  *  (0 .  given  they 
all  equal  /  *  (t)  at  the  sample  instants.  The  function  /  *  (t)  is  identically  equal  to  zero 
between  the  sample  times  and  equal  to  impulses  at  the  sample  times.  As  such,  it  is  a 
discrete  representation  of  a  continuous  function,  but  can  also  be  thought  of  as  a 
continuous  function.  It  is  also  an  impulse  modulator — sampling  /  *  (r)  gives  simply 
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/  *  (0  >  which  means  that  it  is  also  one  of  the  infinite  continuous  functions  we  could 
sample  to  generate 

Taking  the  Laplace  transform  of  Equation  A.  14  gives 


Jfe=0 


[A.15] 


or  alternatively, 


■^S  «=“Oo 


[A.  16] 


The  function  F*(s)  can  be  thought  of  as  the  sampled  version  of  F(s) ,  the  Laplace 
transform  of  /(t) .  Kuo  [7:22-23]  and  Ragazzini  &  Zadeh  [32]  define  several  important 
properties  of  F  *  (j) ; 

1 .  From  Equation  A.  16  it  is  easy  to  see  that  F  *  (s)  is  a  periodic  function  with 
period  joD/. 

F\s  +  pnGi^)  =  F\s)  ;  m  integer.  [A.  17] 
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Because  /  *  (t)  is  sampled  in  time,  F  *  (s)  becomes  periodic  in  frequency. 


2.  It  then  follows  that  because  F  *  (s)  is  periodic,  if  F  *  (5)  has  a  pole  or  a  zero 
at  s  =  Si,  then  F  *  (s)  has  poles  or  zeros  at  5  =  5,  +  ymco^ ;  m  integer  from 

—  to  +  «>  . 

3.  F(s)  =  Gis)R  *  (s)  implies  F*(s)  =  G*  is)R  *  (s) ,  or  mathematically: 

[G(j)F*(5)]*=G*(j)/?*(5)  .  [A.18] 

In  other  words,  because  F  *  (5)  is  already  periodic,  shifting  it  an  integral 
number  of  periods  leaves  it  unchanged,  or  sampling  a  sampled  signal  results  in 
the  same  sampled  signal  as  described  above. 

Taking  the  z-transform  of  Equation  A.  14  gives  [31:2] 

F*iz)=Z{f*(t)}  =  f^fikT,)z-'^  .  [A.19] 

ifc=0 
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Notice  the  similarity  of  Equation  A.  19  to  Equation  A.  15.  Indeed,  as  before  mentioned, 
the  Laplace  transform  of  f*  (t)  is  the  same  as  the  z-transform  of  /  *  (t)  given  that 

z  =  .  [A.20] 

Because  of  this,  the  mapping  of  F(z)  to  F  *  (5')  is  unique,  while  the  mapping  of  F(z)  to 
FCi')  is  not.  This  unique  form  of  F*(s),  then,  is  the  form  we  will  use  in  our  design 
process. 

One  important  aspect  of  sampling  comes  from  Shannon  [1]  [52].  Shannon  states  that  if  a 
signal  contains  no  frequency  content  above  a  certain  frequency,  say  ,  then  that  signal 

can  be  completely  characterized  following  sampling  only  if  the  sampling  frequency,  /^, 
is  greater  than  twice  the  greatest  frequency  contained  in  the  signal.  In  other  words, 

/.  >2/v,  .  [A.21] 

An  example  of  the  effects  of  sampling  follows  for  a  second  order  system  with  a  natural 
frequency  ©„  =  100  rad/sec  and  a  damping  ratio  ^  =  0.01.  For  purposes  of  this  example, 

a  sampling  rate  of  CO2  =  200  rad/sec  was  selected.  Notice,  however,  that  by  so  selecting 
the  sample  rate.  Shannon’s  criterion  for  sampling  was  not  met.  The  effects  of  sampling 
are  shown  well  in  the  frequency  response  function  of  Figure  A.  1 .  If  a  signal  is  sampled  at 


a  rate  lower  than  that  defined  by  Shannon,  a  phenomenon  known  as  aliasing  occurs.  The 
phenomenon  is  indeed  observable  in  Figure  A.  1 . 


Frequency  Response  Function  of  Second  Order  System 
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Figure  A.l  The  Effects  of  Sampling  on  a  Second-Order  System 
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A.7  The  Zero-Order  Hold  (ZOH) 


Following  any  digital  controller,  the  output  signal  must  be  converted  to  a  continuous 
signal  for  input  to  the  continuous  plant.  This  is  accomplished  in  the  digital-to-analog 
converter  (DAC).  The  DAC  converts  the  binary  output  of  the  controller  to  an  analog 
voltage.  Since  the  output  consists  of  values  at  discrete  times,  a  method  for  “holding”  the 
analog  value  from  one  discrete  time  to  the  next  must  be  incorporated.  The  most  common 
means  of  accomplishing  this  is  through  the  Zero-Order-Hold  (ZOH). 

The  ideal  ZOH  simply  holds  the  converted  output  voltage  constant  from  one  time  step  to 
the  next.  In  effect,  then,  the  ZOH  converts  the  series  of  binary  impulses  from  the 
controller  into  a  series  of  rectangular  voltage  pulses.  These  pulses  have  a  height  equal  to 
the  value  of  the  binary  output  and  a  width  equal  to  the  sample  period. 

Assuming  the  binary  impulse  input  to  the  ZOH  is  a  unit  impulse,  the  impulse  response  for 
the  ZOH  is  given  by 

zoh{t)  =  u{t)-  u{t  -  r, )  ,  [A.22] 


where  u{t)  is  the  unit  step  function  and  is  the  sample  period.  Taking  the  Laplace 
transform  of  Equation  A.22  gives 
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[A.23] 


1-e  * 

ZOH(s)  = - 

s 

Letting  s  =  jG3,  and  after  some  mathematical  manipulation  [1:165-166]  or  [7:28],  we 
have 


ZOH(j(0)  =  sine 


(ofTA  -J 


V  2  7 


[A.24] 


where  the  sine  function  is  defined  by 


sinc(jc)  = 


sin(x) 


[A.25] 


Examination  of  Equation  A.24  shows  the  ZOH  effectively  introduces  a  phase  shift  of 

cor  f  (aT  ^ 

— -  and  a  change  in  magnitude  equal  to  7^  sine  — -  .  Plotting  the  ZOH,  as  in 

2  V  ^  j 

Figure  A.2,  also  shows  this. 
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The  Zero-Order  Hold 


Figure  A.2  The  Zero-Order  Hold 


To  show  the  effect  of  the  ZOH  applied  to  a  signal,  the  lightly  damped  system  from  above 
(cOn  =  100  rad/sec,  ^  =  0.01)  is  pre-multiplied  by  a  ZOH.  The  sampling  frequency  is  set  to 
twice  the  natural  frequency.  The  results  are  shown  below  in  Figure  A.3. 

One  thing  to  notice  in  Figure  A.3  is  the  wrapping  of  the  phase.  This  will  be  the  case  in 
many  of  the  plots  to  follow,  as  it  aids  in  a  scale  for  the  phase  plots  which  allows  closer 
examination  of  the  phase  relationships. 
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TheZOH  Effects 


Figure  A.3  Zero-Order  Hold  Effects  on  Second-Order  System 


A.8  Appendix  A  Summary 


Here,  the  equations  for  the  Laplace-  and  Z-transforms  and  the  natural  logarithm  of  a 
complex  number  were  introduced  and  reviewed.  The  concepts  of  root  loci  and  frequency 
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response  function  plots  were  discussed.  The  definition  of  an  ideal  sampler  and  its 
mathematical  model  was  given,  and  the  effects  of  sampling  were  examined.  Finally,  the 
zero-order  hold  and  its  effect  on  a  continuous  system  were  reviewed.  These  results  are 
used  extensively  in  the  development  of  this  research. 
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Appendix  B.  Second-Order  Model  Implementation 
Using  Matlab  (m-file) 


%  mlOlst.m  —  m-file  to  examine  design  of  2nd-order  discrete-to- 
%  continuous  system  of  the  form: 

% 

%  1 

%  Gp(s)  =  -  , 

%  (s  +  pi) (s  +  p2) 

% 

%  with  proportional/ integral  control 
% 

%  z  Kp{z  -  1)  +  Ki*z  K*z  -  Kp 

%  Gc(z)  =  Kp  +  Ki* -  = - = - , 

%  z-1  z-1  z-1 

% 

%  where 
% 

%  K  =  Kp  +  Ki  . 

% 

%  This  will  plot  z-plane  and  s* -plane  root  loci,  bode  plot,  impulse  and 
%  step  responses,  and  response  to  sine-wave  inputs  at  two  different 
%  frequencies , . * 

% 

%  David  S.  Hansen,  20  April  2000 
clear  all 

%  User  defined  parameters 
wn  =  1; 
zeta  =  0.1; 

wd  =  wn*sqrt(l  -  zeta^2) ; 
pi  =  wn*{zeta  -  sqrt(zeta^2  -  1)); 
p2  =  conj(pl); 

Kp  =  0.5; 

Ki  =  0; 

K  =  Kp  +  Ki; 

FR  =  2/4; 

ws  =  FR*wd; 
frequency 
Ts  =  (2*pi)/ws; 
max time  =  150; 
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%  natural  frequency 
%  damping  ratio 
%  damped  natural  frequency 
%  plant  poles 

%  proportional  gain 
%  Integral  gain 

%  frequency  ratio  (equal  sampling 
%  freq.  over  damped  nat.  freq.) 

%  sample 

%  sample  time 

%  Maximum  time  of  simulation 


istepsl  =  40; 


%  Number  of  intermediate  steps 
%  between  each  sample  times 

%  size  of  step  for  root  locus  calcs. 


isteps2  =  40; 
stepsize  =  .05; 


%  ***  Plot  roots  and  poles  in  both  z-  and  s*-planes  for  defined 

%  Find  zeros  and  poles  of  H(s)  =  ZOH(s)*Gp(s) 
polelc  =  -pi;  %  from  Gp 

pole2c  =  -p2; 

polelc  =0;  %  from  ZOH 

%  Using  partial  fraction  expansion,  Gp{s)/s  is 
% 

%  Gp(s)  aO  al  a2 

%  s  s  s+pl  s+p2 

% 

%  where  aO,  al,  and  a2  are: 
aO  =  l/(pl*p2); 
al  =  l/(pl*(pl-p2)); 
a2  =  -l/{p2*(pl-p2)); 

%  Taking  z-transform  of  H(s)  gives 
% 

%  cl*z  +  c2 

%  H (z)  =  - 

%  (z  +  bl) (z  +  b2) 

% 

%  where  bl,  b2 ,  cl  and  c2  are: 
bl  =  -exp (-pl*Ts) ; 
b2  =  -exp {-p2*Ts) ; 

cl  =  real {wn^2* (aO  -  al*bl  -  a2*b2)); 

c2  =  real (wn^2* (a0*bl*b2  -  al*b2  -  a2*bl) ) ; 

%  Continuous  zeros  and  poles  "in"  z-domain: 

poleld  =  exp (polelc*Ts) ; %  plant  poles  from  both  H(z)  and  H(s)  -> 
pole2d  =  exp (pole2c*Ts) ; 
zerold  =  0; 
zero2d  =  -c2/cl; 
poleld  =  1; 
zerold  =  1; 

%  Recall, 

% 

%  (Gcnum)  (Hden)  (K*z  -  Kp)  (z  +  bl)  (z 

%  F{z)  - -  =  - 

%  (Gcnum)  (Hnum)  +  (Gcden)  (Hden)  z^'l  +  dl*z''2  +  d2*z  + 

% 

%  where 

dl  =  K*cl  +  bl  +  b2  -  1; 

d2  =  K*c2  -  Kp*cl  +  bl*b2  -  bl  -  b2; 

dl  =  -Kp*c2  -  bl*b2; 


%  plant  zeros  (2)  from  H(s)  ->  z-plane 
%  plant  zero  from  H(z) 

%  ZOH  pole  and  zero 


Kp  *** 


z -plane 


+  b2) 
d3 
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poles  =  roots ( [1  dl  d2  d3] ) ; 


%  Now  find  values  for  specified  gains... 
zerolz  =  -Kp/K;  %  Controller  zero 

zero2z  =  -bl;  %  Hden  zeros 

zero3z  =  -b2; 

polelz  =  poles (1);  %  F(z)  poles 

pole2z  =  poles (2 ) ; 
pole3z  =  poles (3 ); 

%  Now  transform  to  s* -plane  and  repeat 

zls  =  l/Ts*log (zerolz) ;  %  Gc(z)  zero 

z2s  =  l/Ts*log (zero2z) ;  %  Zeros  from  Hden{z) 

z3s  =  l/Ts*log{zero3z) ; 

z4s  =  l/Ts*log(zero2d) ;  %  H{z)zero 

pis  =  l/Ts*log(polelz) ;  %  F(z)  poles 

p2s  =  l/Ts*log(pole2z) ; 
p3s  =  l/Ts*log(pole3z) ; 
for  loop  =  1:5 

zerols(loop)  =  zls  +  (loop  -  3)*j*ws;  %  combined  zeros 

zero2s(loop)  =  z2s  +  (loop  -  3)*j*ws; 

zero3s(loop)  =  z3s  +  (loop  -  3)*j*ws; 

zero4s(loop)  =  z4s  +  (loop  -  3)*j*ws;  %  "discrete"  zero 

polels(loop)  =  pis  +  (loop  -  3)*j*ws;  %  combined  poles 

pole2s(loop)  =  p2s  +  (loop  -  3)*j*ws; 

pole3s(loop)  =  p3s  +  (loop  -  3)*j*ws; 

zerolc(loop)  =  j*(loop  -  3)*ws;  %  zeros  from  ZOH(s) 

end 

%  Plot  results 
figured);  clf; 

%  z-domain  root  locus 
subplot (221) ;  hold  on; 

plot (real (poleld) , imag (poleld) ,'*') ;  %  "continuous"  poles 

plot (real (pole2d) , imag (pole2d) ,'*'); 

plot (real (zerold) , imag (zerold) , 's') ;  %  "continuous"  zeros 

plot (real (zero2d) , imag (zero2d) d' ) ;  %  discrete  zero 

plot (real (zero3d) , imag (zero3d) ,  "^' ) ;  %  zoh  zero 

plot (real (pole3d) , imag (pole3d) ,'+') ;  %  zoh  pole 

plot (real (polelz) , imag (polelz) , 'x' ) ;  %  combined  poles 

plot (real (pole2z) , imag (pole2z) , 'x' ) ; 

plot (real (pole3z) , imag (pole3z) ,'x'); 

plot (real (zerolz) , imag (zerolz) o' ) ;  %  combined  zeros 

plot ( real ( zero2  z ) , imag ( zero2  z ) , ' o ' ) ; 
plot ( real ( zero3  z ) , imag ( zero3  z ) , ' o ' ) ; 

hold  off;  zgrid;  axis  square;  axis ([-1.1  1.1  -1.1  1.1]); 
phrase  =  ['Z-Plane  Root  Locus;  FR  =  ' , num2str (FR) , 

';  Kp  =  ' ,num2str (Kp) , ' ;  Ki  =  ' , num2str (Ki) ] ; 
title (phrase) 

%  s'^-domain  root  locus 

vectl  =  [-1.5  1];  %  vectors  for  plotting  axes  and 

strip  boundaries 

vect2  =  [-5*ws/2  -5*ws/2] ; 

vect3  =  [-3*ws/2  -3*ws/2] ; 
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vect4  =  [-ws/2  -ws/2]; 
vectS  =  [0  0] ; 
vect6  =  [-2.5  2.5]; 
subplot (223) ;  hold  on; 

plot(real(polelc) ,imag(polelc) , ;  %  plant  poles 

plot (real (pole2c) , imag(pole2c) 

plot  (real  (poleSc)  ,iinag(pole3c)  ;  %  zoh  pole 

for  loop  =  1:5 

plot (real (polels (loop) ), imag (polels (loop) ), 'x' ) ;  %  combined  poles 

plot (real (pole2s (loop) ) , imag (pole2s (loop) ) ,’x') 
plot (real (pole3s (loop) ) , imag {pole3s (loop) ) , 'x' ) 
plot (real (zerols (loop) ) , imag (zerols (loop) ) , 'o' ) 
plot (real (zero2s (loop) ) , imag (zero2s (loop) ) , 'o' ) 
plot (real (zero3s (loop) ) , imag (zero3s (loop) ) , 'o') 
plot (real (zero4s (loop) ) , imag (zero4s (loop) ) , 'd') 
plot ( real (zerolc( loop) ) , imag (zerolc (loop) ),"'') 
end 

plot ( vectl , vect5 ,':');  plot ( vect5 , vect6 ,':'); 
plot(vectl,vect2,  '  —  plot  (vectl, vect3  ,  '--'); 
plot  (vectl ,  vect4 ,'  —  ');  plot  (vectl ,  -vect2 ,  '  — ' )  ; 
plot (vectl , -vect3 , ' — ' ) ;  plot (vectl , -vect4 , ' — ' ) ; 
hold  off;  axis ([-1.5  1  -2.5  2.5]); 

phrase  =  ['S*-Plane  Root  Locus;  FR  =  ' ,num2str (FR) , 

';  Kp  =  ' ,num2str(Kp) , ' ;  Ki  =  ' ,num2str (Ki) ] ; 
title (phrase) 


%  ***  Plot  System  Bode  Plot  for  defined  Ki  *** 

nximpoints  =  2000;  %  Number  of  frequency  points  for  FRF 

num  =  wn''2; 

den  =  [1  pl+p2  pl*p2] ; 

Gps  =  tf (num, den);  %  Continuous  plant 

Hz  =  c2d(Gps,Ts, 'zoh') ;  %  H*(z):  Sampled  Plant  with  ZOH 

w  —  linspace ( 0 . 001 , 50 , numpoints ) ;  %  Linearly  spaced  frequency  vector 

normw  =  w./wd;  %  Normalized  frequency  vector  for  plotting 

%  Porportional/ Integral  Controller 
Gcznum  =  [K  -Kp] ; 

Gczden  =  [1  -1] ; 

Gcz  =  tf(Gcznum,Gczden,Ts); 

[magGcz,phsGcz]  =  bode(Gcz,w) ; 
magGcz  =  reshape (magGcz, 1 , numpoints) ; 
phsGcz  =  reshape (phsGcz , 1 , numpoints ) ; 

phsGcz  =  phsGcz*pi/180 ;  %  Convert  phase  from  degrees  to  radians. 

%  Continuous  System 
[magGps,phsGps]  =  bode (Gps, w); 
magGps  =  reshape (magGps, 1 , numpoints ) ; 
phsGps  =  reshape (phsGps,!, numpoints ) ; 

phsGps  =  phsGps *pi/l 80;  %  Convert  phase  from  degrees  to  radians. 

%  Discrete  System 

[magHz,phsHz]  =  bode(Hz,w); 

magHz  =  reshape (magHz, 1, numpoints); 


%  combined  zeros 


%  "discrete"  zeros 
%  zoh  zeros 


phsHz  =  reshape (phsHz, l,numpoints) ; 

phsHz  =  phsHz*pi/180;  %  Convert  phase  from  degrees  to  radians. 

%  Calculate  magnitude  and  phase  at  discrete  frequency  points 
for  loop  =  l:numpoints 

ZOH  =  (l-exp(-j*w(loop) *Ts) ) / ( j*w(loop) *Ts) ; 

Gpsl  =  magGps (loop) *exp (j *phsGps (loop) ) ; 

Gczl  =  magGcz (loop) *exp(j*phsGcz (loop) ) ; 

Hpz  =  magHz (loop) * exp (j*phsHz (loop) ) ; 
numcl  =  Gczl*ZOH*Gpsl ; 
dencl  =  1  +  (Gczl*Hpz) ; 
cltf(loop)  =  numcl/dencl; 

end 

magFRF  =  20*logl0 (abs (cltf ) ) ;  %  Magnitude  (in  dB)  of  FRF 

phsFRF  =  angle(cltf ) *180/pi;  %  Phase  (in  degrees)  of  FRF 

a  =  max  (magFRF)  ; 

plotmax  =  a  +  5; 

plotmin  =  plotmax  -  100; 

subplot (222) ; 

plot (no rmw, magFRF) ;  grid; 

phrase  =  ['Frequency  Ratio  Ws/Wd  =  ' ,num2str (FR) , 

';  Kp  =  ' , num2str (Kp) , ' ;  Ki  =  ' , num2str (Ki) ] ; 
title (phrase) ; 

xlabel ( 'Normalized  Frequency  (W/Wd) ' ) ; 

ylabel ( 'Magnitude  (dB) ' ) ; 

axis ( [0  4.5  plotmin  plotmax] ) ; 

subplot (224) ; 

plot (normw,phsFRF) ;  grid; 

axis( [0  4.5  -190  190]) ; 

xlabel ( 'Normalized  Frequency  (W/Wd) ' ) ; 

ylabel ( ' Phase  (deg) ' ) ; 


%  ***  Plot  Impulse  Response  for  defined  Ki  *** 

%  ====  =  =  ===:  =  =  =:  =  =  =  ===:==  =  =  =  ===  =  =  =  =  =  =  =  =  =  =  ===  ===:=  = 

Tn  =  Ts/istepsl;  %  Intermediate  time  steps 

[F,G,H,J]  =  tf 2ss (num, den) ; 

[phi, gamma]  =  c2d(F,G,Tn);  %  Continuous  approximation  in 

%  discrete  domain 

%  Define  initial  conditions 
X  =  [  0 ; 0 ] ;  uthen  =  0 ;  errorthen  =  0 ; 
tcon  =  0;  ycon  =  0; 
tdis  =  0;  ydis(l)  =  0; 

%  Step  through  calculations  for  "maxtime"  seconds 
for  countl  =  0 : round (maxtime/Ts) ; 

%  Define  unit  impulse: 
if  countl  ==  0 
r  =  1; 
else 

r  =  0; 

end 

errornow  =  r  -  ydis (countl  +  1); 


257 


u  =  uthen  +  K*errornow  -  Kp * err or then; 


%  Control  effort 


%  Calculate  continuous  output 
for  count2  =  liistepsl; 

tcon  =  [tcon; (countl*istepsl  +  count2)*Tn]; 

X  =  phi*x  +  gainma*u; 
ycon  =  [ycon; (H*x  +  J*u) ] ; 
end 

index  =  istepsl* (countl  +  1)  +  1; 
tdis  =  [ tdis ; tcon ( index) ] ; 
ydis  =  [ydis ;ycon (index) ] ; 
uthen  =  u; 

err or then  =  errornow; 

end 

%  Plot  Results 
a  =  max  (ycon)  ; 
b  =  min (ycon) ; 
plotmax  =  a  +  0.2*a; 
plotmin  =  b  +  0.2*b; 
figure (2);  clf; 
subplot (221)  ; 

plot {tcon. /Ts, ycon, tdis . /Ts , ydis , 'r-- tdis . /Ts, ydis , 'r*') ;  grid; 

xlabel ( 'Normalized  Time  (t/Ts) ' ) ;  ylabel ( 'Response' ) ; 

axis([0  maxtime/Ts  plotmin  plotmax]); 

phrase  =  ['Impulse  Response  at  Ws/Wd  =  ' ,num2str (FR) , 

';  Kp  =  ' ,num2str (Kp) , ' ;  Ki  =  ' ,num2str (Ki) ] ; 
title (phrase) ; 


%  ***  Plot  Step  Response  for  defined  Kp  *** 

clear  x  uthen  tcon  ycon  tdis  ydis 
%  Define  initial  conditions 
X  =  [0;0];  uthen  =  0;  errorthen  =  0; 
tcon  =  0;  ycon  =  0; 
tdis  =  0;  ydis(l)  =  0; 

2:  =  1;  %  Unit  step 

%  Step  through  calculations  for  "maxtime"  seconds 
for  countl  =  0 : round (maxtime/Ts ) ; 
errornow  =  r  -  ydis (countl  +  1) ; 

u  =  uthen  +  K*errornow  -  Kp*errorthen;  %  Control  effort 

%  Calculate  continuous  output 
for  count2  =  1 : istepsl; 

tcon  =  [tcon; (countl*istepsl  +  count2)*Tn]; 
x  =  phi*x  +  gamma *u; 
ycon  =  [ycon; (H*x  +  J*u) ] ; 
end 

index  =  istepsl* (countl  +  1)  +  1; 
tdis  =  [tdis ; tcon (index) ] ; 
ydis  =  [ydis ; ycon ( index) ] ; 
uthen  =  u; 

errorthen  =  errornow; 
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end 


%  Plot  Results 
a  =  inax{ycon)  ; 
b  =  rain{ycon)  ; 
plotmax  =  a  +  0.2*a; 
plotmin  =  b  +  0.2*b; 
subplot (223 ) ; 

plot (tcon. /Ts,ycon, , tdis . /Ts,ydis , 'r — ' , tdis . /Ts ,ydis , 'r*') ;  grid; 

xlabel ( 'Normalized  Time  (t/Ts) ' ) ;  ylabel ( 'Response' ) ; 

axis  ( [0  maxtime/Ts  plotmin  plotmax] )  ; 

phrase  =  ['Step  Response  at  Ws/Wd  =  ' , num2str (FR) 

Kp  =  ' ,num2str (Kp) , ' ;  Ki  =  ' , num2str (Ki) ] ; 
title (phrase) ; 


%  ***  Plot  time  response  for  sine-wave  input  *** 


clear  x  uthen  tcon  ycon  tdis  ydis 

ww  =  0.05;  %  Input  sine  wave  frequency 

maxtime  =  200;  %  Maximum  time  for  simulation  (in 

seconds) 

Tn  =  Ts/isteps2;  %  Intermediate  time  steps 

%  2nd-Order  State-space  System  paramters 

[phi, gamma]  =  c2d(F,G,Tn);  %  Continuous  approximation  in 

%  discrete  domain 

%  Define  initial  conditions 

X  =  [  0 ; 0 ] ;  uthen  =  0 ;  errothen  =  0 ; 

tcon  =  0;  ycon  =  0; 

tdis(l)  =  0;  ydis(l)  =  0; 

index  =  0;  indexl  =  0; 


%  Step  through  calculations  until  "maxtime"  seconds 
for  countl  =  0 : round (maxtime/Ts ) ; 

%  Define  sinusoidal  input  and  control  as  a  result... 

time  =  countl*Ts; 

r  =  sin(ww*time)  ; 

errornow  =  r  -  ydis (countl  +  1) ; 

u  =  uthen  +  K*errornow  -  Kp*errorthen; 

%  Calculate  continuous  output 
for  count2  =  l:isteps2 
indexl  =  indexl  +  1 ; 

t (indexl)  =  (countl*isteps2  +  count2)*Tn; 
y(indexl)  =  sin (ww*t ( indexl) ) ; 

tcon  =  [tcon; (countl*isteps2  +  count2)*Tn]; 

X  =  phi*x  +  gamma*u; 

ycon  =  [ycon; (H*x  +  J*u) ] ; 

end 

index  =  isteps2* (countl  +  1)  +  1; 
tdis  =  [tdis ; tcon (index) ] ; 
ydis  =  [ydis ;ycon( index) ] ; 
uthen  =  u; 
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errothen  =  errornow; 


end 

%  Plot  Time  Response 
a  =  max(ycon) ; 
b  =  min  (ycon)  ; 
if  a  >  1 

plotmax  =  a  +  0.2*a; 
else 

plotmax  =  1; 

end 
if  b  < 

plotmin  =  b  +  0.2*b; 
else 

plotmin  =  -1; 

end 

subplot (222) ;  grid; 

plot ( tcon. /Ts, ycon, ' - ' , tdis . /Ts,ydis, '*',t./Ts,y, ' — ') ; 
xlabel ( 'Normalized  Time  (t/Ts) ' ) ;  ylabel ( 'Response' ) ; 
phrase  =  ['Response  to  sin( ' ,num2str (ww) , ' *t)  input;  FR  =  ' , 

num2str(FR) , Kp  =  ' ,num2str (Kp) , ' ;  Ki  =  ' ,num2str (Ki) ] ; 
title (phrase) ; 

axis([0  maxtime/Ts  plotmin  plotmax]); 
clear  x  uthen  tcon  ycon  tdis  ydis 

ww  -  0.09;  %  Input  sine  wave  frequency 

%  Define  initial  conditions 

X  =  [  0 ; 0 ] ;  uthen  =  0 ;  errorthen  =  0 ; 

tcon  =  0;  ycon  =  0; 

tdis  =  0;  ydis(l)  =  0; 

index  =  0;  indexl  =  0; 

%  Step  through  calculations  until  "maxtime"  seconds 
for  countl  =  0 : round (maxtime/Ts ) ; 

%  Define  sinusoidal  input  and  control  as  a  result... 

time  =  countl*Ts; 

r  =  sin(ww*time) ; 

errornow  =  r  -  ydis (countl  +  1) ; 

u  =  uthen  +  K*errornow  -  Kp*errorthen; 

%  Calculate  continuous  output 
for  count2  =  l;isteps2 
indexl  =  indexl  +  1 ; 

t(indexl)  =  (countl*isteps2  +  count2)*Tn; 

y(indexl)  =  sin (ww*t (indexl) ) ; 

tcon  =  [tcon; (countl*isteps2  +  count2)*Tn]; 

X  =  phi*x  +  gamma*u; 
ycon  =  [ycon; (H*x  +  J*u) ] ; 

end 

index  =  isteps2* (countl  +1)  +  1; 
tdis  =  [tdis; tcon (index) ] ; 
ydis  =  [ydis ; ycon ( index) ] ; 
uthen  =  u; 

errorthen  =  errornow; 

end 
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%  Plot  Time  Response 
a  =  max(ycon); 
b  =  min  (ycon)  ; 
if  a  >  1 

plotmax  =  a  +  0.2*a; 
else 

plotmax  =1; 

end 

if  b  <  -1 

plotmin  =  b  +  0.2*b; 
else 

plotmin  =  -1; 
end 

subplot (224) ;  grid; 

plot  (toon.  /Ts ,ycon,  '  - '  ,  tdis .  /Ts,ydis,  '  * '  ,  t .  /Ts  ,y, 
xlabel ( 'Normalized  Time  (t/Ts)');  ylabel {' Response ') ; 
phrase  =  ['Response  to  sin{ '  ,niim2str  (ww)  ,  '  *t)  input;  FR  =  '  , 

num2str (FR) , ' ;  Kp  =  ' , num2str (Kp) , ' ;  Ki  =  ' ,num2str (Ki) ] ; 
title (phrase) ; 

axis ( [0  maxtime/Ts  plotmin  plotmax] ) ; 
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Appendix  C.  Fourth-Order  Model  Implementation 
Using  Matlab  (m-file) 


%  fourth. in  m-file  to  examine  design  of  a  4th- order  system  of  the 
%  form: 

% 

%  Kg 

%  Gp(s)  =  - 

%  (s+pl) {s+p2) (s+p3) (s+p4) 

% 

% 

%  =  K*wnl^2*wn2^2 

%  -  ^ 

%  (s^2+2*zetal*wnl*s+wnl'^2 )  (s'^2+2*zeta2*wn2*s+wn2^2 ) 

% 

%  with  PID  control 
% 

%  z-1  z 

%  Gc(z)  =  Kp  +  Kd* -  +  Ki* - , 

%  z  z-1 

% 

%  This  will  plot  z-plane  and  s*-plane  root  loci,  bode  plots,  and  impulse 
%  and  step  responses 
% 

%  David  S.  Hansen,  17  May  2000 
clear  all; 

%  User-defined  parameters 
Kp  =  0.7; 

Kd  =  1.4; 

Ki  =  0.6; 

Kl  =  Kp  +  Kd  +  Ki; 

K2  =  -Kp  -  2*Kd; 

FR  =  1; 

K  =  1.0; 
maxtime  =  15 ; 
istepsl  =  100; 

%  Non-aliased  mode 
wnl  =  1.0; 
zetal  =  0.5; 

pi  =  wnl* (zetal-sqrt (zetal^2-l) ) ; 


%  Proportional  gain 
%  Derivative  gain 
%  Integral  gain 


%  Frequency  ratio  (ws/wd) 

%  System  gain 

%  maximum  time  for  simulations 
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p2  =  conj(pl); 

wdl  =  wnl*sgrt (l-zetal^2) ; 

%  Aliased  mode 
wn2  =  10; 
zeta2  =  0.01; 

p3  =  wn2* (zeta2-sqrt (zeta2^2-l) ) ; 

p4  =  con j (p3 ) ; 

wd2  =  wn2*sqrt(l-zeta2^2) ; 

Kg  =  K*wnl^2*vm2^2; 

%  Sampling  Parameters 
ws  =  FR*wd2 ; 

Ts  =  (2  *pi ) /ws ; 


%  ***  Plot  Root  Loci  *** 


%  Continuous  Plant  Poles 
polelc  =  -pi; 
pole2c  =  -p2 ; 
pole3c  =  -p3; 
pole4c  =  -p4; 

%  Expand  Gp(s)/s  using  partial  fractions 
% 

%  Gp{s)  aO  al  a2  a3  a4 

^  ~  + - + - + -  (all  multiplied  by  Kg) 

%  s  s  s+pl  s+p2  s+p3  s+p4 

% 

%  where  the  constants  are  defined  below. . . 

p5  =  (Pl*p2)+(pl*p3)+(pl*p4)+(p2*p3)+(p2*p4)+(p3*p4) ; 
p6  =  (p2*p3) + (p2*p4) + {p3*p4) ; 
p7  =  (pl*p3)+(pl*p4)+(p3*p4) ; 
p8  =  (pl*p2)+(pl*p4)+{p2*p4) ; 
p9  =  (pl*p2)+(pl*p3)+(p2*p3) ; 

plO  =  (Pl*p2*p3)+(pl*p2*p4)+(pl*p3*p4)+(p2*p3*p4) ; 

pll  =  p2*p3*p4; 

pl2  =  pl*p3*p4; 

pl3  =  pl*p2*p4; 

pl4  =  pl*p2*p3; 

do  =  1/ (pl*p2*p3*p4) ; 
dl  =  -d0*pl/ (pl-p2) ; 
d2  =  (p3-pl)/(pl-p2) ; 
d3  =  (p4-pl) / (pl-p2) ; 
d5  =  -dO-dl; 
d6  =  l+d2; 
d7  =  l+d3; 

d8  =  (d0*p5)+(d5*p6)+{dl*p7) ; 
d9  =  (d2*p7)+p8- (d6*p6) ; 
dlO  =  (d3*p7)+p9-(d7*p6) ; 
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dll  =  -d8/d9; 
dl2  =  -dl0/d9; 
dl3  =  d5-{d6*dll); 
dl4  =  -d7-{d6*dl2)  ; 
dl5  =  dl+(d2*dll); 
dl6  =  d3+(d2*dl2); 

%  ...aahh,  here  they  are 
aO  =  dO; 

a4  =  -( (d0*pl0)+(dl3*pll)+(dl5*pl2)+{dll*pl3) ) 
a4  =  a4/ ( {dl4*pll)+(dl6*pl2)+(dl2*pl3)+pl4) ; 
a3  =  dll+(dl2*a4) ; 
a2  =  dl5+(dl6*a4) ; 
al  =  dl3+(dl4*a4) ; 

%  Taking  the  z~ transform  of  H(s) 

% 

%  Kg*(f3*z^3  +  f2*z'^2  +  fl*z  +  fO) 

%  H(z)  =  - 

%  (z  +  bl)  (z  +  b2)  (z  +  b3)  (z  +  b4) 

% 

%  where 

bl  =  -exp (-pl*Ts) ; 
b2  =  -exp (-p2*Ts) ; 
b3  =  -exp (-p3*Ts) ; 
b4  =  -exp (-p4*Ts) ; 
b5  =  bl+b2+b3+b4; 

b6  =  {bl*b2)+(bl*b3)+(bl*b4)+(b2*b3)+(b2*b4)+{b3*b4) ; 
b7  =  {bl*b2*b3)+{bl*b2*b4)+{bl*b3*b4)+(b2*b3*b4) ; 
b8  =  bl*b2*b3*b4; 


%  and 

fO  =  (a0*b8) - (al*b2*b3*b4) - (a2*bl*b3 *b4) - (a3*bl*b2*b4) - (a4*bl*b2*b3 ) ; 
fl  =  a0*b7; 

fl  =  fl+(al*((b2*b3*b4)-(b2*b3)-(b2*b4)-(b3*b4))); 
fl  =  fl+(a2*( (bl*b3*b4)-(bl*b3)-(bl*b4)-(b3*b4) ) ) ; 
fl  =  fl+(a3*( (bl*b2*b4)-(bl*b2)-(bl*b4)-(b2*b4) ) ) ; 
fl  =  fl+{a4*{ (bl*b2*b3)-(bl*b2)-(bl*b3)-{b2*b3) ) ) ; 
f2  =  a0*b6; 

f2  =  f2+{al*{ (b2*b3)+(b2*b4)+(b3*b4)-b2-b3-b4) ) ; 
f2  =  f2+{a2*( (bl*b3)+(bl*b4)+(b3*b4)-bl-b3-b4) ) ; 
f2  =  f2+{a3*{ (bl*b2)+{bl*b4)+(b2*b4)-bl-b2-b4) ) ; 
f2  =  f2+{a4*{ (bl*b2)+(bl*b3)+(b2*b3)-bl-b2-b3) ) ; 
f3  =  a0-(al*bl)-(a2*b2)-(a3*b3)-(a4*b4) ; 
fO  =  real{fO) ; 
fl  =  real(fl) ; 
f2  =  real(f2) ; 
f3  =  real ( f3 ) ; 

%  This  causes  the  mixed  term  F(Z)  to  be 
% 

%  (Kl*z^2  +  K2*z  +  Kd) (z  +  bl) (z  +  b2) (z  +  b3) (z  +  b4) 

%  F(z)  =  - 

%  Kg* (Kl*z"2+K2*z+Kd) ( f3*z^3+f2*z^2+f l*z+f 0) 

%  +  z(z-l) (z+bl) (z+b2) (z+b3) {z+b4) 
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% 

%  giving  zeros  from  numerator  and  poles  from  denominator. 

%  Zeros  of  F(z) 
zerold  =  -bl; 
zero2d  =  -b2 ; 
zero3d  =  -bS ; 
zero4d  =  -b4; 

zeros  =  roots {[Kl  K2  Kd] ) ; 
zeroSd  =  zeros (1); 
zero6d  =  zeros (2 ); 

%  Poles  of  F(z) 

gl  =  Kl*Kg*f3+b5-l; 

g2  =  Kg*{Kl*f2+K2*f3)+b6-b5; 

g3  =  Kg* (Kl*fl+K2*f2+Kd*f3)+b7-b6; 

g4  =  Kg* (Kl*f0+K2*fl+Kd*f2)+b8~b7; 

g5  =  Kg*{K2*fO+Kd*fl)-b8; 

g6  =  Kg*Kd*fO; 

poles  =  roots ([1  gl  g2  g3  g4  g5  g6]); 

poleld  =  poles (1); 

pole2d  =  poles{2); 

pole3d  =  poles (3 ); 

pole4d  =  poles (4); 

poleSd  =  poles (5); 

pole6d  =  poles (6 ); 

%  Location  of  continuous  zeros  and  poles  if  they  were  "in"  the  z-plane 
zerolcd  =0;  %  3 -poles  located  at  negative 

infinity  in  s-plane. , . 

polelcd  =  exp(polelc*Ts) ;  %  Plant  poles 

pole2cd  =  exp {pole2c*Ts) ; 
pole3cd  =  exp (pole3c*Ts) ; 
pole4cd  =  exp (pole4c*Ts) ; 

%  Zero-order  hold  zero  and  pole 
zero7d  =  1; 
pole7d  =  1; 

%  Plot  z-plane  root  locus 
figured);  clf; 
subplot (211) ;  hold  on; 
plot (real (zerolcd) , imag (zerolcd) , 'd' ) ; 
plot (real (polelcd) , imag (polelcd) ,'+'); 
plot (real (pole2cd) , imag (pole2cd) ,'+'); 
plot (real (polelcd) , imag (polelcd) ,'+'); 
plot (real (pole4cd) , imag (pole4cd) ,'+'); 
plot (real (zerold) , imag (zerold) , 'o' ) ; 
plot (real (zero2d) , imag(zero2d) , 'o' ) ; 
plot (real (zero3d) , imag (zerold) , 'o') ; 
plot (real (zero4d) , imag (zero4d) , 'o') ; 
plot (real (zeroSd) , imag (zeroSd) , 'o' ) ; 
plot (real (zero6d) , imag (zero6d) , 'o') ; 
plot (real (poleld) , imag (poleld) , 'x' ) ; 
plot (real (pole2d) , imag (pole2d) , 'x') ; 


%  Gp(s)  zero 
%  Gp(s)  poles 


%  F(z)  zeros 


%  F(z)  poles 
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plot (real (pole3d) , imag(pole3d) , 'x' ) ; 
plot (real (pole4d) , imag{pole4d) , 'x' ) ; 
plot (real (poleSd) , imag(pole5d) , 'x' ) ; 
plot (real (pole6d) , imag(pole6d) , 'x' ) ; 

plot (real (zero7d) , imag(zero7d) ^ ;  %  ZOH  zero  and  pole 

plot (real (pole7d) , imag (pole7d) 

hold  off;  zgrid;  axis  square;  axis([”l  1  -1  1] ) ; 
phrase  =  ['Z-Plane  Root  Locus,  FR  =  '  ,nuin2str (FR)  ]  ; 
title (phrase) ; 

phrase  =  ['Kp  =  ' , num2str (Kp) , ' ;  Ki  =  ' , num2str (Ki) , ' ; 

Kd  =  ' ,num2str (Kd) ] ; 
xlabel (phrase) ; 


%  Now  convert  to  and  plot  in  the  s* -plane. . . 

zldc  =  l/Ts*log (zerold) ; 

z2dc  =  l/Ts*log(zero2d) ; 

z3dc  =  l/Ts*log(zero3d) ; 

z4dc  =  l/Ts*log{zero4d) ; 

z5dc  =  l/Ts*log { zeroSd) ; 

z6dc  =  l/Ts*log{zero6d) ; 

z7dc  =  l/Ts*log (zero7d) ; 

pldc  =  l/Ts*log (poleld) ; 

p2dc  =  l/Ts*log (pole2d) ; 

p3dc  =  l/Ts*log (pole3d) ; 

p4dc  =  l/Ts*log (pole4d) ; 

p5dc  =  l/Ts*log (poleSd) ; 

p6dc  =  l/Ts*log(pole6d) ; 

for  loop  =  1:4, 


zeroldc (loop) 

= 

zldc 

+ 

(loop  - 

2) *j*ws; 

zero2dc (loop) 

z2dc 

+ 

(loop  - 

2) *j*ws; 

zero3dc (loop) 

z3dc 

+ 

(loop  - 

2) *j*ws; 

zero4dc (loop) 

= 

z4dc 

(loop  - 

2)  *j*ws; 

zeroSdc (loop) 

= 

z5dc 

+ 

(loop  “ 

2) *j*ws; 

zero6dc (loop) 

= 

z6dc 

+ 

(loop  - 

2) * j  *ws; 

zero7dc (loop) 

= 

z7dc 

+ 

(loop  “ 

2 ) * j  *ws ; 

poleldc (loop) 

= 

pldc 

+ 

(loop  - 

2) *j*ws; 

pole2dc (loop) 

= 

p2dc 

+ 

(loop  “ 

2) * j  *ws ; 

pole3dc (loop) 

= 

p3dc 

+ 

(loop  - 

2) *j  *ws; 

pole4dc (loop) 

= 

p4dc 

+ 

(loop  -- 

2) *j*ws; 

poleSdc (loop) 

= 

p5dc 

+ 

(loop  - 

2) *j*ws; 

pole6dc (loop) 

= 

p6dc 

+ 

(loop  - 

2) *j*ws; 

end; 

pole7dc  =  1/Ts*log{pole7d) ; 

vectl  =  [-2  2]; 

vect2  =  [-3*ws/2  -3*ws/2]; 

vect3  =  [-ws/2  -ws/2]; 

vect4  =  [0  0] ; 

subplot (212) ;  hold  on; 

plot (vectl, vect4, ' : ' ,vect4, [-3*ws/2  3*ws/2] ,':'); 

plot (vectl, vect2 , ' ' , vectl , -vect2 

plot (vectl, vect3 , ' ' , vectl, -vect3 

plot (real (polelc) , imag (polelc) , '  +  ' )  ; 

plot (real (pole2c) , imag (pole2c) ,'+'); 
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plot (real (poleSc) , imag{pole3c) 
plot (real (pole4c) , imag(pole4c) 
for  loop  =  1:4, 

plot (real (zeroldc) , imag (zeroldc) , 'o' ) ; 
plot (real (zero2dc) , imag(zero2dc) , 'o' ) ; 
plot (real (zero3dc) , imag (zero3dc) , 'o' ) ; 
plot ( real ( zero4dc ) , imag ( zero4dc ),' o' ) ; 
plot (real (zeroSdc) , imag(zero5dc) , 'o' ) ; 
plot (real (zero6dc) , imag(zero6dc) , 'o' )  ; 
plot (real (zero7dc) , imag(zero7dc) ,  "'* ) ; 
plot (real (poleldc) , imag (poleldc) , 'x') ; 
plot (real (pole2dc) , imag(pole2dc) , 'x' ) ; 
plot (real (poleldc) , imag (poleldc) , 'x' ) ; 
plot (real (pole4dc) , imag(pole4dc) , 'x' ) ; 
plot (real (poleBdc) , imag(pole5dc) ,'x'); 
plot (real (pole6dc) , imag (poleSdc) ,'x'); 
end; 

plot (real (pole7dc) , imag(pole7dc) ,'*'); 
hold  off;  axis([-1.25  0.1  -3*ws/2  3*ws/2]); 
title ( 'S*-Plane  Root  Locus'); 

print  -deps  -tiff  figSll 


%  ***  Plot  System  Bode  Plot  for  defined  Kp  *** 

numpoints  =  2000;  %  Number  of  frequency  points  for  FRF 

num  =  K*vml''2*wn2''2  ; 

deni  =  [1  2*zetal*wnl  wnl''2]  ; 

den2  =  [1  2*zeta2*wn2  wn2''2]  ; 

den  =  conv ( deni , den2 ) ; 

Gps  =  tf (num, den);  %  Continuous  plant 

Hz  =  c2d(Gps,Ts, 'zoh') ;  %  H*(z):  Sampled  Plant  with  ZOH 

w  —  linspace ( 0 . 001 , 5 0 , numpoints ) ;  %  Linearly  spaced  frequency  vector 

normw  =  w./wd2;  %  Normalized  frequency  vector  for  plotting 

%  PID  Controller 

Gcznum  =  [Kp+Ki+Kd  -Kp-2*Kd  Kd] ; 

Gczden  =  [1  -1  0] ; 

Gcz  =  tf  (Gczniim,  Gczden)  ; 

[magGcz,phsGcz]  =  bode(Gcz,w); 
magGcz  =  reshape (magGcz, 1, numpoints ) ; 
phsGcz  =  reshape (phsGcz , 1 , numpoints) ; 

phsGcz  =  phsGcz*pi/180 ;  %  Convert  phase  from  degrees  to  radians. 

%  Continuous  System 
[magGps,phsGps]  =  bode (Gps, w); 
magGps  =  reshape (magGps , 1 , numpoints ) ; 
phsGps  =  reshape (phsGps,!, numpoints ) ; 

phsGps  =  phsGps*pi/180;  %  Convert  phase  from  degrees  to  radians. 


%  Discrete  System 

[magHz , phsHz ]  =  bode (Hz, w); 
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magHz  =  reshape (magHz , 1 , numpoints) ; 
phsHz  =  reshape (phsHz , 1 , numpoints ) ; 

phsHz  =  phsHz*pi/180 ;  %  Convert  phase  from  degrees  to  radians. 

%  Calculate  magnitude  and  phase  at  discrete  frequency  points 
for  count  =  1: numpoints 

ZOH  =  (l-exp (-j*w (count ) *Ts) )/ (j*w( count ) *Ts) ; 

Gpsl  =  magGps (count) *exp ( j  *phsGps (count) ) ; 

Gczl  =  magGcz (count )* exp (j *phsGcz (count) ) ; 

Hpz  =  magHz (count) *exp(j*phsHz (count) ) ; 
numcl  =  Gczl*ZOH*Gpsl; 
dencl  =  1  +  (Gczl*Hpz) ; 
cltf (count)  =  numcl/dencl; 

end 

magFRF  =  20*logl0 (abs (cltf ) ) ; 
phsFRF  =  angle (cltf ) *18 0/pi; 

a  =  max(magFRF); 
plotmax  =  a  +  5; 
plotmin  =  plotmax  -  160; 

figure (2);  clf; 
subplot (211)  ; 
plot (w, magFRF) ;  grid; 
phrase  =  ['Kp  =  ' , num2str (Kp) , ' ;  Ki  =  ' ,num2str (Ki) , ' ; 

Kd  =  ' ,num2str (Kd) ] ; 
title (phrase) ; 

phrase  =  ['wnl  =  ' , num2str (wnl) , ' ;  zetal  =  ' ,num2str (zetal) , ' ; 

wn2  =  ' ,num2str (wn2) , ' ;  zeta2  =  ' ,num2str (zeta2) ] ; 
xlabel (phrase) ; 
y label ( 'Magnitude  (dB) ' ) ; 
axis ( [0  50  plotmin  plotmax] ) ; 
subplot (212) ; 
plot (w,phsFRF) ;  grid; 
axis ( [0  50  -190  190] ) ; 
xlabel ( 'Frequency' ) ; 
ylabel ( ' Phase  (deg) ' ) ; 

%print  -deps  -tiff  fig68 


%  Magnitude  (in  dB)  of  FRF 
%  Phase  (in  degrees)  of  FRF 


[F,  G, H, J] =tf2ss (num, den) ; 

%  ***  Plot  Impulse  Response  for  defined  Kp  *** 

%  =  =  =  =  =  =  =  =  =  =  =  =  =:  =  =  =  =  =  =  =  =:==:=  =  =  =  =  =  =  =  =:  ===  =  =  =  =:  =  =  =  =  = 

Tn  =  Ts/istepsl;  %  Intermediate  time  steps 

[phi, gamma]  =  c2d(F,G,Tn);  %  Continuous  approximation  in  discrete  domain 

%  Define  initial  conditions 

x  =  [  0 ;  0 ;  0  ;  0  ]  ; 

tcon  =  0;  ycon  =  0; 

tdis  =  0;  ydis  =  0; 

uthen  =  0;  ethen  =  0;  ethen2  =  0; 
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%  step  through  calculations  for  "maxtime”  seconds 
for  countl  =  0 : round (max time/ Ts) ; 

%  Define  unit  impulse: 
if  countl  ==  0 
r  =  1; 
else 

r  =  0; 

end 

enow  =  r  “  ydis (countl  +  1); 

u  =  uthen  +  Kl*enow  +  K2*ethen  +  Kd*ethen2;  %  Control  effort 

%  Calculate  continuous  output 
for  count2  =  liistepsl; 

tcon  =  [tcon; (countl*istepsl  +  count2)*Tn]; 
x  =  phi*x  +  gamma*u; 
ycon  =  [ycon;real (H*x  +  J*u) ] ; 
end 

index  =  istepsl* (countl  +  1)  +  1; 
tdis  =  [tdis; tcon (index) ] ; 
ydis  =  [ydis ;ycon( index) ] ; 

%  Age  data 
uthen  =  u; 
ethen2  =  ethen; 
ethen  =  enow; 

end 

a  =  max  (ycon)  ; 
b  =  min  (ycon)  ; 
plotmax  =  a  +  0.2*a; 
plotmin  =  b  +  0.2*b; 

%  Plot  Results 
figure (3);  clf; 
subplot (211)  ; 

plot ( tcon . /Ts , ycon, ' - ' , tdis . /Ts , ydis r— tdis . /Ts , ydis r* ') ;  grid; 
%plot { tdis . /Ts , ydis ; 

phrase  =  ['wnl  =  ' ,nuin2str (wnl) , ' ;  zetal  =  ' ,num2str (zetal) , ' ; 

^^2  =  ' , nuin2str  (wn2 )  ,  '  ;  zeta2  =  '  ,nuin2str  (zeta2)  ]  ; 
xlabel (phrase) ; 

axis([0  maxtime/Ts  plotmin  plotmax]); 

phrase  =  ['Impulse  Response,  FR  =  ' ,num2str (FR) ] ; 

title (phrase) ; 


%  ***  Plot  Step  Response  for  defined  Kp  *** 

%  Define  initial  conditions 

X  =  [0;0;0;0] ; 

tcon  =  0;  ycon  =  0; 

tdis  =  0;  ydis  =  0; 

uthen  =  0;  ethen  =  0;  ethen2  =  0; 

^  =  1;  %  Unit  step 
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%  step  through  calculations  for  "maxtime"  seconds 
for  countl  =  0 : round (maxt ime/Ts) ; 
enow  =  r  -  ydis (countl  +  1) ; 

u  =  uthen  +  Kl*enow  +  K2*ethen  +  Kd*ethen2;  %  Control  effort 

%  Calculate  continuous  output 
for  count2  =  l:istepsl; 

tcon  =  [tcon; (countl*istepsl  +  count2)*Tn]; 

X  =  phi*x  +  gamma*u; 

ycon  =  [ycon;real (H*x  +  J*u) ] ; 

end 

index  =  istepsl* (countl  +  1)  +  1; 
tdis  =  [tdis; tcon (index) ] ; 
ydis  =  [ydis ; ycon (index) ] ; 

%  Age  data 
uthen  =  u; 
ethen2  =  ethen; 
ethen  =  enow; 

end 

a  =  max  (ycon)  ; 
b  =  min  (ycon)  ; 
plotmax  =  a  +  0.2*a; 
plotmin  =  b  +  0.2*b; 

%  Plot  Results 
subplot (212) ; 

plot (tcon. /Ts ,ycon, ' - ' , tdis . /Ts ,ydis, 'r — ' , tdis . /Ts,ydis, 'r* ' ) ;  grid; 
%plot (tdis . /Ts ,ydis ,'-'); 

xlabel ( 'Normalized  Time  (t/Ts) ' ) ;  ylabel ( 'Response' ) ; 
axis  ( [0  maxt  ime/Ts  plotmin  plotmax] )  ; 

phrase  =  ['Step  Response;  Kp  =  ' , num2str (Kp) , ' ;  Ki  =  ' ,num2str (Ki) , ' ; 

Kd  =  ' ,num2str (Kd) ] ; 
title (phrase) ; 

print  -deps  -tiff  fig612 
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